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1.0  INTRODUCTION  AND  SUMMARY 

A  two-degree-of-freedom  (2-DOF)  digital  computer  program  called  ETO  (Earth-to-Orbit) 
has  been  written.  Five  equations  of  notion  describing  weight,  velocity,  flight 
path  angle,  altitude,  and  range  of  a  vehicle  are  integrated  by  a  second-order 
routine.  The  program  assumes  a  spherical,  nonrotating  earth.  The  primary  purpose 
of  the  program  is  to  evaluate  horizontal  takeoff  earth-to-orbit  vehicles  of  single- 
or  two-stage  configurations.  Options  are  available  in  the  program  to  evaluate 
cruise  vehicles  of  various  types,  aircraft,  and  ground  or  air  launched  missiles. 

Decause  the  vehicle  motion  is  in  a  single  plane,  i.e.,  2-DOF,  simplified 
aerodynamics  are  used  in  the  form  of  user-input  coefficients  for  zero-lift  drag,  the 
lift  curve  slope,  drag  due  to  lift  factor,  and  friction  drag.  Airbreathing 
propulsion  performance  is  user  input  as  specific  impulse  vs  fuel  equivalence  ratio 
and  flight  Mach  number.  These  data  along  with  input  values  of  inlet  air  capture 
area  are  sufficient  to  define  airbreathing  thrust  at  every  integration  point. 

Rocket  thrust  is  user  input  as  specific  impulse  vs  altitude  and  propellant  flow  rate 
vs  Mach  number.  Airbreathing  and  rocket  thrust  may  be  used  separately  or 
simultaneously  in  any  phase  of  the  flight. 

A  number  of  options  are  available  for  tailoring  earth-to-orbit  missions  or  cruise 
missions.  For  earth-to-orbit  missions,  options  exist  to  climb  to  orbit  at  a 
commanded  constant  dynamic  pressure,  along  a  varying  dynamic  pressure  path,  or  a 
ballistic  ascent  beginning  at  a  prescribed  velocity.  For  cruise  aircraft,  constant 
altitude  cruise  or  maximum  L/D  cruise  may  be  specified. 

For  ground-launched  missiles,  launch  at  a  prescribed  initial  angle  is  used.  For  air 
iaunches,  an  initial  velocity  and  angle  are  prescribed.  Missiles  may  be  single 
s^age  or  two  stage. 

A  major  objective  of  the  program  v/as  to  achieve  reasonable  accuracy  and  very 
short  run  times  on  microcomputers .  Run  times  of  1  to  3  minutes  were  sought. 
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2.0  DEVELOPMENT  AND  DESCRIPTION  OF  PROGRAM  ETO 


2.1  Program  Language  and  Compilers 

★ 

Initially,  the  program  was  written  in  Microsoft  GW-BASIC  with  strict  adherence  to 
syntax  that  would  permit,  it  to  run  under  the  BASICA  interpreter  for  debugging.  Once 
a  suitable  version  was  obtained,  it  was  compiled  with  the  Microsoft  BASIC  Compiler, 
BASCOM.  A  little  later  in  the  program  development,  the  Microway  87  BASIC  compiler 
was  obtained.  It  supported  the  Intel  8087/80287  floating  point  math  chip, and 
program  execution  speed  was  improved  greatly.  Near  the  end  of  the  code  development 
cycle  fa  refinement  phase),  the  Microsoft  Quick  BASIC  4.0  compiler  and  the  Borland 
Turbo  BASIC  1.0  compilers  became  available.  They  both  support  the  8087/80287  math 
chip,  do  not  require  line  numbers  in  the  codes,  and  support  more  structured 
programming.  The  code  as  listed  in  this  report  now  conforms  to  the  syntax 
requirements  of  the  Microsoft  Quick  BASIC  4.0  compiler.  It  will  not  run  under  the 
BASIC  interpreter.  The  integrated  editing,  debugging,  and  run  environments  of  new 
Microsoft  or  Borland  compilers  greatly  ease  code  development.  They  are  available  at 
prices  as  low  as  $65. 

The  source  code  is  about  45,000  bytes,  and  a  stand-alone  executable  file  of  the  code 
is  about  85,000  bytes.  On  -XT  or  -AT  compatible  machines  with  the  8087/80287 
floating  point  math  chip,  very  short  run  tines  are  possible.  The  following  gives 
some  results  obtained  for  the  example  problem  TEST1.DAT,  an  earth-to-orbit  flight: 


Computer 

Run 

Time,  secs 

Zenith  Z248  (AT  class) 

8  MHz  80286/5  MHz  80287 

38 

secs 

8  MHz  80286/no  math  chip 

Sanyo  550  (Turbo  XT  class) 

210 

secs 

7.2  MHz  V-20  /8  MHz  8087 

60 

secs 

ACER  1100 

16  MHz  80386/16  MHz  80387 

Vendex  Turbo  XT 

12 

secs 

8  MHz  8088/no  math  chip 

916 

secs 

*  At  time  of  publication,  a  Fortran  version  is  now  available.  Contact 
authors  for  detail s  . 
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2.2  Equations  of  Motion  in  a  Plane 


Figure  2-1*  is  a  free  body  diagram  of  a  flight,  vehicle  in  notion  about  a  spherical, 
nonrotating  earth.  The  vehicle,  of  mass  W/gc ,  is  at  a  distance  (R  +  hQ)  from  the 
center  of  the  earth  of  mass  Mg.  The  vehicle  is  moving  along  a  path  at  angle  y 
relative  the  earth.  Vehicle  thrust,  T,  acts  along  the  axis  of  the  vehicle.  The 
vehicle  is  inclined  to  its  axis  of  motion  by  its  angle  of  attack,  a.  Lift,  L,  is 
exerted  in  a  direction  normal  to  the  axis  of  notion,  and  drag,  D,  is  exerted  along 
the  axis  of  motion.  The  weight,  U,  and  the  centrifugal  force,  F  ,  act  oppositely 
and  along  the  radius  to  the  earth. 

The  motion  of  the  vehicle  is  described  by  its  acceleration.  It  is  most  convenient 
to  resolve  the  acceleration  along  the  motion  axis  (the  velocity  axis)  and  along  a 
direction  normal  to  the  motion  axis.  The  acceleration  along  the  motion  axis  will  be 
defined  from  Newton's  law  for  acceleration: 


Resolving  all  of  the  forces  to  their  components  along  the  notion  axis  gives 

W  dv  _  _  „  q 

g  dt  ~  T  cos  a  -  D  -  W  J  sin  y  +  F  sin  y 

c  yc  c 

The  acceleration  normal  to  the  axis  of  motion  will  be  given  by  Newton's  law  for 
accel eration : 


W 

9 


c 


dt 


E  F 


n 


The  rate  of  change  of  flight  path  angle  is 

dY  .  dvn  dY 
dt  dt  ’  dv 


♦Section  2  figures  begin  on  page  11. 
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Since 


dY  _  1 


He  have 


jv  dv 
dY  _  _ n 

dt  ”  dt 


1 


v 


o 


and 


So  that  Mewton's  law  gives 


dv 
_ n 

dt 


dY 

dt 


W_  v  dy  = 
gc  o  dt 


Z  F 


n 


The  summation  of  normal  forces  gives 

w  dY  ,  T  .  ,,  a 

vo  dt  L  +  Tsma-W|—  cos  y  +  Fc  cos  y 

The  above  two  expressions,  one  for  the  axial  acceleration  and  one  for  the  normal 
acceleration,  are  sufficient  to  describe  the  time-dependent  motion  of  the  vehicle 
but  are  not  yet  sufficient  to  allow  solution. 

An  expression  for  the  time  rate  of  change  of  weight  is  needed,  and  since  drag  and 
lift  are  also  dependent  on  altitude,  an  expression  for  the  time  rate  of  change  of 
altitude  is  needed.  The  time  rate  of  change  of  range  is  not  required  but  is  of 
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interest  for  cruising  flight.  The  final  system  of  five  equations  to  be  integrated 
becomes 


Weight  change: 


dW 

dt 


=  -  W 


'elocity  change: 


dvo .  fc  r. 

dt  ‘ 


T  cos  a  -  D  -  W 


2_ 

9. 


sin  a 


+  F_  sin  Y 


] 


Flight  path  angle  change: 


Altitude  change: 


dY 

dt 


L  +  T  sin  ot  +  Fc  cos  Y  -  W  cos  Y 
o  ^c 


dh 

dt 


vQ  sin  y 


Range  increment: 

dR 

dt  =  vo  cos  Y 

The  first  four  equations  are  coupled  and  must  he  integrated  by  employing  an 
algorithm  for  integrating  first-order  systems  of  differential  equations.  As  was 
stated  above,  the  fifth  equation,  which  tracks  the  range,  is  of  interest  for 
cruising  flights.  We  decided  to  employ  a  simple  two-step  integration  scheme  (second- 
order  Runge-Kutta,  sometimes  called  Heun's  method),  and  we  hoped  that  its  accuracy 
would  be  compatible  with  the  choice  of  large  time-integration  steps,  i.e.,  2  to  5 
seconds,  and  the  desire  to  keep  program  execution  time  short. 
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The  overall  strategy  for  solving  the  problem  of  a  flight  trajectory  is  to  integrate 
the  coupled  equations  of  motion  forward  in  time,  enforcing  constraints/limits  on 
values  of  the  variables  until  some  predefined  value  of  one  of  the  variables  is 
reached  and  triggers  a  stop  to  the  problem.  The  solution  to  the  trajectory  is  the 
time-ordered  history  of  values  of  all  of  the  relevant  variables  up  to  the  stopping 
condition . 

At  each  integration  time  step,  current  values  are  needed  for 

W  =  rate  of  weight  change  (usually  propellant  flow) 

W  =  weight 
T  =  thrust 
D  =  drag 
L  =  lift 

Fr  =  centrifugal  force 
a  =  vehicle  angle  of  attack 
Y  =  vehicle  flight  path  angle 
v0  =  flight  velocity 

g  =  acceleration  due  to  gravity  (varies  with  altitude) 

Constraints  are  applied  by  defining  limits  for  such  variables  as  a  and  y.  Thrust  ma 
be  limited  by  commanded  axial  acceleration  or  bv  maximum  fuel  flow.  Other 
constraints  such  as  maximum  and  minimum  flight  dynamic  pressure  impose  adjustments 
on  a  and  thrust.  Basically,  a  trajectory  can  be  flown  by  using  two  commands, 
angle-of-attack  and  axial  acceleration  ("stick  and  throttle").  The  trajectory 
proceeds  through  a  series  of  flight  phases  (up  to  five  are  used  in  ETO)  which 
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enforce  the  problem  constraints.  The  main  program  module  enforces  the  constraints 
appropriate  to  each  flight  phase.  At  the  end  of  each  integration  step,  the  program 
returns  to  the  beginning  of  this  main  module  and  "falls  through"  to  the  current 
flight  phase  to  obtain  the  appropriate  angl e-of-attack  and  acceleration  commands. 

The  flight  proceeds  through  phases  based  on  computed  events  and  on  inputs  from  the 
vehicle  input  data  file. 

The  flow  diagram  of  Figure  2-2  shows  the  basic  program  strategy  in  which  lift,  drag, 
and  the  flight  conditions  ChQ,  vQ)  are  defined;  then  the  main  module  is  entered  to 
obtain  angl e-of-attack  and  acceleration  commands,  followed  by  generation  of  the 
current  thrust  value  and  finally  going  to  the  integration  module.  A  runtime  output 
step  is  taken,  and  the  program  returns  to  update  lift,  drag,  and  the  flight 
condition  (hQ,  vQ).  This  process  is  repeated  until  the  problem  is  completed. 

The  remaining  parts  of  Section  2  give  more  specifics  on  how  aerodynamics, 
propulsion,  flight  phases,  and  stopping  conditions  are  handled. 

2.3  Aerodynamics 

Because  the  program  is  two-degree-of-freedom,  i.e.,  motion  in  a  plane,  only  trim 
aerodynamics  are  used.  The  term  "trim  aerodynamics"  means  use  of  lift  and  drag 
coefficient  values  at  which  the  sum  of  all  moments  about  the  vehicle  are  zero.  A 
data  table  of  Cnn  vs  Mach  number,  a  data  table  of  C.  and  K  vs  Mach  number  and  a 
table  of  friction  increment  to  C^q  vs  altitude  defines  the  aerodynamics. 

At  a  given  Mach  and  altitude,  the  data  tables  are  consulted  for  C^q,  C^ft,  K  and  ACp 
then  used  in  the  following  equations: 


L  =  CL  qQ  A  Ref 

CD  =  CD0  +  K  CL  +  ACD 
0  =  CD  qo  A  Ref 
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2.4  Propulsion 


Airbreathing  propulsion  is  described  by  tables  of  Isp  vs  Mach  number  and  equivalence 
ratio  (<f>),  and  a  schedule  of  maximum  vs  Mach  number.  This  information,  along  with 
tables  of  inlet  capture  area  ratio  {ho/br)  vs  Mach  and  angle  of  attack,  the  inlet 
cowl  area  f Ac ) ,  and  stoichiometric  fuel -air  ratio  are  sufficient  to  establish 
thrust,  fuel  flow,  and  airflow  at  each  integration  time  step.  At  each  time  step, 
the  following  steps  are  taken  to  obtain  thrust  and  fuel  flow: 


Ao 


p  v  A 
ooo 


at  a  specified  <j> ; 
f/a  =  4>  (f/a)  STOICH 
Wf  ■  (f/a) 

T  =  Isp 


An  interation  scheme  is  used  to  adjust  by  testing  if  T  is  greater  or  less  than 
required  to  produce  the  commanded  acceleration. 

Rocket  propulsion  is  described  by  tables  of  Isp  vs  altitude  and  by  propellant  flow 

rate  (W  )  vs  Mach  number: 

P 

T  *  Wp  Isp 

Rocket  thrust  is  nodulated  in  response  to  the  acceleration  command.  If  the 
airbreather  and  the  rocket  are  operating  simultaneously,  the  rocket  is  throttled 
first  to  reduce  thrust  to  the  commanded  value. 

2.5  FI ight  Phases 

Figure  2-3  shows  a  typical  succession  of  flight  phases  as  discussed  in  the  following 
paragraphs. 
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Phase  I  is  the  takeoff  roll.  The  vehicle  accelerates  to  the  takeoff/rotation 
velocity  at  zero  lift  and  zero  commanded  climb.  Once  the  rotation  velocity  is 
achieved,  control  passes  to  Phase  II. 

Phase  II  is  a  climb  at  a  commanded  pull-up  load  factor  and  axial  acceleration. 

After  the  rotation  velocity  is  exceeded  in  Phase  I,  the  vehicle  is  permitted  to 
pitch  up  to  a  maximum  angle  of  attack  or  a  commanded  load  factor,  whichever  limit  is 
controlling.  The  vehicle  climbs  and  accelerates  to  a  commanded  maximum  subsonic 
velocity  whi^h  cannot  be  exceeded  until  a  specified  altitude  is  reached.  Once  this 
llach  and  altitude  threshold  is  reached,  angle  of  attack  is  commanded  to  zero,  and 
the  vehicle  accelerates  on  a  ballistic  path  (no  induced  drag)  through  the  transonic 
region  until  a  commanded  dynamic  pressure  is  reached.  This  completes  Phase  II. 

Phase  III  is  a  climb  and  acceleration  at  a  commanded  dynamic  pressure,  q  com.  This 
value  of  q  nay  be  constant  or  may  be  made  linearly  variable  with  velocity,  beginning 
at  a  specified  velocity.  A  maximum  and  a  minimum  value  of  q  is  specified.  A 
control  algorithm  generates  angl e-of-attack  and  acceleration  commands  at  each 
integration  step  to  keep  the  vehicle  within  the  minimum  and  maximum  q  limits.  The 
acceleration  command  is  generated  by  examining  whether  q  is  greater  or  less  than  the 
commanded  value  and  whether,  in  either  case,  it  is  getting  larger  or  smaller.  The 
algorithm  for  implementing  this  control  is  shown  in  Figure  2-4.  A  more  detailed 
discussion  is  in  Section  3.4.  The  output  of  this  algorithm  is  a  new  value  of  a, 
angle  of  attack,  to  be  used  in  the  next  integration  step.  With  proper  choices  for 
the  adjustment  coefficients,  a  relatively  smooth  ascent  path  can  be  achieved.  Once 
the  a  command  is  generated,  an  acceleration  command  is  generated.  The  algorithm  for 
this  can  be  visualized  from  Figure  2-5.  Basically,  the  strategy  is  that  if  q  is 
greater  than  the  maximum,  then  no  acceleration  is  permitted.  If  q  is  less  than  the 
connanded  value,  then  the  maximum  available  acceleration  is  commanded.  In  the 
region  where  q  is  greater  than  commanded  but  less  than  maximum,  a  proportional 
reduction  in  acceleration  is  commanded  as  suggested  by  Figure  2-5.  Acceleration 
along  the  commanded  q  path  continues  until  a  cruise  command  or  a  ballistic  path 
command  is  made. 

Phase  IV  generates  a  commanded  zero  acceleration  to  maintain  a  specified  velocity  at 
a  commanded  q  (altitude)  or  to  maintain  velocity  at  maximum  L/D .  Maximum  L/D  cruise 
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is  accomplished  by  finding  the  angle  of  attack  that  produces  maximum  L/D  and  then 
maintaining  it  to  fuel  exhaustion. 

Phase  V  is  a  ballistic  ascent  and  is  normally  commanded  after  having  proceeded  along 
a  Phase  III  commanded  q  path.  It  is  used  primarily  for  earth-to-orbit  flights.  At 
a  commanded  velocity,  the  vehicle  pulls  up  to  a  specified  angle  of  attack, 
maintaining  it  until  a  specified  flight  path  angle  is  achieved.  Once  the  flight 
path  angle  is  reached,  the  angle  of  attack  is  commanded  to  zero  and  the  flight  is 
ballistic,  i.e.,  angle  of  attack  equals  zero  from  this  point  on.  During  Phase  V  an 
acceleration  command  different  from  that  applied  in  Phase  III  may  be  used. 

2.6  Stopping  Conditions 

Three  stopping  conditions  are  utilized:  fuel  exhaustion,  achievement  of  commanded 
velocity,  and  achievement  of  the  orbital  condition.  When  the  vehicle  weight  falls 
to  the  specified  final  value,  the  program  stops  and  issues  a  message  to  that  effect. 
For  two-stage  vehicles,  the  first  stage  is  usually  flown  to  a  prescribed  velocity  or 
to  fuel  exhaustion.  When  either  occurs,  the  program  halts,  issues  a  message  that 
the  first-stage  flight  has  been  completed  and  asks  if  you  care  to  continue  the 
problem  with  the  second  stage. 

The  orbital  condition  is  defined  by  the  centrifugal  force  becoming  equal  to  the 
weight.  When  this  occurs  the  program  stops  and  advises  that  orbit  has  been 
achieved.  Because  the  code  is  assuming  a  nonrotating  earth,  the  orbital  velocity 
required  is  always  equivalent  to  be  that  of  a  polar  orbit.  Normally  this  will  be 
slightly  under  25,900  fps.  The  exact  value  depends  on  altitude.  An  Eastern 
equatorial  orbit  would  be  the  minimum  required  orbital  velocity  and  would  be  on  the 
order  of  24,500  fps.  The  user  could  specify  a  final  velocity  less  than  the  polar 
orbit  value  that  is  appropriate  to  the  orbit  inclination  of  interest  to  him. 


10 


12 


14 


Figure  2-5.  Acceleration  Command 


3.0  RUNNING  PROGRAM  ETO 


3.1  Hardware  Requirements 

The  distribution  disk  is  a  360K  I BM/MS-DOS  format  5  *s-inch  floppy  disk  which 
contains  the  following  files: 

ET013.BAS  (See  Note  Below) 

ET013.BAC 

ET013.EXE 

TEST1.DAT 

STG1.DAT 

STG2.DAT 

STG2BAL.DAT 

ETO. TAB 

GUID14.BAS 

ET0123.BAS 

ET0123.EXE 

ET013.BAS  is  the  program  source  code.  ET013.3AC  is  a  backup  copy  of  the  source  code 
and  ET013.EXE  is  a  stand-alone  executable  version  of  the  program.  The  designation 
ET013  reflects  the  version  number  of  program  ETO  available  at  the  time  of  this 
report.  Significant  upgrades  will  have  higher  numbers.  The  files  with  a  .DAT 
suffix  are  sample  input  data  tiles.  Each  data  file  illustrates  a  significant 
feature  of  the  program.  Each  of  these  .DAT  files  is  discussed  in  Section  3.6.  The  file 
ETO. TAB  is  an  alphabetical  listing  with  description  of  all  the  program  variables. 

The  file  GUID14.BAS  is  an  optional  angl e-of-attack  module  for  flight  Phase  III  which 
the  user  may  substitute  for  GUID13  in  the  trajectory  program  (Section  2.5). 

ET0123.BAS  and  ET0123.EXE  are  the  source  code  and  executable  file,  respectively, 
which  may  be  used  to  generate  a  file  readable  by  LOTUS  123  of  the  program  output 
file.  This  is  included  for  users  who  may  wish  to  use  LOTUS  123  to  generate  plots  of 
the  program  output. 

NOTE:  Since  completion  of  this  report  FORTRAN  versions  are  also  available.  Contact 
the  authors  for  more  information. 
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The  mininum  hardware  equipment  for  running  the  stand-alone  .EXE  version  of  ETO  is  an 
IBM  or  compatible  XT  or  AT  microcomputer  using  MS .DOS  2.1  or  higher.  DOS  versions 
earlier  than  2.1  will  not  work  properly.  One  hundred  percent  IBM  compatibility  does 
not  appear  to  be  required,  however.  A  single  disk  drive  (360K,  5^-inch  floppy  disk) 
and  256K  RAM  are  needed.  It  may  be  possible  to  run  the  program  with  less  than  256K 
RAM  but  this  has  not  been  demonstrated.  A  math  coprocessor,  8087/80287,  is  highly 
desirable  for  reducing  run  times  but  is  not  required.  A  printer  is  required  for 
printed  output  of  the  time  history  of  a  flight.  An  Epson  FX  or  compatible  dot 
matrix  printer  is  the  best  choice.  With  this  printer  letter-width  paper  can  be  used 
because  this  printer  accepts  the  program  commands  to  set  compressed  print.  If  your 
letter-width  printer  cannot  respond  to  the  program  command  for  compressed  print  and 
cannot  be  set  by  a  hardware  switch,  you  will  need  to  employ  a  wide-carriage  printer. 

The  most  convenient  environment  for  running  the  program  is  from  within  the  Microsoft 
Quick  BASIC  4.0  editor/compiler.  The  computer  will  have  to  be  100  percent  IBM 
compatible.  The  program  was  developed  using  this  compiler,  and  the  stand-alone  .EXE 
version  was  generated  by  it.  MS-DOS  versions  2.11,  3.1,  and  3.2  were  used  in  the 
development.  It  does  require  more  compute:  hardware  resources,  however.  Microsoft 
advises  that  320K  RAM  is  needed.  The  mininum  RAM  used  during  development  was  512K. 
It  was  satisfactory,  but  it  is  not  known  how  close  this  is  to  the  minimum  needed  for 
this  program.  Two  360K  floppy  disk  drives  a.o  needed.  A  hard  drive  is  much  more 
convenient  and  permits  quicker  compilation.  A  RAM  drive  is  also  convenient.  Quick 
BASIC  4.0  senses  if  a  math  coprocessor  [8087/80287)  is  present  and  uses  it  if 
available.  Version  4.0  of  Quick  BASIC  is  distinctly  different  from  earlier 
versions.  These  earlier  versions  are  not  recommended. 

Some  limited  development  work  was  done  using  the  Borland  Turbo  BASIC  1.0 
editor/compiler.  It  does  not  require  100  percent  IBM  compatibility.  It  provides  an 
integrated  program  development  and  compiler  environment.  Its  choice  over  the 
Microsoft  product  is  a  matter  of  personal  preference.  To  use  the  Borland  compiler, 
the  program  variable  AT  would  have  to  be  changed  to  some  other  name.  AT  is  a 
reserved  work  in  Turbo  BASIC.  Some  work  would  be  needed  in  formatting  the  screen 
and  printer  output.  Either  the  Microsoft  or  the  Borland  compiler  will  be  required 
if  the  user  desires  to  modify  the  code  and  recompile  it. 
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3.?.  Using  the  Stand-Alone  .EXE  Version  of  ETO 

The  first  step  should  be  to  make  a  backup  copy  of  the  distribution  disk.  The  backup 
copy  should  be  used  for  running  the  program.  We  recommended  that  you  copy  your 
version  of  the  MS-DOS  COMMAND.COM  file  to  the  backup  disk.  The  reason  for  this  is 
that  the  program  uses  the  BASIC  "SHELL"  command  to  temporarily  exit  to  DOS  to  copy 
the  input  data  file  to  the  printer  when  the  program  prints  out  results.  The  program 
needs  access  to  COMMAND.COM  on  the  default  drive  in  order  for  the  "SHELL"  command  to 
work  properly.  MS-DOS  2.1  or  greater  is  required. 

Once  the  above  has  been  accomplished,  the  program  is  ready  to  run.  At  the  DOS 
prompt,  type  ET013  and  press  RETURN.  The  screen  clears,  and  the  following  will  be 
displ ayed: 

IS  THIS  A  SINGLE-STAGE  (1),  OR  A  TWO-STAGE  VEHICLE  (2)’ 

Respond  by  entering  a  1  or  a  2  and  pressing  the  ENTER  key.  If  you  responded  to  the 
above  with  a  1,  you  will  be  prompted  with: 

SINGLE-STAGE  FILE  NAME? 


Enter  the  name  of  the  input  file  following  the  DOS  file  name  convention,  DRIVE: 
NAME. EXT.  The  input  file  does  not  need  to  be  on  the  default  drive. 

If  you  had  responded  with  a  2  to  the  first  question,  you  would  be  asked  the 
following: 

FIRST-STAGE  FILE  NAME? 

Again  enter  input  file  name  as  DRIVE:  NAME. EXT  and  press  RETURN  key. 

You  will  then  be  prompted  with: 

SECOND-STAGE  FILE  NAME? 

Respond  as  above. 
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Next  you  will  be  asked: 


"Do  you  want  a  listing  of  tabular  input  data  (Y/N)?" 

Normally  you  will  not  want  this  and  will  respond  with  a  "no"  ( N ) . 

If  you  respond  with  "yes"  (Y),  then  the  nine  data  tables  in  the  first/single-stage 
input  file  are  printed  out.  The  printout  is  convenient  for  examining  the  data 
tables  for  accuracy. 

The  next  prompt  will  be: 

"Do  you  want  printer  output  during  run  (Y/N)?" 

As  above,  you  will  usually  not  want  this.  The  printer  output  format  is  not 
particularly  convenient  and  the  program  runs  much  more  slowly  in  this  mode.  The 
principal  use  for  this  printout  is  diagnosis  of  problems  early  in  the  flight.  The 
program  can  be  stopped  after  the  desired  flight  tine  has  passed.  The  printer  output 
interval  is  the  same  as  the  screen  print  interval,  DELPRINT. 

Whether  you  respond  with  "Y"  or  "N"  the  next  screen  message  displays  the  value  of 
the  last  item  in  the  input  file.  If  the  screen  shows  a  value  different  from  what 
you  intended,  then  the  input  file  is  in  error.  Also  you  will  be  prompted  with  the 
fol 1  owing : 

"Do  you  want  a  pi  otter/printer  file  created  (Y/N)?" 

If  you  respond  with  "Y",  then  you  will  be  prompted  with: 

"Enter  the  pi  otter/printer  file  name  (DRIVE:  NAME. EXT)". 
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After  entering  the  desired  drive  and  file  nane,  press  RETURN  and  the  following  will 
appear: 


At  the  proipt  belos, enter  a  value  for  FREQ , the  interval  of  plotter/ 
printer  output. 

A  value  of  1  for  FFEQ  leans  output  at  each  tiie  increient 
A  value  of  2  leans  every  other  tiie  increient, etc.  As  a  guide  to 
the  aiount  of  disk  space  required, use  the  following : 

BTTES  REQUIRED  :  (200)  I  (FLIGHT  TIME) 


(IHTEGRATIOS  I ETERI AL)  I  (FREQ) 

Very  large  files  exceeding  200E  lay  be  generated. Exit  the 
prograi  by  pressing  CTRL  C  if  disk  space  is  insufficient, 
otherwise  enter  value  for  FREQ  below: 


In  response,  enter  a  value  that  is  some  integer  multiple  of  the  integration 
interval,  e.g.,  a  FREQ  value  of  20  means  that  output  will  be  written  to  the  disk 
every  20  integration  intervals.  If  the  integration  interval  is  5  seconds,  then  the 
disk  file  would  contain  the  flight  history  in  100-second  increments.  After  entering 
the  desired  value  for  FREQ  and  pressing  ENTER,  the  screen  clears  and  the  problem 
solution  begins.  As  the  solution  proceeds,  the  screen  will  display  the  following 
type  data: 


T-2000.0  PRASE  3 , CLIMB  AT  C08SAEDED  Q umtttttt 


V0-15628 . 3  B0-15 

.1  ALPHA:  1.99 

QDOTF:  20,05 

BES1: 

0.0 

H0=  129013 . 5 

GAH8A:  0.10 

KDOTA:  431.83 

EDOTP: 

0.0 

R:  2304.32 

Q0: 1045. 90 

PHI:  1.590 

ISP:  1357.00 

ALPHAHAX:  5.00 

FT:  71861,9 

L/ti  : 

0.544 

ISPEFF-  958.7 

ALPHAUIE--  -5.00 

FC--  25796.4 

ISPHEAK= 1431 . 8 

ACCCOH  :  0.27 

LIFT:  39107.9 

L/D  : 

4.967 

1SPAIG:  2842.7 

ACCEL-G'S:  0.27 

DRAG:  7873.4 

1SPA/B--  1357.0 

T8BBST:  27209.2 

T/D  : 

3.456 

1SP80C:  0.0 

THROSTCOH:  27209.2 

THROSTSAI:  42899.2 

T/K  : 

0.379 
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The  screen  is  updated  at  the  frequency  DELPRINT,  a  user  input  to  the  input  data 
file.  A  DELPRINT  value  of  20  to  40  v/i  1 1  insure  that  the  speed  of  the  solution  is 
not  limited  by  writing  to  the  screen.  The  user  should  experiment  with  values  of 
DELPRINT  to  find  the  values  that  are  frequent  enough  to  track  the  calculation 
progress  without  impeding  the  solution  speed.  The  pi  otter /printer  file  that  is 
generated  has  values  of  the  following  25  variables  recorded  every  FREQ  X  DT  seconds 
of  the  fl ight: 


able  # 

Variabl  e 

Units 

1 

TIMEX 

secs 

2 

V0 

ft/sec 

3 

H0 

feet 

4 

RANGE/6076 

nautical  miles 

5 

ISP 

secs 

6 

ISPEFF 

secs 

7 

ISPHEAN 

secs 

8 

WDOTF 

Ibs/sec 

q 

WDOTA 

lbs/sec 

10 

PHI 

fraction  of 

stoichiometric 

11 

THRUST 

lbs 

12 

THRUSTCOM 

lbs 

13 

THRUSTMAX 

lbs 

14 

00 

lbs/ft2 

15 

ALPHA 

degrees 

16 

GAMMA 

degrees 

17 

WEIGHT 

lbs 

18 

LIFT 

lbs 

19 

DRAG 

lbs 

20 

WDOTP 

Ibs/sec 

21 

RES1 

not  used 

22 

ISPAVG 

secs 

23 

ACC 

g's,  accel/GA 

24 

ACCCOM 

g's 

25 

THRUST/DRAG 

lbs/lbs 
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One  use  of  this  file  is  for  generating  plots  of  different  variables  after  run 
completion.  The  spreadsheet  program  LOTUS  1-2-3  can  be  used  to  generate  graphics 
plots.  The  program  ET0123  on  the  distribution  disk  can  be  used  to  build  a  file  from 
this  data  that  can  be  read  by  LOTUS. 

At  the  end  of  the  flight,  you  will  be  prompted  for  a  printout  of  the  plotter/printer 
file  (presuming  you  responded  with  "Y"  at  start  of  problem).  If  you  respond  with 
"Y",  then  you  will  get  an  output  giving  the  final  run  conditions,  a  copy  of  the 
input  data  files  (first-stage  and  second-stage)  and  a  time-ordered  history  of  all  of 
the  above  variables.  This  is  the  only  opportunity  you  will  have  to  generate  this 
information  from  within  the  program.  See  the  sample  output  at  the  end  of  this 
section.  The  program  was  written  so  that  all  printed  output  could  fit  on  standard 
8'2-inch-wide  sprocket-feed  computer  paper  if  compressed  print  (17  characters  per 
inch)  is  supported  by  the  printer.  The  code  to  turn  on  and  turn  off  the  compressed 
print  is  compatible  with  the  Epson  FX  series  dot  matrix  printers.  If  your  printer 
cannot  generate  the  17-character/inch  compressed  print,  you  will  have  to  use  a 
standard  wide  carriage  printer  and  wide  paper.  If  your  printer  will  not  respond  to 
the  compressed  print  commands  generated  by  the  program  but  can  be  set  by  a  printer 
switch  to  compressed  print,  you  will  be  able  to  use  8?j-inch  paper. 

If  you  are  working  with  a  two-stage  vehicle,  you  will  be  prompted  at  completion  of 
its  first-stage  flight  that  the  second-stage  input  file  has  been  read  and  the  value 
that  was  read  from  the  input  file.  If  the  file  was  read  correctly  and  the  stopping 
conditions  of  the  first-stage  were  satisfactory,  then  the  second-stage  run  can 
proceed.  If  you  want  to  stop,  press  CTRL-C  and  the  program  will  exit  to  DOS. 
Presuming  you  wish  to  proceed,  you  will  need  to  respond  to  the  prompt  for  a  listing 
of  the  second-stage  tabular  data  (nine  input  data  tables).  lf*the  response  is  "Y" 
the  tables  will  be  listed  and  the  program  will  then  run  to  completion.  An  "N" 
response  bypasses  this  printout  and  starts  the  second-stage  problem.  At  the  end  of 
a  second-stage  run,  you  will  be  prompted  for  a  printout  of  the  plotter/printer  file 
if  you  had  requested  one  at  the  start  of  the  problem. 

Any  time  during  the  run,  pressing  CTRL-S  or  the  BREAK  key  will  suspend  program 
execution.  The  program  can  be  resumed  from  this  point  by  pressing  the  spacebar.  If 
you  wish  to  terminate  the  program  after  suspending  it,  press  the  spacebar  first  to 
resume  it,  then  press  CTRL-C  to  terminate.  The  program  may  be  terminated  at  any 
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tine  by  pressing  CTRL-C.  If  you  had  requested  a  plotter/printer  file  at  the 
beginning  of  the  run. you  will  be  prompted  for  a  printout  of  the  run  up  to  the  point 
of  termination. 

3.3  Running  ETO  within  the  QUICKBASIC  4.0  Environment 

The  program  was  developed  and  tested  using  this  compiler.  It  places  the  user  in 
direct  contact  with  the  source  code,  allowing  modifications  to  be  made  even  within  a 
run.  Also,  a  "scratch  space"  can  be  utilized  within  the  source  code  to  temporarily 
overwrite/replace  variables  that  were  read  in  from  the  input  data  file.  For 
example,  between  runs,  fuel  weight,  launch  weight,  or  other  variables  night  be 
adjusted  if  program  run  results  indicate  the  need.  Values  written  into  the  "scratch 
space"  later  can  be  moved  to  the  input  data  file.  The  Quick  BASIC  environment 
allows  other  files,  such  as  an  input  data  file  to  be  loaded,  modified,  saved  and 
remain  resident,  while  the  main  program  is  still  in  memory. 

To  run  the  program  in  the  Quick  BASIC  editor/compiler  environment,  you  need  to 
follow  the  compiler  instructions  for  loading  Quick  BASIC  and  ET013.BAS.  After 
selecting  the  RUN  option  and  pressing  ENTER,  the  prompts  and  required  responses  are 
exactly  like  those  described  previously  in  this  section.  Screen  output  and  printer 
output  are  identical  also.  Be  sure  to  use  a  copy  of  the  distribution  disk,  not  the 
original.  Also,  the  MS-DOS  COMMAND.COM  for  your  computer  should  be  on  the  disk  or 
default  drive.  We  have  found  that  100  percent  IBM  compatibil ity  is  needed  to  use 
the  Quick  BASIC  editor/compiler  environment. 

3.4  Input  Data  Files  for  a  Single-  or  First-Stage  Vehicle 

All  of  the  ^formation  about  the  vehicle  to  be  simulated  is  contained  in  an  input 
data  file  (two  files  for  a  tv/o-stage  vehicle)  which  is/are  opened  by  the  program  and 
the  information  read  and  stored  in  memory.  The  input  file(s)  consists  of  nine  data 
tables  plus  additional  vehicle  information.  An  input  file  is  normally  less  than 
?,000  bytes,  (about  three  quarters  of  a  page). 
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Input  File  Information: 

The  first,  set  of  information  is  two  lines  of  descriptive  text  about  the  data  file, 
such  as  name,  type  of  vehicle,  etc.  The  program  reads  these  two  lines  as  string 
variables,  so  any  ASCII  characters  may  he  used.  Each  line  should  be  limited  to  80 
characters.  Next  follows  nine  data  tables  (Tables  1  through  9).  See  Figure  3-1* 
and  the  input  files  at  the  end  of  Section  3.6. 

Table  1.  Airbreather  Specific  Impulse 

This  table  is  the  airbreather  specific  impulse,  ISPA.  This  table  must  be 
constructed  as  follows: 

ISPA,  2 

n  (number  of  Mach  args),  m  (number  of  phi  args) 

nls  ri2,  ri3,  . —  fin 

$1,  02,  03 ,  - 0m 

Ispj,  Ispp,  Isp3,  .  Ispn  at  fij 

Ispi,  Ispo,  Isp3,  - . Ispn  at  02 

I sp1 ,  Isp?,  Isp3, .  I spn  at  0m 

The  above  example  is  for  a  2-D  table.  For  a  1-0  table,  use  the  following  format: 
ISPA,  1 

n,  number  of  Mach  args. 


Ispj,  Ispo,  I sp3 ,  .  Ispn 

The  program  uses  the  INPUT#  statement  so  it  exnects  all  values  to  be  separated  by 
commas.  If  more  than  one  line  is  required  to  list  a  group  of  variables,  then  the 
data  can  be  continued  on  the  next  line.  0o  not  use  a  comma  as  the  last  character  in 
a  line  to  be  continued. 


*See  page  35. 
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Table  ?..  Airbreathor  Maximum  Fuel  Equivalence  Ratio 


The  next  data  field  is  the  airbreathor  maximum  allowable  fuel  equivalence  ratio  as  a 
function  of  Mach  number: 

This  is  a  1-D  table: 

PHIMAX,  1 

n,  number  of  Mach  args. 

^1»  f12’  [13’  f1n 

ClmaXp  0max2,  0max3, - 0nax 

Following  the  sane  format  as  the  above, the  remainder  of  the  tabular  data  follows: 
Table  3.  Zero  Lift  Drag  Coefficient 
CD0 ,  1 

n  (number  of  Mach  args) 

Mp  Mp  M3,  . Mn 

Table  4.  Friction  Drag  Increment  as  a  Function  of  Altitude 
DELCD ,  1 

n  (number  of  altitude  args) 
ho, ,  ho„,  ho-,,  . ho„ 

i  t.  o  n 

ACDP  ACq2,  ACd3,  .  . . AC0n 

Table  5.  Lift  Curve  Slope  as  a  Function  of  Mach  Number 
CLALPHA,  1 

n  (number  of  Mach  args) 

Mp  11,,  M3,  . Mn 

CLal  ’  CLa2’  CLa3’  . CLan 
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Table  6.  Value  of  Drag  Due  to  Lift  Factor  K,  as  a  Function  of  Mach  Number 


K,  1 

n  (number  of  Mach  arcs) 


ft,,  ft.,  ft 


'1 


?.’  "3’ 

Ko,  K,, 


fin 

Kn 


Table  7,  Airbreather  Inlet  Capture  Area  Ratio  vs  Mach  Nunber  and  Angle  of 
Attack 

Th'S  is  most  often  a  3-D  table.  The  example  data  file  TEST1.DAT,  however,  shows  how 
it  is  made  effectively  a  1-D  table  by  making  the  Ao/Ac  values  invariant  with  alpha. 

ApAC,  2 

n  (number  of  Mach  args),  n  (number  of  o  args) 

ft.,  f12,  M3,  — . fin 

°r  a2’  a3 . am 

AOACj,  A0AC2,  A0AC3,  - A0AC  at  a1 

A0AC j ,  A0AC2,  A0AC3,  . A0AC  at  a2 


A0AC^ ,  A0AC2,  A0AC3,  - A0AC  at  am 

Table  A.  Rocket  Specific  Impulse  vs  Altitude 
ISPR,  1 

n  (nunber  of  altitudes) 


hop  ho2,  ho3,  - . hon 

IspP.l’  IspR2’  IspR3’ . .  IspRn 


Th's  table  allows  the  rocket  performance  to  increase  with  altitude  in  accordance 
wi*h  the  increasing  nozzle  pressure  ratio. 

Table  9.  Maximum  Rocket.  Propellant  Flow  Rate  vs  Mach  Humber 
’/nOTPMAX ,  1 

n  (number  of  Mach  args) 
rij,  n2,  n3,  . fin 

WDOTPMAX. ,  WDOTPMAX,,  WDOTPMAX,,  .  WDOTPMAX 

1  c.  3  n 
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Note  that  the  data  tables  do  not  require  uniform  increments  for  the  arguments.  Note 
also  that  a  step  change  in  *he  value  of  a  variable  cm  be  used,  e.g.,  a  rocket  can 
he  turned  on  at  Mach  ft  with  a  propellant,  flow  of  600  lbs/sec  by  having  WDOTPMAX  go 
from  0  to  600  between  Mach  8  and  ft .01.  Linear  interpolation  is  used  between  data 
points . 

Following  the  tabular  input  values  the  following  data  items  must  be  entered: 

AC,  AREF 

AC  ts  the  airbreather  inlet  capture  area  in  sc  fee* 

AREF  is  the  vehicle  reference  area  in  sq  feet. 

WLAUNCH ,  V.'FU EL ,  WFIIIAL,  VFII1AL ,  STAGE 
liLADNCH  is  launch  weight,  lbs 
WF1IEL  is  fuel  plus  rocket  propellant,  lbs 
WFIf.’AL  is  burnout  weight  of  stage,  lbs 

VFItlAL  is  the  desired  fina1  velocity  for  this  stage.  For  a  single  stage  to  orbit 
vehicle  use  a  value  *«'"■  d  the  polar  orhit  velocity,  say  27,000  fps  unless  a  smaller 
velocity  is  desired  for  a  lower  inclination  crb’>. 

STAGE  is  1  fcr  a  single-stage  vehicV  or  the  first-stage  of  a  two-stage  vehicle. 

110,  H0,  GAMMA,  ALPHA,  FT,  DELPRINT,  TIMEX 

M0  is  the  starting  Mach  number,  a  value  of  0.001  is  suggested  instead  of  zero  for  a 
stationary  launch.  For  an  air  launch,  use  the  launch  Mach  number.  For  a  low  angle 
ground  launch,  use  the  Mach  number  at  the  end  of  the  rail  (assuming  some  cold  launch 
system  is  used  for  the  initial  acceleration).  H0  is  the  launch  altitude  in  feet. 

Use  0.001  instead  of  a  zero  value. 

GAMMA  is  the  initial  flight  path  angle,  degrees.  A  value  of  0.0  may  be  used  for  a 
horizontal  launch. 

ALPHA  is  the  initial  angle  of  attack,  degrees.  A  value  of  0.0  may  be  used  for  a 
horizontal  launch. 
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DT  is  the  integration  interval,  secs.  The  prngran  is  designed  to  use  large  values, 
i.e.,  2  to  5  secs.  If  the  flight  seems  to  exhibit  large,  erratic  angle  of  attack 
excursions  in  Phase  III,  then  t.rv  smaller  values  ot  DT. 

DELPRINT  is  the  frequency  at  which  integration  steps  are  displayed  to  the  screen.  A 
value  of  20  to  40  insures  that  program  execution  is  not  being  controlled  by  the  rate 
that  the  screen  can  be  updated. 

TIMEX  is  the  starting  value  for  the  program  tine  in  seconds.  8.0  may  be  used  to 
start  the  problem. 

QCOHD ,  Q0HAX,  P0MIH,  ALPHAMAX,  ALPHAHIN 

QCOMD  is  the  commanded  dynamic  pressure  during  Phase  III  flight,  psf. 

O0MAX  anG  Q0MIM  are  values  of  q  greater  and  lesser,  respectively,  than  QCOMD  that 
the  vehicle  is  limited  to. 

ALPHAMAX  and  ALPHAHIN  are  maximum  and  minimum  values  that  are  permitted  during 
Phases  III  and  IV. 

VTAKEOFF,  ALPHA2HAX 


VTAKEOFF  is  the  takeoff  rotation  speed  in  foet/'sec. 

ALPHA2HAX  is  the  maximum  angle  of  attack  permitted  during  Phase  II,  the  interval 
from  rotation  to  Phase  III  flight.  Normally  this  angle  has  its  maximum  just  as 
Phase  II  begins.  It  becomes  the  vehicle  rotation  angle. 

LOADFAC,  ACCCOMD,  GAMMAMAX 


LOADFAC  is  the  maximum  nornal  load  factor  permitted  during  the  Phase  II  climb,  in 
g's.  The  angle  of  attack  is  permitted  to  become  as  large  as  ALPHA2MAX  in  order  to 
generate  the  maximum  LOADFAC. 
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ACCCOMD  is  the  commanded  maximum  axial  acceleration,  in  g's,  during  Phases  I,  II,  and 
III. 

GAM'IAIIAX  is  maximum  flight  path  angle,  in  degrees,  permitted  during  any  flight  phase 
except  Phases  IV  and  V. 

FASTOIC 

FASTOIC  is  the  stoichiometric  fuel-air  ratio  for  the  airbreath’ng  propulsion. 

V0TEHP,  Q0FINAL 

V0TE11P  is  the  flight  velocity,  in  ft/sec  during  Phase  III  at  which  the  vehicle 
begins  to  command  to  a  final  dynamic  pressure  O0F IMAL .  This  variation  is  made 
linearly  between  V0TEMP  and  VFIMAL.  A  combination  of  Y0TEMP  and  Q0FIMAL  is  usually 
selected  so  the  vehicle  will  ascend  along  an  approximation  of  a  constant  vehicle 
skin  temperature  paxh.  Enter  a  value  of  0.0  for  V0TEMP  and  the  same  value  as  QCOMD 
for  Q0FIUAL  if  constant  g  climb  to  VFIHAL  is  desired. 

V0CRUISE,  SWITCH4,  GAMitMMAX 

V0CRUISE  is  the  velocity,  in  ft./soc,  that  nay  be  selected  for  a  constant  velocity 
cruise  to  fuel  exhaustion.  This  is  Phase  IV.  Enter  a  zero  if  there  is  no  cruise 
phase. 

SWITCH^  has  a  value  of  either  0  or  1.  Zero  selects  cruise  at.  constant  altitude. 

The  vehicle  maintains  the  altitude  it  attained  during  the  Phase  III  climb.  A  value 
of  1  selects  cruise  at  maximum  lift-to-drag  ratio.  Enter  a  zero  if  there  is  no 
cruise. 

GAf  1MA4MAX  is  the  maximum  flight  path  angle  (in  decrees)  permitted  during  a  maximum 
l/D  cruise  (SWITCH4=1).  A  small  value,  irt  the  range  of  0.?0  to  0.50  degrees  is 
suggested.  The  vaiue  should  be  chospn  so  that  the  vehicle  performs  a  fairly  smooth 
climb/cruise  only  occasionally  encountering  the  limiting  value  after  experiencing 
some  flight  path  oscillations  early  in  the  cruise.  This  variable  has  no 


signif icance  for  a  constant  altitude  cruise  (SWITCH4=0)  and  a  value  of  zero  should 
be  entered. 


FLAGS,  Vg5,  ALPHA5MAX,  GAMMAS ,  ACCC0MD5 

These  five  variables  are  used  to  prescribe  a  pull  up  and  ballistic  ascent  after 
completing  part  of  the  ascent  after  another  phase.  This  is  Phase  V. 

FLAG5  can  have  a  value  of  0  or  1.  Zero  m<=ans  that  a  ballistic  ascent  is  not  go ’no 
to  be  used.  In  this  case,  enter  zeros  for  all  the  remaining  variables  on  this  line. 
If  FLAGS  is  set  to  1,  then  enter  values  for  the  following: 

V0S  is  the  velocity,  ft/sec,  at  which  pull  up  is  to  begin.  ALPHACMAX  is  the 
commanded  angle  of  attack,  in  degrees,  to  be  employed  during  pull  up. 

GAMMAS  is  the  flight  path  angle,  in  degrees,  to  be  achieved  during  the  pull  up.  As 
seen  as  this  flight  path  angle  is  achieved,  annle  of  attack  is  commanded  to  zero  to 
permit  the  flight  to  continue  ball istically. 

ACCC0MD5  is  the  axial  acceleration,  g's,  commanded  during  Phase  V.  Thrust  is  reduced 
if  required  to  avoid  exceeding  this  commanded  acceleration. 

Al,  A2,  Bl,  Cl,  CP,  Dl,  DELH,  V0P,  H0P 


The  first  seven  variables  are  ad;:uctment  coefficients  used  in  the  Phase  III 
navigation  algorithm,  GPID13  (see  Figure  ?-4  for  flow  chart  of  control  logic). 

These  coefficients  are  used  to  set  the  sensitivity  of  angle  of  attack  (a)  while  the 
vehicle  is  trying  to  fo1 1  ow  a  constant  d'/nanic  pressure  (q)  pa+h.  Suggested  values 
are,  in  order: 

0.1,  0.1,  0.009,  0.1,  2.0,  0.1,  50.0 
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As  a  further  suggestion,  set  DELQ  at  about  5  percent  of  the  commanded  Phase  III 
dynamic  pressure,  QCOMD.  Some  work  has  been  done  to  implement  a  simpler  Phase  III 
algorithm,  GUID14,  (included  on  distribution  disk  as  GUID14.BAS),  which  has 
constants  "hardwired"  into  the  code.  Whether  you  are  using  GUID13  or  GUID14,  this 
is  an  area  where  some  experimentation  may  be  needed  to  obtain  smooth  constant  q 
trajectory  paths. 

If  the  vehicle  deviates  from  a  smooth  constant  q  trajectory,  the  following 
information  will  assist  you  in  making  the  necessary  changes.  There  are  six  possible 
deviations  that  may  occur  as  shown  in  the  following  sketch: 

_ _  QCOMO  -  oao 

ALTITUDE 


VELOCITY 

Note  that  the  value  DELQ  is  used  to  determine  if  the  vehicle  reaches  deviations  2 
and  5.  In  order  to  correct  each  of  these  deviations  (1  to  6),  a  must  be  changed 
such  that  the  vehicle  returns  to  the  desired  q  path.  For  each  different  deviation, 
there  is  a  different  a  sensitivity  that  is  required  to  return  to  the  desired  dynamic 
pressure  path.  If  the  vehicle  is  having  difficulty  recovering  from  any  of  the 
deviations,  identify  the  deviation  condition  that  exists  and  slightly  vary  the 
corresponding  adjustment  coefficient. 

The  adjustment  coefficients  affect  the  angle  of  attack  in  the  following  manner. 
DESCRIPTION  "q  is"  REACTION  CALCULATION 

(1)  Higher  than  QCOMD  i  Increasing  Increase  a 

(2)  Higher  than  QCOMD  &  Outside  DELQ  Increase  a 


|  a  =  a 
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DESCRIPTION  "q  is" 


(3)  Higher  than  QCOMD  &  Decreasing 

(4)  Lower  than  QCOMD  &  Decreasing 

(5)  Lower  than  QCOMD  &  Outside  DELQ 

(6)  Lower  than  QCOMD  &  Increasing 

Note  that  in  conditions  3  and  6  the  control  logic  is  attempting  to  prevent 
overshooting  QCOMD. 

V0 2  and  H02  are  a  velocity  (ft/sec)  and  altitude  (feet)  pair  which  defines  the  end 
of  the  Phase  II  constant  load  factor  climb.  The  vehicle  is  not  permitted  to  exceed 
V02  until  the  altitude,  H02  has  been  exceeded.  For  an  aircraft,  V02  values  of  700 
to  900  ft/sec  are  usually  chosen.  H02  can  range  from  10,000  feet  to  30,000  feet.  A 
combination  of  H02,  V02  values  is  usually  chosen  that  would  be  consistent  with  an 
optimum  subsonic-transonic  climb,  such  as  would  be  prescribed  by  a  Rutowski  path 
optimization  method. 

3.5  Input  Data  Files  for  a  Second-Stage  Vehicle 

The  nine  data  tables  described  in  the  previous  section  are  required  and  are  used  in 
the  same  manner  for  a  second-stage  vehicle. 

The  remaining  values  required  are 
AC,  AREF 

WEIGHT,  WFUEL,  WFINAL ,  VFINAL ,  STAGE 
DT,  DELPRINT 

Q0COMD ,  Q0MAX,  Q0MIN,  ALPHAMAX,  ALPHAMIN 
LOADFAC ,  ACCCOMD,  GAMMAMAX 
FAST IOC 

V0TEMP,  Q0FINAL 

V0CRUISE,  SWITCH4,  GAMMA4MAX 

FLAG5,  V05,  ALPHA5MAX ,  GAMMA5 ,  ACCC0MD5 

Al,  A2,  Bl,  Cl,  C2,  Dl,  DELO,  V02,  H02 


REACTION 


CALCULATION 


Decrease  a  a  =  a 
Decrease  a) 
Decrease  aj 
Increase  a 


Q-QCOMD 

QCOMD 

C2  '  3t 
Q-QCOMD 
QCOMD 
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Note  that  the  variables  VTAKEOFF  and  ALPHA2MAX  are  not  used,  and  since  the  flight  is 
continuing  from  the  previous  stage,  starting  values  of  M0,  HJD,  GAMMA,  ALPHA,  and 
TIMEX  are  not  input.  When  the  second-stage  data  file  is  loaded,  its  values 
overwrite  the  first-stage  data  and  current  values  of  appropriate  program  variables. 
As  an  example,  the  current  value  of  WEIGHT  is  replaced  by  the  value  of  WEIGHT  read 
from  the  second-stage  file. 

3.6  Sample  Input  Data  Files  and  Program  Output 

This  section  describes  the  example  input  data  files  contained  on  the  distribution 
disk.  Each  example  illustrates  features  typically  used  in  combination.  The  example 
input  data  files  are  completely  arbitrary  and  are  not  intended  to  represent  any 
known  or  projected  vehicle  or  level  of  technology.  No  inference  should  be  made  that 
the  example  vehicles  typify  current  or  projected  capabilities. 

The  file  TEST1.DAT  is  the  simplest  example.  It  is  a  highly  speculative 
single-stage-to-orbit  horizontal  takeoff  vehicle  having  all -airbreathing  propulsion 
to  orbital  velocity.  The  vehicle  is  commanded  to  fly  to  the  orbital  condition  at  a 

O 

constant  1,000  lb/ft1'  dynamic  pressure.  The  vehicle  has  a  gross  takeoff  weight  of 
100,000  lbs  and  has  60,000  lbs  of  fuel. 

The  file  STG1.DAT  and  STG2.DAT  define  a  two-stage  vehicle  intended  to  fly  to  orbit 

from  a  horizontal  takeoff.  The  first  stage  uses  airbreathing  propulsion  up  to  the 

Mach  8  staging  point.  The  second  stage  operates  with  rocket  propulsion  from  staging 

to  orbit.  The  gross  takeoff  weight  of  the  vehicle  is  400,000  lbs,  of  which  150,000 

lbs  is  first-stage  fuel.  The  second  stage  weighs  150,000  lbs  at  staging,  of  which 

120,000  lbs  is  propellant.  The  empty  weight  of  the  first  stage  is  100,000  lbs.  The 

2 

first  stage  accelerates  along  a  constant  1,000  Ib/ft  dynamic  pressure  trajectory  to 
the  staging  condition. 

The  file  STG2BAL.DAT  may  be  used  as  an  alternative  to  STG2.DAT,  above,  as  a  second 
stage.  This  vehicle  is  identical  to  STG2.DAT  but  employs  a  pull  up  after  staging  and 
ballistic  ascent  to  orbit. 

The  file  CRUZ2.DAT  is  a  Mach  4  cruise  aircraft.  It  has  a  gross  takeoff  weight  of 
400,000  lbs,  of  which  170,000  lbs  is  fuel.  It  employs  airbreathing  propulsion  in 
all  flight  phases.  It  cruises  at  maximum  L/D. 
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A  maximum  L/D  cruise  may  experience  oscillations  in  altitude  and  flight  path  angle 
if  the  vehicle  does  not  start  the  cruise  phase  at  lift  nearly  equal  to  weight  and  at 
a  very  low  flight  path  angle.  The  code  includes  an  algorithm  that  smooths  the 
climb/cruise  flight  path  by  reducing  angle  of  attack  whenever  GAM11A4MAX  is  exceeded 
and  whenever  lift  becomes  less  than  weight.  Also,  if  the  flight  path  angle  becomes 
negative,  angle  of  attack  is  increased.  This  algorithm  seems  to  produce  fairly 
smooth  climb/cruise  paths  once  it  has  forced  the  vehicle  to  a  fairly  steady  lift 
equal  weight  condition. 

The  program  tests  if  the  vehicle  input  data  files  have  been  read  correctly.  In 
starting  the  program,  the  following  message  will  appear  on  the  screen  just  above  the 
prompt  for  creation  of  a  plotter/printer  file: 

1ST/SINGLE  (OR  2ND)  STAGE  INPUT  FILE  HAS  BEEN  READ.  LAST  ITEM  IN  INPUT  FILE  IS  H02, 
VALUE  READ  FROM  FILE  IS  18000.  IF  THIS  VALUE  IS  INCORRECT,  INPUT  FILE  IS  IN  ERROR. 

In  all  of  the  example  vehicle  files,  the  value  of  H02  is  18,000.  The  screen  message 
above  indicates  that  the  program  has  read  all  of  the  data  correctly.  If  you  see  an 
unexpected  value,  this  means  that  the  number  of  data  items  in  the  file  exceeds  the 
number  that  the  program  was  directed  to  read.  Check  to  see  that  the  number  of 
entries  in  the  nine  data  tables  correspond  correctly  with  the  declared  number  and 
that  the  data  entered  after  the  tabular  input  is  correct.  The  problem  will  be  found 
to  be  too  many  entries.  A  duplicate  line  of  data  or  extra  commas  is  the  usual 
cause.  An  error  of  too  many  data  items  will  not  prevent  the  program  from  running, 
but  erroneous  results  will  immediately  be  apparent.  The  other  type  of  error,  not 
enough  data  items,  will  not  permit  the  program  to  get  as  far  as  the  prompt  above. 

An  error  message  "Input  past  end"  and  program  termination  will  occur  if  this  error 
is  present.  In  this  case  the  program  is  looking  for  data  which  is  not  present.  A 
line  of  data  items  omitted,  or  a  lacking  comma  are  the  usual  causes.  The  program  is 
quite  exacting  in  its  requirement  that  the  input  file  correspond  precisely  to  the 
declared  number  of  items  to  be  read.  Values,  separated  by  commas,  must  appear 
exactly  in  the  prescribed  format.  If  an  option  is  not  to  be  employed,  or  a  zero 
value  is  planned  for  a  variable,  a  numerical  entry  must  be  present  (even  a  zero). 
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PROGRAM  "ETO 


Sample  Input  Data  Files 

INPUT  DATA  FILE  "TEST1.DAT"  1  SEPT88 
A/B  TO  ORBIT  AT  CONSTANT  Q=1  000 
ISPA.2 
12  3 

0.0  ,  .5  ,.9  ,1.5  ,2.0  ,2.5  ,3.0  ,4.0  ,8.0  ,12.0  ,20.0  ,28.0 

0.1  ,1.0  ,10.0 

4000,  4000,3500,3300,3650,3750,3700,3500,3000,1800,800,600 

3000,  3000,3000,3300,3650,3750,3700,3500,3000,1800,800,600 

500,  500,  500,  500,  500,  500,  500,  500,  600,500,400,400 

PHIMAX, 1 
10 

0.0  ,0.8  ,1.0  ,1.4  ,2. 0,8.0  ,16.0  ,18.0  ,25.0  ,28.0 
1.0  ,1.0  ,1.0  ,1.0  ,1.0, 1.0  ,  3.0  ,  3.0  ,  5.0  ,  5.0 
CD0, 1 
15 

0.0,  .4  ,  .8  ,  .9  ,1.0  ,1.1  ,1.2  ,1.5  ,2.0  ,3.0  ,5.0  ,10.0  ,15.0  ,20.0  ,25.0 
0.0,. 010,. 010,. 010,. 023,. 022,. 018,. 016,. 013,. 008,. 004,. 002,. 0019,. 0017,. 0015 
DELCD, 1 
5 

0.0,50000. ,100000. ,150000. ,300000. 

0.0, .0006  ,.002  ,.001  ,.0030 

CLALPHA , 1 
10 


0.0, 

1.5 

,2.0  , 

3.0  ,  4 

•0  , 

5.0 

,6.0  ,8.0  ,15.0  ,25.0 

0.040, 

.040, 

.039  , 

.024,  . 

020, 

.018, 

.017, .016,  .010  ,.010 

K,1 

0.0, 

1.5 

,2.0 

,  3.0 

,4.0 

,5.0 

,6.0  ,8.0  ,15.0  ,25.0 

0.400 

,  0.400 

,.440 

,.730  , 

.870, 

,.970, 

1.025,1.090  ,1.745  ,1.745 

A0AC.2 

10,3 

0.0  , 

0.01,  0. 

1  , 

0.5  , 

9.8  , 

1.0  ,  10.0  ,  15,  15.01 

00 

CM 

.0 

-10.0 

,0.0  ,10.0 

80.0, 

2.00,  0 

.50, 

0.15, 

.15, 

.15  ,  .50,  1.0 

,1.0,1 

.0 

80.0, 

2.00,  0 

.50, 

0.15, 

.15, 

.15  ,  .50,  1.0 

,1.0,1 

.0 

80.0, 

2.00,  0 

•  50, 

0.15, 

.15, 

.15  ,  .50,  1.0 

,1.0,1 

.0 

ISPR, 1 
4 

0.0  ,  30000.0  ,  80000.0  ,350000.0 

390  ,440  ,445  ,460 

WDOTPMAX , 1 
4 

0.0  ,  8.0  ,8.01  ,28.0 

0.0  ,  0.0  ,0.0  ,0.0 

100.0  ,  1875.0 

100000., 60000., 40000.  ,27000  ,1 
0.001  ,0.001  ,0.0  ,0.0  ,  5.0  ,  10.00  ,  0.0 

1000.0,1100.0,900.0,  5.0  ,  -5.0 
375.0  ,5.0 
2.00,  0.50  ,  8.0 

0.0292 
0  ,1000 
0  ,0  ,0 
0  0  0  0  0 

.1  ,’.r,  ! 009 ,  .1,  2.0  ,  .1  ,50  ,900  ,18000 
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INPUT  DATA  FILE  ”  STG1 .DAT"/26JUNE88/1530/CONST  Q  TO  ORBIT  / PROGRAM" ET02" 
ALL  AIRBREATHING  PROPULSION  TO  MACH  8  /FIRST  STAGE 
ISPA,  2 
13,3 

0.0  ,  0.5  ,1.0,  1.5  ,2.0  ,2.5  ,3.0  ,4.0,  5.0  ,  5.1  ,  8.0  ,  15.0  ,  26.0 

0.1  ,1.0  ,10.0 

4500,  4400,4200,4000,3700,3500,3200,2800,2400,  2400,  2300,  1000,  600 

4500,  4400,4200,4000,3700,3500,3200,2800,2400,  2400,  2300,  1000,  500 

1000,  1000,1000,1000,1000,1000,1000,  800,  600,  600,  600,  450,  450 

PHIMAX , 1 
12 

0.0  ,0.8  ,1.0  ,1.4  ,2.0,  5.0,  5.1,  8.0  ,16.0  ,18.0  ,25.0  ,28.0 
1.0  ,1.0  ,1.0  ,1.0  ,1.0,  1.0,  1.0  ,  1.0,  3.0  ,  5.0  ,  5.0  ,  5.0 
CD0, 1 
15 

0.0, .4  ,.8  ,.9  ,1.0  ,1.1  ,1.2  ,1.5  ,2.0  ,3.0  ,5.0  ,10.0  ,15.0  ,20.0  ,25.0 
0.0,. 009,. 009,. 009,. 023,. 022,. 018,. 016,. 013,. 008,. 004,. 002,. 0019,. 0017,. 0015 
DELCD, 1 
5 

0.0, 50000 . , 100000 . , 1 50000 . , 200000 . 

0.0, .0006  ,.002  ,.001  ,.0020 
CLALPHA, 1 
10 


0.0, 

1.5  ,2.0 

,3.0  ,  4.0  ,  5.0  ,6.0  ,8.0 

,15.0  ,25.0 

0.043, 

.043, .0391, .024,  .020,  .018 , .017 , . 016 , 

.010  ,.010 

K,1 

10 

0.0, 

1.5  ,2.0 

,3.0  ,  4.0  ,  5.0  ,6.0  ,8.0 

,15.0  ,25.0 

0.407 

,  0.407, .447 

,.728, .873  ,.970  ,1.026,1.091 

,1.745,1.745 

A0AC , 2 
9,3 

0.0  , 

0.01,  0.1  , 

0.5  ,  0.8  ,  1.0  , 

10.0  ,  25.0 

,28.0 

-10.0 

80.0, 

,0.0  ,10.0 

2.00,  0.50, 

0.15,  .15,  .15 

,  .51,  1.0 

,1.0 

80.0, 

2.00,  0.50, 

2.00,  0.50, 

0.15,  .15,  .15 

,  .51,  1.0 

,1.0 

80.0, 

0.15,  .15,  .15 

,  .51,  1.0 

,1.0 

ISPR ,  1 


4 


*  , 

30000, 

80000, 

350000 

390  , 

440, 

445, 

460 

WDOTPMAX , 1 

/. 

0.0  , 

8.0  , 

8.01  , 

28 

0.0  , 

0.0  , 

0.0  , 

0.0 

400.0 

,  5500 

.0 

400000 

,150000,250000 

,8000  ,  1 

0.001 

,0.001 

,0.0  ,0.0 

,  2.0, 

10.0  ,  0.0 

1500.0 

,1600.0 

,1400.0, 

5.0  ,  -5.0 

390.0 

,15.0 

2.00, 

•  50  , 

8.0 

0.0292 

0  ,  1000 
0  ,0  ,0 

0  ,0  ,0  ,0  ,0 

•1  ,.l  ,  *009  ,  .1  ,  2.0  ,  .1  ,  50.0,900  ,  18000 


INPUT  DATA  FILE  "STG2 . DAT" /26JUNE88/ 1530/ CONSTANT  Q  TO  ORBIT  / PROGRAM” ET02” 

ALL  ROCKET  SECOND  STAGE 

ISPA.2 

13,3 

0.0  ,  0.5  ,1.0,  1.5  ,2.0  ,2.5  ,3.0  ,4.0,  5.0  ,  5.1  ,  8.0  ,  15.0  ,  26.0 

0.1  ,1.0  ,10.0 

4500,  4400,4200,4000,3700,3500,3200,2800,2400,  2400,  2300,  1000,  600 

4500,  4400,4200,4000,3700,3500,3200,2800,2400,  2400,  2300,  1000,  500 

1000,  1000,1000,1000,1000,1000,1000,  800,  600,  600,  600,  450,  450 

PHIMAX, 1 
12 

0.0  ,0.8  ,1.0  ,1.4  ,2.0,  5.0,  5.1,  8.0  ,16.0  ,18.0  ,25.0  ,28.0 
0.0  ,0.0  ,0.0  ,0.0  ,0.0,  0.0,  0.0  ,  0.0,  0.0  ,  0.0  ,  0.0  ,  0.0 
CD0, 1 
15 

0.0, .4  ,.8  ,.9  ,1.0  ,1.1  ,1.2  ,1.5  ,2.0  ,3.0  ,5.0  ,10.0  ,15.0  ,20.0  ,25.0 
0.0,. 009,. 009,. 009,. 023,. 022,. 018,. 016,. 013,. 008,. 004,. 002,. 0019,. 0017,. 0015 

DELCD, 1 
5 

0.0, 50000 . , 100000 . , 1 50000 . , 200000 . 

0.0, .0006  ,.002  ,.001  ,.0020 
CLALPHA , 1 
10 

0.0,  1-5  ,2.0  ,3.0  ,  4.0  ,  5.0  ,6.0  ,8.0  ,15.0  ,25.0 

0.043,  .043, .0391, .024,  .020,  . 018 , .017 , . 016 ,  .010  ,.010 

K,1 
10 


0.0,  1.5  ,2.0 

,3.0  ,  4.0  ,  5.0 

,6.0  ,8.0 

,15.0 

,25.0 

0.407  ,  0.407, .447 

A0AC.2 

9,3 

,.728, .873  ,.970 

,1.026,1.091 

,1.745,1.745 

0.0  ,  0.01,  0.1  , 
-10.0  ,0.0  ,10.0 

oo 

1.0  , 

10.0  ,  25.0 

,28.0 

80.0,  2.00,  0.50, 

0.15,  .15, 

.15 

,  .51, 

1.0 

,1.0 

80.0,  2.00,  0.50, 

0.15,  .15, 

.15 

,  .51, 

1.0 

,1.0 

80.0,  2.00,  0.50, 

0.15,  .15, 

.15 

,  .51, 

1.0 

,1.0 

ISPR.l 

4 

0  , 

30000, 

80000, 

350000 

390  , 

440, 

445, 

460 

WDOTPMAX 

4 

0.0  , 

,1 

8.0  , 

8.01  , 

28 

600  , 
100.0  , 

600  , 

600  , 

600 

1500.0 

150000  ,120000,  30000.0  ,27000  ,  2 

2.0  10.0 

1000.'0,1100. 0,900.0,  5.0  ,  -5.0 
2.00,  1.00  ,  8.0 
0.0292 

0  ,  1000 

0  ,0  ,0 

0  ,0  ,0  ,0  ,0 

.1  ,.l  ,  .009  ,  .1  ,  2.0  ,  .1  ,  50.0,900  ,  18000 
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INPUT  DATA  PILE  "STG2BAL.DAT"/29DEC88/BALLISTIC  ASCENT  TO  ORBIT/"ETO“ 

ALL  ROCKET  SECOND  STAGE 

ISPA.2 

13,3 

0.0  ,  0.5  ,1.0,  1.5  ,2.0  ,2.5  ,3.0  ,4.0,  5.0  ,  5.1  ,  8.0  ,  15.0  ,  26.0 

01  10  10  0 

4500’  4400,4200,4000,3700,3500,3200,2800,2400,  2400,  2300,  1000,  600 

4500,  4400,4200,4000,3700,3500,3200,2800,2400,  2400,  2300,  1000,  500 

1000,  1000,1000, 1000, 1000, 1000,1000,  800,  600,  600,  600,  450,  450 

PHIMAX, 1 
12 

0.0  ,0.8  ,1.0  ,1.4  ,2.0,  5.0,  5.1,  8.0  ,16.0  ,18.0  ,25.0  ,28.0 
0.0  ,0.0  ,0.0  ,0.0  ,0.0,  0.0,  0.0  ,  0-0,  0-0  ,  0.0  ,  0.0  ,  0.0 
CD0, 1 
15 

0.0, .4  ,.8  ,.9  ,1.0  ,1.1  ,1.2  ,1.5  ,2.0  ,3.0  ,5.0  ,10.0  ,15.0  ,20.0  ,25.0 
0.0,. 009,. 009,. 009,. 023,. 022,. 018,. 016,. 013,. 008,. 004,. 002,. 0019,. 0017,. 0015 
DELCD, 1 
5 

0.0, 50000., 100000., 150000., 200000. 

0.0,. 0006  ,.002  ,.001  ,.0020 
CLALPHA , 1 
10 


0.0, 

1.5  ,2.0 

,3.0  ,  4.0 

,  5.0 

,6.0  ,8.0 

,15.0  ,25.0 

0.043, 

.043, .0391, .024,  .020 

,  .018 

,.017, .016, 

.010  ,.010 

K,1 

10 

0.0, 

1.5  ,2.0 

,3.0  ,  4.0 

,  5.0 

,6.0  ,8.0 

,15.0  ,25.0 

0.407 
A0AC , 2 
9,3 

,  0.407, .447 

,.728, .873 

,.970 

,1.026,1.091 

,1.745,1.745 

0.0  , 
-10.0 

0.01,  0.1  , 
,0.0  ,10.0 

0-5  , 

0.8  , 

1.0  , 

10.0  ,  25.0 

80.0, 

2.00,  0.50, 

0.15, 

.15, 

.15 

,  .51,  1.0 

80.0, 

2.00,  0.50, 

0.15, 

•  15, 

.15 

,  .51,  1.0 

80.0, 

2.00,  0.50, 

0.15, 

.15, 

.15 

,  .51,  1.0 

ISPR.l 

4 

0  , 

30000, 

80000 

9 

350000 

390  , 

440, 

445, 

460 

WDOTPMAX, 1 

/. 

0.0  , 

8.0  , 

8.01 

9 

28 

600  , 

600  , 

600 

9 

600 

100.0 

,  1500.0 

150000 

,120000,  30000.0  ,27000 

,  2 

1.0, 

5.0 

1000.0 

,1100.0,900.0, 

5.0  ,  -5.0 

2.00, 

1.00  ,  8.0 

0.0292 

0  ,  1000 
0  ,0  ,0 

1  ,8000  ,10  ,1.5  ,2 

.1  ,.l  ,  .009  ,  .1  ,  2.0  ,  .1  ,  50.0,900  ,  18000 
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INPUT  DATA  FILE  "  CRUZ2 . DAT”/24AUG87/0045/M  4  CRUISE  PROGRAM  "ETO" 

AB  SINGLE  STAGE 

ISPA.2 

14.3 

0.0  ,  0.5  ,1.0,  1.5  ,2.0  ,2.5  ,3.0  ,4.0,  5.0  ,  5.1  ,  8.0  ,  15.0,15.01,28.0 

0.1  ,1.0  ,10.0 

4500,  4400,4200,4000,3700,3500,3200,2800,2400,  2400,  2300,  1000,  0.0  ,0.0 

4500,  4400,4200,4000,3700,3500,3200,2800,2400,  2400,  2300,  1000,  0.0  ,0.0 

1000,  1000,1000,1000,1000,1000,1000,  800,  600,  600,  600,  450  ,  0.0  ,0.0 

PHIMAX, 1 
11 

0.0  ,0.8  ,1.0  ,1.4  ,2.0,  5.0,  5.1,  8.0  ,15.0  ,15.01,28.0 
1.0  ,1.0  ,1.0  ,1.0  ,1.0,  1.0,  1.0  ,  1.0,  3.0  ,  0.0  ,  0.0 
CD0, 1 
15 

0.0, .4  ,.8  ,.9  ,1.0  ,1.1  ,1.2  ,1.5  ,2.0  ,3.0  ,5.0  ,10.0  ,15.0  ,20.0  ,25.0 
0.0,. 009,. 009,. 009,. 023,. 022,. 018,. 016,. 013,. 008,. 004,. 002,. 0019,. 0017,. 0015 
DELCD, 1 
5 

0.0, 50000., 100000., 150000., 200000. 

0.0, .0006  ,.002  ,.001  ,.0020 
CLALPHA, 1 
10 

0.0,  1.5  ,2.0  ,3.0  ,  4.0  ,  5.0  ,6.0  ,8.0  ,15.0  ,25.0 

0.043,  .043, .0391, .024,  .020,  .018 , .017 ,. 016,  .010  ,.010 

K,1 
10 

0.0,  1.5  ,2.0  ,3.0  ,  4.0  ,  5.0  ,6.0  ,8.0  ,15.0  ,25.0 

0.406  ,  0.406, .446, .727, .873,  .97  ,  1.027,1.091,1.745,1.745 

A0AC.2 

10.3 


0.0  , 

0.01, 

0.1  , 

0.5  , 

0.8  , 

1.0  , 

10.0,  15.0 

,15.01 

• 

00 

CN 

w 

-10.0 

,0.0  ,10.0 

80.0, 

2.00, 

0.50, 

0.15, 

.15, 

.15 

,  .51,1.0 

,0.0, 

0.0 

80.0, 

2.00, 

0.50, 

0.15, 

•15, 

.15 

,  .51,1.0 

,0.0, 

0.0 

80.0, 

2.00, 

0.50, 

0.15, 

•15, 

.15 

,  .51,1.0 

,0.0, 

0.0 

ISPR, 1 
4 

0,  30000,  80000,  350000 

390,  440,  445,  460 

WDOTPMAX , 1 
4 

0.0,  8.00,  8.01,  28.0 

0.0,  0.0,  600.0,  600.0 
400.0  ,  5500.0 

400000  ,170000,230000  ,  27000  ,  1 

0.001  ,0.001  ,0.0  ,0.0  ,  2.0,  20.0  ,  0.0 
800.0,  850.0,  750.0,  5.0  ,  -5.0 
390.0  ,15.0 
2.00,  .50  ,  8.0 

0.0292 

0  ,  1000 

4000  ,1  ,0.5 

0  ,  0  ,  0  ,  0  ,  0 
.1  ,.l  ,  .009  ,  .1  ,  2.0  ,  .1  ,  50.0,900  ,  18000 
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3.7  Sample  Program  Output 


The  following  sample  output  is  for  the  two-stage  vehicle  defined  by  the  data  files 
STG1.DAT  and  STG2.DAT.  The  integration  interval  was  2  seconds.  A  FREQ  value  of  50 
was  used;  thus  the  output  is  for  every  100  seconds  of  flight. 


01-29-1989 


1ST /SINGLE  STAGE  INPUT  DATA  FILE:  STG1.DAT 
2ND  STAGE  INPUT  DATA  FILE  :  STG2.DAT 


PLOTTER/PRINTER  FILE  NAME 


:  G : 2STG „ OUT 


- PROGRAM  ET02 

LAUNCH  WEIGHT  --- 
FINAL  WEIGHT  = 
FUEL.  /  PROPELL  ANT  REM  A I N I NG  -• 


RUN  SUMMARY-- . 

400000.00  LBS. 
38045.15  LBS. 
8045.15  LBS. 


VEHICLE  REACHED  ORBITAL  VELOCITY. 


FLIGHT  TIME 

3094 . 7.1 

SECS . 

FINAL.  VELOCITY  == 

25858.04 

F  I /SEC 

FINAL  ALTITUDE  == 

152301.95 

FT 

FINAL  G0 

1049.85 

F'BF 

AVERAGE  I  SR 

.1024. 1.1 

SEC 

MEAN  ISP 

359,.  19 

SEC 

FINAL  CONDITIONS  ARE  SHOWN  BELOW: 


TIME 

V0 

H0 

R 

ISP 

ISPEFF 

ISPMEAN 

WDOTF 

WDOTA 

PHI 

Q8 

ALPHA 

6AMMA 

WEIGHT 

LIFT 

DRAG 

NDOTP 

RESl 

ISPAV8 

3094.71 

25858.04  152301.95 

4447.27 

449.02 

CD 

►O 

~o 

359.19 

0.00 

0.00 

0.00 

3094.71 

1049.85 

-5.00 

0.08 

38045.15 

-78738.77 

10941.28 

66.56 

0.00 

1024.11 

10  s  17:42 


THRUST  THRUSTCOH  THRUSTHAI 
ACCEL-6 'S  ACCCON  T/D 

298B6.49  29886.49  269489.94 
8.50  0.50  2.73 
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INPUT  DATA  FILE  "  STG1 .DAT" /26JUNE80/1 530/CONST  Q  TO  ORBIT  /PROGRAM "ET02" 
ALL.  AIRBREATHING  PROPULSION  TO  MACH  8  /FIRST  STAGE 
I  SPA ,  2 


0 . 0  , 

0  ..  1  , 

0  .  !:•  1 . 0  ,  1.5 

L.0  ,10.0 

2.0  ,2.5  , 

3 ..  0  ,4.0,  5.0  , 

13 .  1 

8.0  , 

i  5.0  , 

*’6  - 

4500, 

4400,4200,4000, 

3/00,3500 , 

3200 , 2800 , 2400 , 

2400, 

2300 , 

1000, 

600 

4500, 

4400,4200,4000, 

3700,3500, 

3200 , 2800 , 2400 , 

2400, 

2300 , 

1000, 

500 

.1000, 

1000,1000,1000, 

1000,1000, 

.1.000,  800,  600, 

600, 

600, 

450, 

450 

PHI MAX , 1 

12 


0.0  ,0.8  ,1.0  ,1.4  ,2.0,  5.0, 
1.0  ,1.0  ,1.0  ,1.0  ,1.0,  1.0, 
CD0 , 1 

15 

0.0, .4  ,.8  ,  .  9  ,1.0  ,1.1  ,1.2 
0.0,. 009 , . 009 , . 009 , . 023 , . 022 , 
DEL.CD ,  1 


5.1,  8.0  ,16.0  ,18.0  ,25.0  ,28.0 
1.0  ,  1.0,  3.0  ,  5.0  ,  5.0  ,  5.0 


,1.5  ,2.0  ,3.0  ,5.0  ,10.0  ,15.0  ,20.0  ,25.0 
018, .016, ,013, .008, .004, .002, .0019, .0017, .001 


0.0, 50000 . , 100000 . , 1 50000 . , 200000 . 
0,0,. 0006  , . 002  , . 001  ,  . 0020 
CL.  ALPHA,  1 
10 


0 . 0  , 

1.5  ,2.0  , 

3.0  ,  4.0  ,  5 . 0 

,6.0  ,8.0 

,15.0  ,25.0 

0 „ 043 , 

L--  1 

.043  ,  ,039.1 , 

.024,  .020,  .,018, 

.017, .016, 

.010  ,.010 

K  ,  ,L 

10 

0  „  0  , 

1.5  ,2.0  , 

3 .0  ,  4.0  ,  5.0 

,6.0  ,8.0 

,15.0  ,25.0 

0 . 407 
A0AC ,  2 

9 , 3 

,  0.407,. 44 7, 

.728,. 873  ,.970  , 

1 ,026, 1.091 

, 1.745,1.745 

0.0  ,  ( 

3,01,  0.1  , 

0.5,  0.8, 

1  . 0  ,, 

10.0  ,  25.0 

,28 

-10.0  ,( 

3.0  ,10.0 

80 . 0  , 

2.00,  0.50, 

0.15,  .15, 

,  15 

,  .51,  1.0 

,  1 . 0 

80.0, 

2.00,  0 . 50 , 

0.15,  .15, 

.  15 

,  .51,  1.0 

,1.0 

80.0, 

I  SF  E ,  1 

4 

2 . 00 ,  0 . 50 , 

0.15,  . 1 5 , 

.  .1  5 

,  . 51 ,  1.0 

,  1.0 

0 

30000 , 

80000 , 

330000 

390  , 

440, 

445, 

4."-.0 

WDOTPMAX , 1 


0 . 0  , 

8.0,  8.01, 

28 

0.0  , 

0.0,  0.0  , 

0.0 

•'1  00 . 0 

,,  5500.0 

400000 

,150000,250000  ,8000  ,  .1 

0.001  ,0.00.1  ,0.0  .,0.0  ,  2.0,  10.0  ,  0.0 

1 500 . 0 , 1 600 .0,1400.0,  5.0  ,  - 5 . 0 
;.90.0  ,15.0 
2.00,  .50  ,  8.0 

0 . 0292 

0  ,  1 000 

0  0  ,0 

0  ,0  ,0  ,0  ,0 

.1  , -1  ,  .009  ,  .1  ,  2.0  ,  .1  ,  50.0,900  ,  18000 
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INPUT  DATA  FILE  "STG2 . DAT" /2&JUNE88/ 1530/CONSTANT  Q  TO  ORBIT  /PROGRAM "E TO 2" 
ALL  ROCKET  SECOND  STAGE 
I SPA, 2 

•it  t; 

.A 

0.0  ,  0.5  ,1.0,  1.5  ,2.0  ,2.5  ,3.0  ,4.0,  5.0  ,  5.1  ,  8.0  ,  15.0  ,,  2  .0 

0.1.  ,1.0  ,10.0 

4  500  ,  4400 , 4200 , 4000 , 3  700 , 3500 , 3200 , 2800 , 2400 ,  '  400 ,  .?  300 ,  1 000 ..  oW0 

4500,  4400 , 4200 , 4000 , 3700 , 3500 , 3200 , 2800 , 2400  ,  2400,  2300,  ir"  500 

1 000 ,  1000,1 000 , 1 000 , 1 000 , 1 000 , 1 000 ,  800 ,  800 ,  600 ,  600 ,  4 30 ,  4  50 

PH  I MAX , 1 
12 

0.0  ,0.8  ,1.0  ,1.4  ,2.0,  5.0,  5.1,  8-0  ,16.0  ,18.0  ,25.0  ,28.0 
0.0  ,0.0  ,0.0  ,0.0  ,0.0,  0,0,  0.0  ,  0.0,  0.O  ,  0.0  ,  2.0  ,  0.0 

CD0  „  1. 

15 

0.0,. 4  ,.8  ,.9  ,1.0  ,1.1  ,1.2  ,1.5  ,2.0  ,3.0  ,5.0  ,10.0  ,15.0  ,20.0  ,25.0 
0.0, ,009, ,009, .009, .023, .022, .018, .016, .013, .000, .004 , .002, .0019, .0017, .0015 
DELOD , 1 


0,0,50000. ,100000. ,150000. ,200000, 
0,0,  ,  0006  ,  .  002  ,  ,  00 1  ,  ..  0020 

CLALPHA , 1 
10 


0.0, 

1 . 5 

,2.0  , 

3.0  , 

4 . 0 

,  5.0 

, 6 . 0  ,8.0 

,15.0  ,25.0 

0.043, 

.043, 

.0391 , 

.  024 , 

..  020 

,  .018, 

,017, .016, 

.010  , .010 

»**  n 

10 

0 . 0  , 

1 . 5 

,2.0  , 

3.0  , 

4 . 0 

,  5.0 

, 6 . 0  ,8.0 

,15.0  ,25.0 

0 . 407 
A0AC ,  2 

,  0 , 407 

,  „  447 , 

.728, 

..  8  7  3 

,.970 

1 . 026 ,1.091 

,  1 . 745 , 1 . 745 

r , 

0.0  , 
-10,0  . 

0.01,  0 . 
,0.0  ,10.0 

1  , 

0 .,  5 

0.8  , 

1.0  , 

10.0  „  25.0 

,28.0 

80 . 0  , 

2 . 00 ,  0 

.  50 , 

0 . 1  5 

„ 

.  1  5  , 

.  !  6 

,  .51,  1.0 

,  3.  .  0 

80.0, 

2 . 00 ,  0 

.  50 , 

0 . 1.  5 

* 

.  1  5  , 

,  L  3 

,  .51,  1.0 

,  1 . 0 

80 . 0 , 
ISPR,  1 

2.00,  0 

.  50, 

0 ,  1  5 

*i 

.  .1  5  ,, 

1  s 

,  .51,  1.0 

,1.0 

30000 , 
440, 


80000 , 
445, 


390  ,  440,  445, 

NDOTPMAX , 1 
4 

0.0  ,  8.0  ,  8  ..01  , 

600  ,  600  ,  6(10  , 

1 00  0  ,  15 00,0 

150000  ,120000,  30000,0  ,27000 
2.0,  10. 0 

I  000 .0,1100.0, 900 . 0 ,  5.0  ,  5.0 

2.00,  1.00  ,8.0 

0 . 0292 

0  ,  I  000 

0  ,0  ,0 

0  ,0  ,0  ,0  ,0 

.1  ,.,  1  ,  .009  ,  .1  ,  2,0 


350000 

460 


50.0,900  „  38000 
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S.-2ST6.0UT 

01-29-1989 

20:17:47 

TIME 

H0 

V0 

Q0 

R 

WE16HT 

KDQTF 

108.00 

2167.31 

919.81 

943.25 

10.81 

395217.31 

27.84 

200.00 

10325.58 

923.11 

740.38 

25.79 

392033.19 

34.37 

300.00 

25111.56 

1474.98 

1155.32 

42.41 

386615.31 

87.87 

488.00 

48919.79 

2677.66 

1373.80 

76.95 

377626.69 

85.09 

508.00 

63469.52 

3311.88 

1048.32 

126.42 

370790.59 

58.73 

600.00 

63445.13 

3982.33 

1517.49 

185.95 

364502.31 

48.92 

708.00 

59209.52 

4018.83 

1892.10 

252.07 

359204.1? 

58.11 

900.00 

65138.59 

4133.29 

1507.68 

318.38 

353262.19 

76.99 

900.00 

73827.55 

4732.20 

1293.51 

391.68 

346526.59 

62.55 

1000.00 

71043.86 

5061.94 

1694.84 

473.45 

340969.69 

49,88 

1100.00 

73864.21 

5152.98 

1531.05 

556.93 

335650.50 

71.63 

1200.00 

81434.77 

5694.14 

1296.90 

646.38 

329224.69 

58.68 

1308.00 

78853.09 

6063.08 

1664.94 

744.13 

323766.00 

47.49 

1400.00 

81387.32 

6163.81 

1523.14 

844.11 

318647.00 

67.23 

1500.00 

88346.38 

6666.24 

1277.44 

949.98 

312657.1? 

54.91 

1600.00 

84914.17 

7033.21 

1674.76 

1063.51 

307359.59 

46.72 

1700.00 

87016.58 

7083.45 

1536.60 

1179.32 

302497.91 

62.01 

1900.00 

94202.17 

7557.83 

1244.34 

1299.98 

296783.19 

51.52 

1900.00 

90241.23 

7958.53 

1664.04 

1428.22 

291588.31 

45.45 

2000.00 

90952.39 

7961.97 

1610.27 

1559.16 

286874.41 

45.55 

2100.00 

98501.49 

8037.19 

1149.67 

1690.99 

145552.0? 

0.00 

2200.00 

187149.90 

9381.72 

1043.05 

1831.56 

129714.28 

0.00 

2300.00 

114917.40 

11303.33 

1044.63 

2001.82 

112933.60 

0.00 

2400.00 

121578.10 

13201.00 

1045.10 

2203.46 

98387.59 

0.00 

2500.00 

127428.40 

15083.61 

1045.54 

2436.22 

85585.84 

0.00 

2600.00 

132643.59 

16948.62 

1045 . 98 

2699,81 

74562.04 

0.00 

2700.00 

137355.70 

18797.00 

1046.30 

2993.96 

65021.18 

0.00 

2800.90 

141659.59 

20632.69 

1046.72 

3318.43 

56726.34 

0.00 

2900.00 

145592.50 

22444.70 

1047.52 

3672.94 

49525.4? 

8.00 

2000.00 

149177.20 

24219.52 

1048.93 

4856.97 

43287.58 

8.08 

3064.00 

151315.20 

25332.92 

1049.71 

4317.94 

39731.11 

0.00 

3083.80 

152095.20 

25757.05 

1850.59 

4422.21 

38426.70 

0.00 

UDOTA 

PHI 

HDOTP 

RES1 

3959.30 

0.24 

0.00 

0.00 

3096.59 

0.38 

0.00 

0.00 

3389.13 

0.89 

0.00 

0.00 

2914.02 

1.00 

0.00 

0.00 

2011.35 

1.00 

0.00 

0.00 

2693.12 

0.62 

0.00 

0.00 

3345.72 

0.59 

0.00 

0.00 

2636.56 

1.00 

0.00 

8.00 

2142.26 

1.00 

0.00 

0.00 

2744.32 

0.62 

0.00 

0.00 

2468.81 

1.00 

0.00 

0.00 

2089.63 

1.00 

0.00 

0.00 

2532.66 

0.64 

0.00 

0.00 

2302.55 

1.00 

0.00 

0.00 

1880.56 

1.00 

0.00 

0.00 

2432.43 

0.66 

0.00 

8.00 

2225.08 

8.95 

8.00 

0.00 

1764.35 

1.00 

0.00 

0.00 

2332.69 

0.67 

0.00 

0.00 

2256.25 

0.69 

0.00 

0.00 

0.00 

0.00 

48.61 

0.00 

0.00 

0.00 

180.50 

0.00 

0.00 

0.00 

155.76 

0.00 

0.00 

0.00 

136.83 

0.00 

0.00 

0.00 

118.53 

0.00 

0.00 

0.00 

101.71 

0.00 

0.00 

8.00 

88.34 

8.00 

0.00 

0.00 

76.86 

0.00 

8.00 

8.00 

66.85 

0.00 

0.00 

0.00 

57.79 

0.00 

0.00 

0.00 

53.43 

8.00 

0.00 

0.00 

57.71 

0.00 

45 


TIME 

H0 

V0 

!SP 

ISPAVG 

ISPEFF 

ISPMEAN 

THRUST 

THRUSTCOfl 

THRUSTHAX 

iee.ee 

2167.31 

919.81 

4267.90 

4366.57 

-8.81 

2395.35 

118334.78 

118834.70 

493417.81 

200.00 

10325.53 

923.11 

4256.79 

4326.29 

-7.06 

1484.12 

146289.70 

146289.78 

384901.41 

300.00 

25111.56 

1474.98 

4018.36 

4270,74 

2206.63 

1436.25 

353138.81 

353138.81 

397726.50 

400.00 

48919.79 

2677.66 

3340.07 

4022.78 

987.47 

1544.59 

284205.00 

388000.31 

284205.00 

500.00 

63469.52 

3311.83 

3031.27 

3824.78 

1114.59 

1435.55 

178031.30 

297046.31 

178031.30 

600,00 

63445.13 

3982.33 

2754.19 

3661.37 

1809.01 

1396.94 

134734.59 

134734.59 

216586.50 

700.00 

50209.52 

4018.33 

2739.10 

3541.64 

-1.64 

1225.65 

159171.50 

159171.58 

26759a. 91 

300.00 

65138.59 

4133.29 

2691.80 

3438.99 

954.57 

1096.76 

207235.59 

224654.30 

207235.59 

900.00 

73827.55 

4732.20 

2455.13 

3329.78 

885.60 

1075.43 

153573.00 

270389.81 

153578.00 

1000.00 

71043.86 

5061.94 

2396.17 

3242.90 

-1.14 

1035.78 

119516.50 

119516.50 

192015.00 

1100.00 

73864.21 

5152.98 

2393.30 

3172.89 

778.30 

963.33 

171421.41 

171421.41 

171972.00 

1200.00 

81434.77 

5694.14 

2375.17 

3101.29 

804.02 

950. '5 

139378.00 

255637.88 

139378.00 

1300.00 

78853.09 

6063.08 

2361.79 

3048.72 

-0.93 

333.30 

112153.40 

112153.48 

174663.30 

1400.00 

81337.32 

6163.81 

2358.61 

3005.47 

767.07 

885.57 

158579.70 

162870.00 

158579.78 

1500.00 

88346.38 

6666.24 

2342.01 

2960.55 

756.51 

877.86 

128605.48 

242208.30 

128605.40 

1600.00 

84914.17 

7033.21 

2328.56 

2924.71 

-8.80 

867.59 

103781.38 

108781.38 

165390.59 

1700.00 

87016.58 

7083.45 

2327.15 

2394.98 

606.75 

828.21 

144298.00 

144298.00 

151200.41 

1300.00 

94202.17 

7557.33 

2311.78 

2963.13 

714.66 

821.95 

119100.60 

229469.70 

119100.60 

1900.00 

90241.23 

7953.53 

2284.02 

2836.14 

-0.73 

818.16 

103804.70 

103804.70 

158641.00 

2080.00 

90952.39 

7961.97 

2284.09 

2813.14 

-0.71 

784.04 

104843.70 

104043.70 

153434.30 

210e.00 

98501.49 

8037.19 

446.03 

1267.74 

0.00 

343.79 

21682.32 

21632.32 

26761a. a9 

2200.00 

107149.90 

9381.72 

446.51 

1219.60 

403.64 

349.35 

80593.04 

88593.84 

267905.00 

2700.00 

114917.40 

11303.33 

446.94 

1174.43 

396.78 

352.20 

69616.70 

69616.78 

2a8163.c! 

2400.00 

121573.10 

13201.00 

447,31 

1139.17 

389.84 

354.1o 

61286.68 

61206.68 

268385. 5i 

2500.00 

127428.40 

15083.61 

447.63 

1111.13 

383.37 

355.60 

53056.02 

53056.02 

268580.91 

2600.00 

132643.59 

16943.62 

447.92 

1038.71 

390.65 

356.69 

45559.97 

45559.97 

268754.81 

2700.00 

137355.70 

18797.00 

448.19 

1070.4a 

339.59 

357.57 

39592.23 

39592.23 

263911.31 

2300.00 

141659.59 

20632.69 

448.43 

1055,.  43 

387.12 

358.25 

34465.09 

34465.09 

269055.31 

2900.00 

145592.50 

22444.70 

448.64 

1042.96 

382.40 

358.76 

29991.06 

29991.06 

269186.41 

7000.00 

149177.28 

24219.52 

448.84 

1032.57 

376.35 

359.11 

25940.01 

25940.81 

269305.91 

7064.00 

151315.20 

25332,92 

443.86 

1026.80 

368.53 

359.24 

23989.11 

23989.11 

269377.1° 

7088.80 

152095.28 

2575'. 85 

448,01 

1024. ’2 

326.0a 

359.27 

25911.53 

25911.53 

260403. 10 

46 


'I  HE 

H0 

V0 

02 

alpha 

100.00 

2167.31 

919.81 

943.25 

3.54 

200.03 

12325.53 

923.11 

742.38 

4.48 

300.00 

25111.56 

1474.93 

1155.32 

2.22 

400.00 

48919.79 

2677.66 

1373.80 

-5.00 

520.22 

63469.52 

3311.83 

1248.32 

-5.02 

402.20 

63445.13 

3932.33 

1517.49 

0.61 

700.22 

59209.52 

4018.83 

1892.12 

5.02 

322.22 

65133.59 

4133.29 

1507.68 

5.20 

922.22 

73327.55 

4732.20 

1293.51 

-5.00 

1202.22 

71043.86 

5061.94 

1694.84 

5.00 

1122.02 

73964.21 

5152.98 

1531.05 

‘.00 

1222.22 

31434.77 

5694.14 

1296.90 

-5.20 

1722.02 

73353.0? 

6263.28 

1664.94 

5.02 

1420.22 

31387.32 

6163.81 

1523.14 

5.02 

1522.22 

33346.78 

6666 . 24 

1277.44 

-5.00 

1622.22 

84914.17 

7033.21 

1674.7a 

5.00 

1720.00 

87016.53 

7033.45 

1536.62 

5.20 

1320.22 

94202.17 

7557.33 

1244.34 

-5.20 

1920.22 

90241.23 

7958.53 

1664.24 

5.00 

2222.20 

90952.39 

7961.97 

1610.27 

5.00 

2122.22 

98501.49 

6037.19 

1149.67 

5.00 

2222.22 

107149. 90 

93B1.72 

1043.05 

2.12 

2722.22 

114917.42 

11323.33 

1244.63 

1.39 

2422.22 

121578.10 

13201.20 

1245.10 

2.39 

2502.22 

127429.40 

15283.61 

1045.54 

2.26 

2602. 20 

132643.59 

16948.62 

1245.98 

1.49 

2702.22 

137355.70 

18797.20 

1246.32 

2.97 

7322.00 

141659.59 

22632.69 

104a. 72 

0.44 

2900.00 

145592.50 

22444.'0 

1247,52 

-0.14 

7222.22 

149177.20 

24219.52 

i043.93 

2.20 

3Es4.02 

151315.20 

25332.92 

1049.71 

0,87 

7238.82 

152095,22 

25757.25 

1050.59 

3.34 

GAMA 

LIFT 

DRAG 

3.3! 

792434.33 

95841.65 

6.94 

784066.38 

93597.94 

a,  22 

2.22 

125580.92 

4.65 

-1239892.02 

168502.59 

1.31 

-643245.00 

103528.78 

-1.18 

99953.24 

57373.41 

2.12 

1024311.02 

157792.30 

J  .47 

306796.13 

124345.50 

2.22 

-652093.00 

96261.4! 

-8.16 

825129.50 

119989.38 

8. 77 

745473.38 

109474.10 

0,24 

-612717.31 

90220.30 

-0.15 

'73596.38 

112590.32 

2.61 

70577a .6? 

103127.90 

2.13 

-533474.13 

85801.66 

-2.19 

755536. 88 

189393.10 

2.52 

692’38.31 

100652.48 

0.15 

-553352.38 

80989.87 

-2.32 

723801.88 

184898.00 

2.40 

725265.69 

101aS6.50 

2.58 

137266.90 

20138.63 

0.53 

2361.28 

6410.08 

2.3a 

39263.28 

7128.71 

2,27 

44215.3a 

7434.94 

2.2: 

'a?3?.23 

*777  30 

0.17 

2335’. 83 

5622.90 

2.14 

15158.83 

5024.11 

2.11 

6631.42 

4584.90 

2.12 

-2142.31 

4318.83 

0,06 

3122.2: 

4120.82 

0.07 

13643.94 

4274.48 

0.07 

52615.54 

7138.62 

rCCEL-S'S 

ACCCGH 

T/D 

-0.00 

0.88 

1.24 

-0.00 

0.00 

1.48 

0.50 

0.50 

3.34 

0.23 

0.50 

1.6? 

0.18 

0.50 

1 . 72 

0.31 

0.23 

2.35 

-0.88 

0.00 

1.01 

0.21 

8.2a 

1  .a’ 

0.16 

0.50 

1.70 

-0.00 

8.00 

1.80 

0.16 

0.P 

1.56 

0.14 

0.58 

1  ,:i 

-0.00 

0.00 

1,88 

0.16 

8.18 

1.54 

0.13 

0.50 

1.58 

-0.00 

0.00 

8.99 

0.12 

8.13 

1.43 

0.13 

0.58 

1 . 

-0.00 

8.00 

8, 95 

-0.00 

0.00 

1.82 

0.00 

8.00 

1.03 

0.57 

0.57 

12.57 

0.55 

0.55 

3.77 

0.55 

0.55 

8.23 

0.55 

0.54 

7 , 3- 

0.54 

0.54 

8.10 

0.54 

0.54 

7 1 30 

0.53 

0.53 

7,52 

8.53 

0.52 

6.94 

0.51 

8.51 

a, 29 

0.50 

0.50 

5.61 

0.49 

0.49 

7  -7 

47 


4.0  LISTING  OF  PROGRAM  ETO 


In  the  following  program  listing,  the  leftmost  column  of  numbers  are  sequential 
program  line  numbers  that  wore  generated  by  the  file  cross  referencing  code  used  to 
generate  the  program  listing.  They  are  not  part  of  the  nrogran.  These  line  ^umbers 
are  used  with  the  cross  referenced  listing  of  program  variables  in  Section  5.0. 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

If 

11 

12 

13 

14 

15 

16 

17 

18 
19 
2f 
21 
22 

23 

24 

25 

26 

27 

28 
29 
38 

31 

32 

33 

34 

35 

36 

37 

38 

39 
48 

41 

42 

43 

44 

45 

46 

47 

48 

49  468 
58 

51 

52  478 

53 

54 

55  488 

56 

57 

58 

59 
68 
61 
62 

63 

64 

65 

66 

67 

68 
69 

78  638 


1  1615/ 14JAN89  ET013.BAS  (Date  of  last  update) 


*  ET013.BAS 

*  DEVELOPMENT  VERSION  OF  "ETO. BAS" 


WRITTEN  BY: 

JOHN  L.  LE INGANG 
WAYNE  A.  DONALDSON 
KENNETH  A.  WATSON 
2LT  LOUIS  R.  CARREIRO 

AFWAL/POPA 

WPAFB.OHIO 

Written  during  period  January  1986  to  February  1987 

This  Is  a  2-DOF  trajectory  program  for  calculating  the  earth-to-orblt 
performance  of  single  stage  or  two  stage  horizontal  takeoff  vehicles. 
The  program  assumes  a  spherical  non-rotating  earth. A  2nd  order  Integ¬ 
ration  Is  used  for  the  equations  governing  welght.veloclty , flight  path 
angle, altitude, and  range. 

Rocket-only.alrbreathlng-only.and  simultaneous  rocket  and  airbreath¬ 
ing  propulsion  may  be  used  In  any  part  of  the  trajectory. 

Trajectory  options  are: 

*  Cruise  at  constant  Mach  number  and  altitude 

*  Cruise  at  constant  Mach  number  and  constant  L/D 

*  Climb  to  orbit  at  constant  Q 

*  Climb  to  orbit  along  a  velocity-biased  Q  path 

*  Climb  to  orbit  ball Istlcally, beginning  at  a  commanded  velocity 

Algorithms  for  table  loading  and  lookup  developed  by  Ken  Watson. 
Algorithm  for  constant  Q8  path  developed  by  Wayne  Donaldson. 


CLS 

WIDTH  LPRIHT  88 

DIN  I 1ST (58) ,  J1ST(58),  TI(58),  10(58),  AT ( 1000) 

DIM  TG(3)  .  GAO) 

OIM  DELW(4),  OELV ( 4 )  ,  D ELG AM ( 4 ) ,  DELH(4)  ,  DELR(4) 

DIM  RH8B( 18)  .  GLMBt 18)  ,  HB(1B),  TMB(IB) 

IWPtIT  "IS  THIS  A  SINGLE-STAGE!  1 )  .OR  A  2-STAGE  VEHICLE  (2)7  ".  NSTAGES 

IF  NSTAGES  •  2  GOTO  468 

IWPUT  "SINGLE  STAGE  FILE  NAME";  IFNI 

60T0  478 

IMPUT  "FIRST  STAGE  FILE  NAME";  IFNI 
IMPUT  "SECOND  STAGE  FILE  NAME";  IFN2I 

OPEN  "I".  #1,  IFNI:  IL  »  8:  IFILE  *  1 

GOTO  488 

LINE  INPUT  #1.  TITL1I 
LINE  INPUT  #1,  TITL2I 

I 

60SU8  16888:  'Input  alrbreather  isp  vs.  Mach  and  PHI  as  table  #1 

GOSUB  16888:  'Input  alrbreather  PHImax  vs.  Mach  as  table  #2 

GOSUB  16888:  'Input  Cdo  vs. Mach  as  table  13 

GOSUB  16888:  'Input  delta  Cd  vs.  Alt.  as  table  #4 

GOSUB  16888:  'Input  Clalpha  vs.  Mach  as  table  #5 

GOSUB  16888:  'Input  K  vs.  Mach  as  table  46 

GOSUB  16888:  'Input  alrbreather  Ao/Ac  vs.  Mach  and  alpha  as  table  17 

GOSUB  16886:  'Input  rocket  Isp  vs.  Altitude  as  table  48 

GOSUB  16688:  ‘Input  rocket  WOOTPMAX  vs. Mach  as  table  49 

INPUT  "Do  you  want  a  listing  of  tabular  Input  data(Y/N)";  WTIN$ 

IF  WTINt  ■  "Y"  OR  WT INI  •  "y"  THEN  638  ELSE  648 
FOR  ITAB  »  1  TO  9:  GOSUB  16238:  NEXT  ITAB 


71  648 
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72  INPUT  "Oo  you  want  printer  output  during  run(Y/N)";  PRCOMS 

73  IF  PRC0M$  »  "Y"  OR  PRCOM*  *  "y"  THEN  PRCOM  =  1  ELSE  PRCOM  =  0 

74  IF  PRCOM  =  1  THEN  SOTO  680  ELSE  SOTO  700 

75  680  '"  LPRINT  CHR$(27);  "B";  CHR$ ( 3) ;  :  'Select  condensed  print  (SR15  printer) 

76  LPRINT  CHR$( 27) ;  CHR$(33);  CHR*(4);  :  'Select  condensed  pr1nt(Epson  printer) 

77  GOTO  710 

78  700  "*  LPRINT  CHR${27) ;  "B";  CHRS(l):  :  'Select  standard  print  (SR15  printer) 

79  "LPRINT  CHRS ( 27 ) ;  CHR${33);  CHR$(0) ;  :  'Select  standard  print  (Epson  printer) 

•  1  •  »  »  •  I  I  I  t  ■  »  i  #  l  I  I  I  l  I  I  I  l  #  I  I  •  l  l  •  i  i  I  i  i  »  l  •  l  l  l  •  I  I  i  *  l  «  l  l  I  »  l  i  i  l  l  »  l  l  i  t  i  t  l  l  l  l  i  t  t  l 

81  710  INPUT  #1,  AC,  AREF 

82  INPUT  #1,  VLAUNCH,  WFUEL ,  WFINAL,  VFINAL ,  STAGE 

83  INPUT  01,  M0,  H0,  GAMMA,  ALPHA,  DT.  OELPRINT ,  TIMEX 

84  INPUT  #1,  Q0COMD ,  Q0MAX  ,  Q0M1N ,  ALPHAMAX,  ALPHAMIN 

85  INPUT  01,  VTAKEOFF ,  ALPHA2MAX 

86  INPUT  01,  LOAOFAC ,  ACCCOMD,  GAMMAMAX 

87  INPUT  01,  FASTOIC 

88  INPUT  01,  V0TEMP,  Q0FINAL 

89  INPUT  01,  V0CRUISE  ,  SWITCH4 ,  GAMMA4MAX 

90  INPUT  01,  FLAG5 ,  V05,  ALPHA5MAX  ,  GAMMA 5 ,  ACCC0MD5 

91  INPUT  01,  Al,  A2,  Bl,  Cl,  C2 ,  01.  DELQ,  V02,  H02 

92 

93  PRINT 

94  PRINT  "1ST/SINGLE  STAGE  INPUT  FILE  HAS  BEEN  READ." 

95  PRINT  "LAST  ITEM  IN  INPUT  FILE  IS  H02, VALUE  READ  FROM  FILE  IS";  H02 

96  PRINT  "IF  THIS  VALUE  IS  INCORRECT , INPUT  FILE  IS  IN  ERROR." 

97  PRINT 

98  INPUT  "Do  you  want  a  plotter/printer  file  created?  (Y/N)";  PLOTS 

99  IF  PLOTS  *  "N"  OR  PLOTS  ■  "n"  THEN  SOTO  1020 

100  INPUT  "Enter  the  plotter/printer  file  name  (DRIVE: NAME .EXT)  ",  PFNS 

101  OPEN  "0",  #2,  PFNS 

102  VARNUM  =  25 

103  WRITE  02,  VARNUM 

104  PRINT 

105  PRINT 

106  PRINT  "At  the  prompt  below, enter  a  value  for  FREQ, the  Interval  of  plotter/" 

107  PRINT  "printer  output." 

108  PRINT  “A  value  of  1  for  FREQ  means  output  at  each  time  increment  * 

109  PRINT  "A  value  of  2  means  every  other  time  increment, etc.  As  a  guide  to  " 

110  PRINT  "the  amount  of  disk  space  required, use  the  following:" 

111  PRINT 

112  PRINT  "  BYTES  REQUIRED  *  (200)  X  (FLIGHT  TIME)" 

113  PRINT  "  . " 

114  PRINT  "  (INTEGRATION  INTERVAL)  X  (FREQ)" 

115  PRINT 

116  PRINT  "  Very  large  files  exceeding  200X  may  be  generated. Exit  the  * 

117  PRINT  "  program  by  pressing  CTRL  C  If  disk  space  Is  insufficient," 

118  PRINT  "  otherwise  enter  value  for  FREQ  below:" 

119  INPUT  "  ",  FREQ 

120  1020  STARTTIMES  *  TINES:  START  *  TIMER 

121  IF  EOF(l)  THEN  GOTO  1022 

122  1022  CLOSE  01 

123 

124  '  Use  this  area  to  the  line  below  as  a  "scratch"  space  to 

125  '  temporarily  overlde  values  read  from  the  lst/slngle  stage 

126  '  Input  file. 

127 

128 

129  .  OT  *  .1 

130  "  "  OELPRINT  *  1 


131 

132 

133  '  . 

134  CLS  :  'Clears  screen 

135 

136  ' .  EVALUATE  CONSTANTS  . 

-137 

138  GC  *  32.174 

139  MEARTH  »  (5.98E*27)  /  454:  'Earth  mass 

140  REARTH  *  (6.38E+08)  /  (2.54  *  12):  'Earth  radius 

141  BIGG  *  (6.67 E- 11)  *  (1  /  4.4482)  *  (.4536  /  .3048)  "  2:  'Universal  gravitation  const 

141  . 

142  FC  •  (WEIGHT  /  GC)  *  (V0)  "  2  /  (H0  ♦  REARTH):  'Centrifugal  force 
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143  ANG  ■  57.3:  'Degrees  per  radian 

144 

145 

146  ' .  INITIALIZE  THE  PROBLEM  . 

147 

148  IF  TIMEX  *  6!  THE*  WEIGHT  «  WLAUNCH 

149  COUNT  «  0 

150  PRCOUNT  «  0 

151  1-0 

152  ' - Load  data  for  atmosphere  routine  - 

153  FOR  J1  -  1  TO  10 

154 

155  READ  HB(J1) :  READ  RH0B(J1):  READ  TMB(Jl):  READ  GLMB(Jl) 

156  REIT  J1 

157 

158  2540  DATA  0,2. 3769E-3, 518. 7  ,-3. 5662E-3 

159  DATA  36089, 7. 0612E-4, 390,0 

160  DATA  65617  ,1 .7082E-4 ,390  ,5 .4864E-4 

161  DATA  104986,2  . 566 IE-5 ,411  .6,1 .5362E-3 

162  DATA  154199, 2. 7698E-6. 487. 2,0 

163  DATA  170604. 1.4735E-6, 487. 2. -1.0973E-3 

164  DATA  200131  ,4.87 19E-7  ,454.8  ,-2 .1946E-3 

165  DATA  259186, 3. 8826E-8, 325. 2,0 

166  DATA  291153,6 . 1508E-9 ,325 .2 ,1 .6953E-3 

167  DATA  323003,9 .65 11E- 10 ,374.2  ,2  .8343E-3 

168  RESTORE  2540 

169 

170  ' $PAGE 
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171  .•*••*«*.*»***•*••*••***•••***»**••*.*.*.******••»****•*»**•»*••»*»*•»••**• 

172  '  MAIN  PROGRAM  ROUTINE 

173  .**••*•*•*••*••»•**•••••»•*»*.**•••••*••••••****•••*•**•*••••*•**•*•*»**** 

174 

175  ' - Get  current  values  from  input  data: 

176  ITAB  *  3:  XARG  *  MB:  GOSUB  164BB:  COB  <=  VALI:  'Look  up  CdB 

177  ITAB  *  4:  XARG  *  HB:  GOSUB  164BB:  DELCD  *  VALI:  ‘Look  up  delCd  friction 

178  ITAB  *  5:  XARG  =  MB:  GOSUB  164BB:  CLALPHA  *  VALI:  'Look  up  Clalpha 

179  ITAB  *  6:  XARG  *  MB:  GOSUB  164BB:  K  ■=  VALI:  'LOOK  UP  K 

18B  ITAB  *  7:  XARG  *  MB:  YARG  *  ALPHA:  GOSUB  164BB:  ABAC  *  VALI:  'LOOK  UP  ABAC 

181 

182  ' .  CALCULATE  CURRENT  LIFT  AND  DRAG  . 

183 

184  CL  «  ALPHA  *  CLALPHA 

185  LIFT  >  CL  *  OB  *  AREF 

186  CD  *  COB  ♦  K  *  CL  *  2  +  DELCD 

187  DRAG  *  CD  *  QB  *  AREF 

188 
189 

19B  FOR  1=8  TO  1BBBBBB! 

jg J  I  I  »  I  I  « 

192  31BB  AT  ■  IRKEYT:  IF  AT  =  CHR$(3)  THER  GOTO  1326B:  'Intercepts  CTRL-C  entered 

193  1  from  keyboard  to  halt  execution  and  exit  program 

194  QT  *  IRKEYT:  IF  Q!  *  CMR$(19)  THER  GOTO  31B8  ELSE  GOTO  3118:  ‘Intercepts 

195  '  CTRL-S  to  halt  execution  until  space  bar  Is  pressed 

196  31B8  ' 

197  31 IB  QT  -  IRKETT 

198  IF  QT  <>  CHR${ 32 )  THER  GOTO  31  IB :  'Resumes  If  space  bar  Is  pressed 

199 

28B  3118  ' 

281  GOSUB  135BB:  '  Go  to  atmosphere  subroutine  and  return  here 

2B2 

2B3  GA  *  BIGG  *  MEARTH  /  (HB  +  REARTH)  *  2  *  GC:  'Gravitational  acceleration 

2B4  '- . Reduce  Integration  Interval  near  end  of  flight 

2B5  IF  WEIGHT  -  WFINAL  <  .81  *  WF INAL  THER  DT  *  .1 

286  IF  VB  <  24BBB  THER  GOTO  3788 

287  DV  *  CIHT( ( 262BB  -  Vfl)  /  1BBB) 

2B8  DT  *  OV 

2B9  IF  DT  *  B  THER  DT  =  .1 

2IB  3788  ' 

211  GOTO  49BB 

212  49BB  ' 

213 

214  ' .  FLIGHT  PHASE  COMMANDS  . 

215 

216  PHASE  »  1 

217  '  Takeoff  Roll  at  zero  alpha  until  VTAKEOFF  Is  reached 

218  IF  VB  >*  VTAKEOFF  OR  TIMEX  >  T TAKEOFF  THER  65BB:  'Test  to  see  If 

219  '  Phase  1  Is  completed 

228  ALPHA  -  8!:  GAMMA  *  B! :  LIFT  *  WEIGHT:  HB  *  8! 

221  ACCCOM  -  ACCCOMO 

222  GOSUB  14868:  'Go  to  thrust  subroutine  and  return  here 

223  ALPHA  =■  8!  :  GAMMA  =>  8!:  LIFT  *  WEIGHT:  HB  *  B! 

224 

225 

226  TTAKEOFF  *  TTAKEOFF  ♦  OT:  AT  *  "PHASE  l, TAKEOFF  ROLL  **•****•***♦*- 

227  GOSUB  15866:  '  Go  to  Integration  subroutine  then  return  here 

228  GOTO  12886:  '  Go  to  runtime  output  section 

229  6588  ' . 

238  PHASE  *  2:  'Commanded  climb  at  constant  load  factor.  Required  ALPHA  Is 

231  *  calculated, but  Is  limited  to  ALPHAMAX.  Maximum  velocity  allowed 

232  '  Is  V62  fps.  ALPHA  Is  commanded  to  zero  above  HB2  ft.  This  phase 

233  '  Is  maintained  until  V82  fps,H82  ft, and  6.95*Q6C0M0  are  reached. 

234 

235  IF  TIMEX  >  TTAKEOFF  ♦  TQfl  ♦  DT  THER  6868:  'Test  If  Phase  2  Is 

236  '  completed 

237  IF  V8  >  VB2  ARD  HB  >  H62  AHD  Q6  >  .95  *  Q6C0MD  THER  GOTO  6866 

238  ALPHA  -  WEIGHT  •  LOAOFAC  /  (CLALPHA  *  Q 6  *  AREF) 

239  IF  ALPHA  >  ALPHA2MAX  THER  ALPHA  »  ALPHA2MAX 

248  ACCCOM  •  ACCCOMD 

241  IF  V8  >  V82  AHD  H6  <  HB2  THER  ACCCOM  <  6 
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242 

243 

244 

245 

246 

247 

248 

249 

250  6860 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264  7100 

265 

266 

267 

268  7180 

269 

270  7220 

271 

272 

273 

274 

275  7320 

276 

277 

278  7380 

279 

280  7420 

281 
282 

283 

284  7500 

285 

286 

287 

288 

289 

290 

291 

292  7660 

293 

294 

295 

296 

297 

298  7780 

299 

300 

301 

302 

303 

304 

305 

306 

307 

308  7980 

309 

310 

311 

312 

313 


IF  H0  >  H02  THEN  ALPHA  *  0 
TQ0  *  TQ0  ♦  OT 

B$  *  ‘PHASE  2, CLIMB  AT  CONSTANT  l/W***************" 

GOSUB  14000:  ‘Go  to  thrust  subroutine  and  return  here 

GOSUB  15800:  'Go  to  Integration  subroutine  then  return  here 
GOTO  12000:  'Go  to  runtime  output  section 


PHASE  =  3:  'Commanded  climb  at  constant  Q0.  An  angle  of  attack  (ALPHA) 
Is  computed  that  will  correct  the  flight  path  angle  (GAMMA)  to 
seek  the  commanded  Q0COM. 

This  is  GUID13  version.  29SEPT86.  THIS  ALGORITHM  DEVELOPED  BY 
WAYNE  DONALDSON 


Q0COM  *  Q0COMD 

IF  V0  <  V0TEMP  THEN  GOTO  7100 

IF  FLAG5  >  0  AND  TSET5  *  1  THEN  GOTO  8300 
Q0COM  »  Q0COMD  -  (Q0COMD  -  Q0FINAL)  *  (V0  -  V0TEMP)  /  (26200  -  V0TEMP) 

Q0COMD  Is  reduced  for  V0  >  V0TEMP 

IF  Q0  <  Q0COM  THEN  GOTO  7320 
IF  Q0  <  QOLO  THEN  GOTO  7220 

RATEOQ  -  (Q0  -  QOLO)  /  DT 
ALPHA  =  ALPHA  +  A1  +  RATEOQ  *  A2 

GOTO  7500 

RATEDELQ  -  (Q0  -  Q0COM)  /  DT 

IF  Q 0  -  Q0COM  >  DELQ  THEN  RATEOQ  *  (QOLD  -  Q0)  /  DT 
IF  Q0  -  Q0COM  >  DELQ  THEN  GOTO  7180 
ALPHA  »  ALPHA  -  RATEDELQ  *  B1 

GOTO  7500 

IF  Q0  >  QOLO  THEN  GOTO  7420 

RATEDQ  ■  (Q0  -  QOLO)  /  OT 
ALPHA  =■  ALPHA  -  Cl  ♦  RATEDQ  *  C2 

GOTO  7500 

RATEDELQ  =  (Q0  -  Q0COM)  /  OT 

IF  Q0COM  *  Q0  >  DELQ  THEN  RATEDQ  «  (QOLD  -  Q0)  /  DT 
IF  Q0COM  -  Q0  >  OELQ  THEN  GOTO  7380 
ALPHA  =  ALPHA  ♦  RATEDELQ  *  D1 
QOLD  =  Q0 

IF  ALPHA  >  ALPHAMAX  THEN  ALPHA  *  ALPHAMAX 
IF  ALPHA  <  ALPHAMIN  THEN  ALPHA  *  ALPHAMIN 


CALCULATE  ACCELERATION  COMMAND  FOR  PHASE  3 


IF  Q0  <*  Q0COM  GOTO  7660  ELSE  GOTO  7780 
""""'For  Q0  <  Q0COM ,  maximum  available  acceleration  Is  commanded 

I 

ACCCOM  =  ACCCOMD 

GOTO  7980 


""""For  Q0  >  Q0MAX ,  acceleration  Is  reduced  linearly  to  zero  for  Q0 

between  Q0COM  and  Q0MAX  and  Is  held  to  zero  when 
Q0MAX  Is  exceeded 

ACCPERQ  *  ACCCOMD  /  (Q0COM  -  Q0MAX) 

ACCCOM  *  ACCCOMD  -  ACCPERQ  *  (Q0COM  -  Q0) 

IF  Q0  >  Q0MAX  THEN  ACCCOM  •  0 
GOTO  7980 


I  I  I  I  *  I  I  i 


Q0OLO  *  Q0 
RH0OLO  •  RH0 
H0OLD  »  H0 

IF  V0CRUISE  >  0  AND  V0  >•  V0CRUISE  THEN  GOTO  8180 
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314  IF  FLAGS  >  0  AND  V0  >  VBS  THEN  GOTO  83SB 

315  C$  =  "PHASE  3. CLIMB  AT  COMMANDED  Q****************" 

316  G0SU8  14BBB:  'Go  to  THRUST  subroutine  and  return  here 

317  60SUB  1580B:  'Goto  integration  subroutine  and  return  here 

318  GOTO  12000:  'Go  to  runtime  output  section 

319 

320  *  . 

321  8180  PHASE  =  4:  'Cruise  at  commanded  V0CRUISE  and  Q0COM 

322  ACCCOM  *  0 

323  IF  SWITCH4  =  0  THEN  GOTO  8260 

324  IF  FLAGCRUISE  «  1  THEN  GOTO  8250 

325 

326  FOR  ALPHA  *  0  TO  ALPHAMAX  STEP  .1 

327  ITAB  -  3:  XARG  -  M0:  GO SUB  16400:  COB  *  VALI:  'Look  up  CdB 

328  ITAB  «  4:  XARG  =  H0:  60SUB  16400:  DELCD  *  VALI:  'Look  up  delCd  friction 

329  ITAB  =  5:  XARG  *  MB:  GOSUB  16400:  CLALPHA  *  VALI:  'Look  up  Clalpha 

330  ITAB  =  6:  XARG  =  MB:  GOSUB  16400:  K  *  VALI:  'LOOK  UP  K 

331  CL  =  ALPHA  *  CLALPHA 

332  CD  -  CDB  +  K  *  CL  *  2  ♦  OELCD 

333  LOVRD  =  CL  /  CO:  PRINT  ALPHA;  LOVRD 

334  IF  LOVRD  <  LOVRDOLD  THEN  GOTO  8225 

335  LOVRDOLO  *  LOVRD 

336  NEXT  ALPHA 

337 

338  8225  ' . 

339  ALPHACRUISE  »  ALPHA  -  .1 

340 

341  FLAGCRUISE  =  1 

342  8250  ALPHA  =  ALPHACRUISE 

343  IF  GAMMA  >  GAMMA4MAX  THEN  ALPHA  =  0! 

344  IF  LIFT  <  WEIGHT  THEN  ALPHA  =  0! 

345  IF  GAMMA  <  0!  THEN  ALPHA  «  ALPHAMAX 

346 

347  8260  ' . 

348  D$  *  "PHASE  4, CRUISE  AT  CONSTANT  ALTITUDE  ************************ 

349  IF  SWITCH4  »  1  THEN  Dl  *  "PHASE  4, CRUISE  AT  CONSTANT  L/D  ************** 

350  GOSUB  14000:  'Go  to  thrust  subroutine  and  return  here 

351  GOSUB  15800:  'Goto  integration  subroutine  and  return  here 

352  GOTO  12000:  'Go  to  runtime  output  section 

353  8300  ' . 

354  PHASE  =5:  'At  V05  begin  pullup  to  flight  path  angle  GAMMA5  then 

355  '  continue  bal 1 1st  leal ly  at  commanded  acceleration  ACCC0M05 

356 

357  ACCCOM  =  ACCC0MD5 

358  ALPHA  =  ALPHA5MAX :  ALPHAMAX  »  ALPHA5MAX 

359  IF  TSET5  *  1  THEN  ALPHA  *  0 

360  IF  TSET5  *  1  THEN  GOTO  8500 

361  IF  GAMMA  >«  GAMMA5  THEN  TSET5  *  1  ELSE  TSET5  *  0 

362 

363  8500  E$  =  “PHASE  5, PULLUP  AND  BALLISTIC  ASCENT" 

364  GOSUB  14000:  'Go  to  thrust  subroutine  and  return  here 

365  GOSUB  15800:  'Go  to  integration  subroutine 

366  GOTO  12000:  'Go  to  runtime  output  section 

367 

368  12000  ' 

369  ' .  RUNTIME  OUTPUT  SECTION  . 

370 

371 

372  ISPEFF  «  WEIGHT  *  (V0  -  V00)  /  <GC  *  (WEIGHT0  -  WEIGHT)) 

373  ITCUM  *  ITCUM  ♦  ISPEFF  *  (WEIGHT!  -  WEIGHT) 

374  ISPMEAN  *  ITCUM  /  (WLAUNCH  -  WEIGHT) 

375 

376  ITSUM  -  ITSUM  +  ISP  *  (WEIGHT!  -  WEIGHT) 

377  ISPAVG  »  ITSUM  /  (WLAUNCH  -  WEIGHT) 

378  ACC  -  (V0  -  V00)  /  (DT  *  GA) 

379 

380 

381  'First  test  for  weight  or  velocity  stopping  condition: 

382  VORBIT  -  REARTH  *  SqR(GC  /  (REARTH  +  H0)) 

383  IF  V0  >*  VORBIT  GOTO  13020 

384  IF  WEIGHT  <«  WF INAL  GOTO  13140 

385  IF  V0  >-  VFINAL  AMD  STAGE  *  1  THEN  GOTO  13190 
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WDOTA*#### ,## 


PHI-##.### 


WT-#######.# 


HDOTP-## 


RANGE  /  617 


L/W  -### 


FC-####### .#" ;  ISPEFF; 


LIFT-#######.# 


L/D  *### 


DRAG-#######.#";  ISPAVG ; 


THRUST-#######.# 


T/D  -### 


THRUSTCOM-#######.#*;  ISPR; 


THRUSTMAX-#######.# 


T/W  -### 


386 

387 

388  COUNT  «  COUNT  +  1 

389  IF  COUNT  <»  (OELPRINT  -  1)  THE*  GOTO  12980 

390  CLS 

391  IF  PHASE  »  1  THE*  X$  -  A$  ELSE  IF  PHASE  -  2  THE*  XI  -  B$  ELSE  IF  PHASE  -  3  THE* 

391  X$  =  C$  ELSE  IF  PHASE  -  4  THE*  X$  »  0$  ELSE  IF  PHASE  -  5  THE*  X$  -  E$ 

392  PRIHT  USING  "T-####.#  \  \* ;  TIMEX;  X$ 

393  PRINT  USING  "V0-#####.#  M0-##.#  ALPHA-###.##  WOOTF-#### .##  RES1-### 

393  ##.#";  V0;  M0;  ALPHA;  WOOTF;  RES1 

394  PRINT  USIN6  "H0-######.#  GAMMA-###.##  WDOTA-####.##  WDOTP-# 

394  ##.#";  H0;  GAMMA;  WDOTA;  WOOTP 

395  PRINT  USING  "  R*######.##  00-####.##  PHI-##.###  ";  RANGE  /  60 

395  6;  00;  PHI 

396  PRINT  USING  "ISP-#####.##  ALPHAMAX-### .##  WT-#######.#  L/W  -#fi 

396  .###";  ISP;  ALPHAMAX;  WEIGHT;  LIFT  /  WEIGHT 

397  PRINT  USING  "ISPEFF-####.#  ALPHAMIN-### .##  FC-#######.#“;  ISPEFF; 

397  ALPHAMIN;  FC 

398  PRINT  USING  "ISPMEAN-####.#  ACCCOM  =###.##  LIFT-#######.#  L/D  -## 

398  .###";  ISPMEAN;  ACCCOM;  LIFT;  LIFT  /  DRAG 

399  PRINT  USING  “ISPAVG-#####.#  ACCEL-G’S-###.##  DRAG-#######.#";  ISPAVG; 

399  ACC;  DRAG 

400  PRINT  USING  "ISPA/B-#####.#  THRUST-#######.#  T/D  -## 

400  .###";  ISPA;  THRUST;  THRUST  /  DRAG 

401  PRINT  USING  "I SPROC-##### .#  THRUSTCOM-####### .#*;  ISPR; 

401  THRUSTCOM 

402  PRINT  USING  "  THRUSTMAX-#######.#  T/W  «## 

402  .###’;  THRUSTMAX;  THRUST  /  WEIGHT 

403 

404  COUNT  =  0! 

405 

406 

407  IF  PRCOM  <>  1  GOTO  12980 

408 

409  IF  PRCOUNT  <  I  GOTO  12660 

410  IF  2  <-  PRCOUNT  <=  14  GOTO  12780 

411  12660  LPRINT  CHR$(12) 

412  IF  PRCOUNT  <  1  THEN  LPRINT  TAB( S) ;  TITL1S,  TABC 5) ;  TITL2S ,  OATE) ,  TINES 

413  WIDTH  LPRINT  136 

414  LPRINT  TAB(5) ;  “TIME";  TAB( 15) ;  "V0";  TAB(25);  "H0";  TAB(35);  "R";  TABC45);  "ISP"; 
414  TA8( 55) ;  “ISPEFF";  TABC 65) ;  "ISPMEAN";  TABC 75) ;  "WDOTF" ;  TAB(85);  “WDOTA";  TAB(97);  "PHI 

414  ;  TABC  105)  ;  "THRUST";  TABU  12)  ;  "THRUSTCOM";  TAB(122);  "THRUSTMAX" 

415  LPRINT  TABC 15) ;  "00";  TABC 25) ;  "ALPHA";  TABC35);  "GAMMA";  TAB(45);  "WEIGHT";  TAB(5 
415  );  "LIFT";  TABC  65) ;  "DRAG";  TABC  75) ;  "WDOTP" ;  TAB(85);  "RES1";  TAB(94);  "ISPAVG";  TABC10 

415  );  "ACCEL-G'S";  TABC115) ;  "ACCCOM";  TABC 128);  "T/D" 

416 

417  12780  LPRINT 

418  PRCOUNT  -  PRCOUNT  +  1 

419  WIDTH  LPRINT  136 

420 

421  LPRINT  USING  "#######.##“;  TIMEX;  V0;  H0 ;  RANGE  /  6076;  ISP;  ISPEFF;  ISPMEAN; 

421  WDOTF;  WDOTA;  PHI;  THRUST;  THRUSTCOM;  THRUSTMAX 

422  LPRINT  USING  "#######.##";  TIMEX;  Q0;  ALPHA;  GAMMA;  WEIGHT;  LIFT;  DRAG;  WDOTP; 

422  RES1 ;  ISPAVG;  ACC;  ACCCOM;  THRUST  /  DRAG 

423  WIDTH  LPRINT  80 

424 

425  IF  PRCOUNT  >»  16  THEN  PRCOUNT  -  0 

426 

427  12980  NEXT  I 

428 

429 


"GAMMA";  TABC45) ;  "WEIGHT";  TABC55 
"RES1* ;  TABC 94) ;  "ISPAVG";  TABC102 
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430 

431 

432 

433 

433 

434 

434 

435 

436 

436 

437 

437 

438 

439 

439 

440 

440 

441 

442 

443 

444 

445 

446 

447 

448 

449 

450 

451 

452 

453 

454 

455 

455 

456 

457 

458 

459 

460 

461 

462 

463 

464 

465 

466 

467 

468 


13020  * . TEST  FOR  STOPPING  CONDITIONS 

I 


PRINT  USING  "REACHED  ORBITAL  VELOCITY*!#####.#  FPS  AT  TIME=#####.#SECS.  .HEIGHT* 
######.#  L8S.“ ;  V0;  TIMEX;  HEIGHT 

IF  P RCOM  *  1  THEN  LPRINT  USING  "REACHED  ORBITAL  VELOCITY*######.#  FPS  AT  TIME*## 
###.#SECS. .HEIGHT*  ######.#  LBS.";  V0;  TIMEX;  HEIGHT 

GOTO  13220 

13140  PRINT  USING  "FUEL  EXHAUSTED  AT  TIME*  ##### .#SECS. .HEIGHT*  ######.# .VELOCITY*#### 
##.#  FPS.";  TIMEX;  HEIGHT;  V0 

IF  PRCOM  =  1  THEN  LPRINT  USING  "FUEL  EXHAUSTED  AT  TIME*  #####  .#SECS.  .HEIGHT*  ### 
###.#, VELOCITY*######. #FPS.“;  TIMEX;  HEIGHT;  V0 
60T0  13220 

13190  PRINT  USING  "REACHED  1ST  STAGE  VELOCITY*######.#  FPS  AT  TIME*##### .#SECS .  .HEIGH 
T-  ######.#  LBS.";  V0 ;  TIMEX;  HEIGHT 

IF  PRCOM  =  1  THEN  LPRINT  USING  "REACHED  1ST  STAGE  VELOCITY*######.#  FPS  AT  TIME* 
#####. #SECS., HEIGHT*  ######.#  LBS.";  V0;  TIMEX;  HEIGHT 
GOTO  13220 


13220  IF  NSTAGES  *  1  THEN  GOTO  13250 

13240  IF  NSTAGES  *  2  AND  STAGE  *  1  THEN  GOSUB  18000:  ‘Reads  2nd  stage  data 

IF  TIMEX  -  TIMEX2  THEN  GOTO  3100 

13250  IF  PRCOM  *  0  AND  PLOT!  *  "N"  OR  PLOTS  *  "n"  THEN  GOTO  13300 

I  l 

13260  "LPRINT  CHR$(27);  "B";  CHR$(1);  :  'Select  standard  print  (SR15  printer) 

LPRINT  CHR$(27);  CHR$(33);  CHR$(0);  :  ‘Select  standard  print  (Epson  printer) 

UIOTH  LPRINT  80 

IF  PLOTS  *  "N“  OR  PLOTS  =  "n"  THEN  GOTO  13300 
INPUT  "Do  you  want  printer  output  of  run(Y/N)“;  PLOTS 
IF  PLOTS  *  “N“  OR  PLOTS  =  "n"  THEN  GOTO  13300 
CLOSE  2 

IF  PLOTS  *  "Y"  OR  PLOTS  *  "y"  THEN  GOSUB  20000:  ‘Goes  to  printer  output  subroutine 

i 3300 

PRINT  TINES 

PRINT  STARTTIMES:  FINISH  *  TIMER 
PRINT  "PROGRAM  TOOK"  ;  FINISH  -  START; 

PRINT  "SECONDS" 

END 

I 

' - End  of  main  routine  - 
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469  ' 

478  ' ************************************•***************< 

471 

472  '  SUBROUTINES 

473 

474  . ft***************************************************' 

475  .*•****••••**»*•••**•**»»•****»•••**••*•*•**•*****••*< 

476 

477  13566  ***************************************************** 

478  '  SUBROUTINE  TO  DEFINE  ATMOSPHERIC  CONDITIONS 

479  '  USING 

486  '  U.S.  STANOARD  ATMOSPHERE,  1962 

481  '  4FEB87  VERSION 

482  .  *••****•***•****••••**••*••••*****»•***•*•*•••*••*•*< 

483 

484  ‘ . DATA . 

485  'FOR  J1  =  1  TO  16 

486  'READ  HB( J1 ) : READ  RH6  B ( J 1 )  :  RE  AD  TMB(Jl) : READ  GLMB(Jl) 

487  '  NEXT  J1 

488 

489  13546  'DATA  0 ,2 .3769E-3 ,518  .7  .-3.5662E-3 

496  'DATA  36689,7 .6612E-4, 396, 6 

491  'DATA  65617 ,1 . 7082E-4 ,396  .5.4864E-4 

492  'DATA  164986 .2 .5661E-5 .411 .6 .1 .5362E-3 

493  'DATA  154199  ,2 . 7698E-6 , 487  .2 .0 

494  'DATA  170664 ,1 .4735E-6 ,487 . 2  . -1 .697 3E-3 

495  'DATA  266131  ,4 .87 19E-7  ,454 .8  , -2 . 1946E-3 

496  'DATA  259186 ,3. 8826E-8  ,325 .2,0 

497  'DATA  291 153 ,6 . 1508E-9 ,325 .2 . 1 .6953E-3 

498  ‘DATA  323083,9  .6511E-10  ,374.2  .2 .8343E-3 


499  '  RESTORE  13540 

580  ' . SET  CONSTANTS . 

501  N1  -  16 

582  GC  *  32.17405 

503  RE  =  2 .6899E+07 :  '  radius  of  earth 

504  CPCV  =  1.4:  '  ratio  of  specific  heats 

585  '  -  -  -  -  CONVERT  TO  GEOPOTENTIAL  ALTITUDE  -  -  -  - 

586  H  =  (RE  *  H0 )  /  (RE  *86!:  '  geopotential  altitude 

507  ' .  FINDS  ALTITUOE  BASE  POINT  . 


588  FOR  II  -  2  TO  N1 

509  HDEL  =  H  -  HB(I1) 

510  IF  HDEL  <*  0  SOTO  13740 

511  NEXT  II 

512 

513 

514  13746  ' . CALCULATE  TEMPERATURE  AND  DENSITY  -  -  -  - 

515  IF  GLMB{ II  -  1}  <>  0  SOTO  13840:  'checks  to  see  if  temp  changes  with  alt. 

516  R0  -  RH0B( II  -  1)  *  EXP(-(H  -  HB(  II  -  1))  *  GC  /  (1716.483  *  TMBtll  -  1))) 

517  T6  -  TMB( I 1  -  1):  ’  temp  doesnt  change  with  alt 

518  GOTO  13888 

519  13848  T0  «  TMB( II  -  1)  +  GLMB( II  -  1)  *  (H  -  HB(  II  -  1) ) :  'temp  changes  with  alt 

520  R0  *  RH6B ( II  -  1)  *  EXP(-(1!  ♦  GC  /  (1716.483  *  GLMBUl  -  1)))  *  LO€(T0  /  TMB(  II 

520  1))) 

521  13880  ' .  AM8IENT  CONDITIONS  . 

522  RH0  »  R6  *  GC:  '  converts  density  to  lbm/cu.ft 

523  P6  *  RH0  *  T0  *  53.28  /  144:  'pressure  in  psia 

524 

525  C0  *  49.01  *  SQR( T0 ) :  '  ft/sec 


526 

V0  * 

M6  * 

C0 

527 

Q8  * 

RH8 

*  ( V0 

*  2)  / 

(2 

*  GC): 

'  lbf/sqft 

528 

FC  ■ 

(WEIGHT  / 

GC)  » 

(V0 

’  2)  / 

(H6  ♦  REARTH) 

529 

G0  » 

(P0 

*  V6) 

/  (53. 

28  * 

T0 )  : 

lbm/(sec*sqin) 

536 

RETURN 

531 

532  ' 4PAGE 
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533  14060  •*********.*.***.**.*.*****.******♦*.*.*...*.***..*********«****.********. 

534  *  THRUST  CALCULATION  SUBROUTINE 

536  ‘First  update  lift  and  drag: 

537 

538  ITAB  =  3:  XARG  =  MB:  GOSUB  16400 :  CDS  =  VALI:  'Look  up  CdB 

539  ITAB  =  4:  XARG  =  H0:  GOSUB  16400 :  OELCO  =  VALI:  ’Look  up  delCd  friction 

54B  ITAB  =  5:  XARG  =  MB:  GOSUB  164B6:  CLALPHA  =  VALI:  'Look  up  Clalpha 

541  ITAB  *  6:  XARG  =  MB:  GOSUB  16400 :  K  «  VALI:  'LOOK  UP  K 

542  CL  *  ALPHA  *  CLALPHA 

543  LIFT  »  CL  *  QB  *  AREF 

544  CD  =  CDB  +  K  *  CL  *  2  +  DELCO 

545  DRAG  =  CD  *  QB  *  AREF 

546 

547  ' . . 

548  ’  Calculate  commanded  THRUST  value 

549 

55B  THRUSTCOM  =  (WEIGHT  /  (GC  *  COStALPHA  /  ANG)))  *  (ACCCOM  *  GA  +  GA  *  SIN( GAMMA  / 

55B  ANG))  +  DRAG  /  (COS( ALPHA  /  ANG)) 

551  '  . 

552  ’ ' AIRBREATHER  THRUST . 

553  ITAB  =  7:  XARG  =  MB:  YARG  =  ALPHA:  GOSUB  16466:  ABAC  *  VALI:  'Look  up  Afl/AC 

554  ITAB  =  2:  XARG  =  MB:  GOSUB  16400 :  PHIMAX  =  VALI:  'Look  up  PHImax 

555  PHI  =  PHIMAX 

556  ITAB  =  1:  XARG  =  MB:  YARG  =  PHI:  GOSUB  1646B:  ISPA  -  VALI:  'Look  up  airbreathing  I 

556  sp 

557 

558  AB  =  ABAC  *  AC 

559  IF  PHI  =  6  THEN  AB  =  B 

566  WDOTA  =  GB  *  AB  *  144 

561  FAMAX  *  FASTOIC  *  PHI 

562 

563  WDOTFMAX  =  FAMAX  *  WDOTA 

564  THRUSTAMAX  *  ISPA  *  WDOTFMAX 

565  "ROCKET  THRUST . 

566 

567  ITAB  =  8:  XARG  =  HB :  GOSUB  16400 :  ISPR  =  VALI:  'Look  up  rocket  Isp 

568  ITAB  =  9:  XARG  =  MB:  GOSUB  16466:  WDOTPMAX  =  VALI:  'Rocket  propellant  flow 

569 

578  THRUSTRMAX  =  ISPR  *  WDOTPMAX 

571  THRUSTMAX  =  THRUSTAMAX  ♦  THRUSTRMAX 

572  THRUST  =  THRUSTMAX 

573  WOOTF  =  WOOTFMAX 

574  WDOTP  =  WDOTPMAX 

575  WDOT  =  WOOTFMAX  +  WDOTPMAX 

576  ISP  =  THRUST  /  WDOT 

577 

578  IF  THRUST  <  THRUSTCOH  GOTO  15668:  ' . 

579 

58B  ’  . 

581  ' * '  For  THRUSTMAX  >  THRUSTCOM  reduce  rocket  to  zero  if  needed  and 

582  '  reduce  airbreathing  thrust  as  needed  to  obtain  THRUSTCOM 

583 

584  IF  WDOTPMAX  =■  B  THEN  GOTO  15260:  ' . 

585 

586  OELTHRUST  *  THRUSTMAX  -  THRUSTCOM 

587  IF  OELTHRUST  >  THRUSTRMAX  THEN  WDOTP  *  6 

588  IF  WDOTP  *  0  THEN  GOTO  15260:  ' . 

589  '"'Throttle  rocket  to  reduce  thrust 

596  THRUSTR  «  THRUSTCOM  -  THRUSTAMAX 

591  WDOTP  =  THRUSTR  /  ISPR 

592  WOOTF  »  WOOTFMAX 

593  WDOT  *  WOOTF  ♦  WOOTP 

594  THRUST  »  THRUSTAMAX  ♦  THRUSTR 

595  ISP  -  THRUST  /  WDOT 

596  GOTO  15660  :  ' . 

597  15266  IF  THRUST  *  0!  THEN  THRUST  *  I! 

598  IF  ABAC  «  0  OR  PHIMAX  *  B  OR  ISPA  *  0  THEN  GOTO  15660:  ' . 

599  PHI  *  PHIMAX 

660 

601  15346  ' 
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602  ITAB  =  1:  XARG  =  M0 :  YARG  *  PHI:  GOSUB  16400:  ISPA  *  VALI:  'Look  up  Isp 

603 

604  THRUST  *  ISPA  *  WDOTF 

605  WDOTF  *  WDOTF  *  THRUSTCOM  /  THRUST 

606 

607  PHI  =  (WDOTF  /  WOOTA)  /  FASTOIC 

608  ' 

609  IF  ABS( THRUST  -  THRUSTCOM)  <  200  60T0  15560 

610 

611  GOTO  15340:  '  Returns  to  try  new  PHI  value 

612  15560  ' 

613  WDOTF  =  PHI  *  FASTOIC  *  WOOTA 

614  THRUST  =  ISPA  *  WDOTF 

615  WDOT  *  WDOTF 

616  ISP  «  ISPA 

617  15660  ' 

618  IF  WDOTF  *  0  THEN  ISPA  -  0 

619  IF  WDOTP  *  0  THEN  ISPR  =  0 

620 

621  RETURN 

622 

623 

624  ' $PAGE 
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625  15860  ************************************************************************** 

626  '  INTEGRATION  SUBROUTINE 

627  *  2ND  ORDER  R-K  INTEGRATION  (HEUN'S  METHOD) 

628  1.*.«*,*o.***«««..*.*...*.****..*.*..*.*.*.#.**.***»**.**«***.***»*.*. 

629  IF  PLOTS  =  "N"  OR  PLOTS  =  "n"  THE*  60T0  15926 

630  '  . CREATE/AOD  TO  PLOTTER  FILE . 

631  IF  FREQC  <  FREQ  THE*  SOTO  15910 

632  FREQC  =  0 

633  WRITE  #2,  TIMEX,  V0,  H0 ,  RANGE  /  6076,  ISP,  ISPEFF,  ISPMEAN,  WOOTF,  WDOTA,  PHI, 

633  THRUST.  THRUSTCOM,  THRUSTMAX,  Q0,  ALPHA.  GAMMA.  WEIGHT,  LIFT,  DRAG,  WDOTP,  RES1 ,  ISPAVG, 

633  ACC.  ACCCOM,  THRUST  /  DRAG 

634  15910  FREQC  =  FREQC  +  1 

635 

636  15920  1 

637 

638  ' - Set  “old"  values  of  Integration  variables - 

639  TIMEX0  =  TIMEX 

640  WEIGHT0  =  WEIGHT 

641  V00  =>  V0 

642  GAMMA0  =  GAMMA 

643  H00  =  H0 

644  RANGE0  =  RANGE 

645 

646  1 - Start  integration - 

647  FOR  N  =  1  TO  2 

648 

649 

650  ' . . --Weight  increment - - - 

651  OELW(N)  =  (-WOOT)  *  OT  /  1 

652  ' . Velocity  Increment . 

653  IF  TIMEX  <=  DT  THEN  WEIGHT  *  WLAUNCH 

654  DELV(N)  =  (GC  7  WEIGHT)  *  ((THRUST  *  COS(ALPHA  /  ANG)  -  DRAG  -  WEIGHT  *  (GA  /  GC)  * 

654  SIN( GAMMA  /  ANG)))  *  DT  /  1 

655  1 - Flight  path  angle  increment- . - - - 

656  IF  V0  <  DT  THEN  GOTO  15940 

657  DELGAM(N)  «  (GC  /  (WEIGHT  *  V0))  *  ((LIFT  ♦  THRUST  *  SIN(ALPHA  /  ANG)  ♦  FC  *  (GA  7 

657  GC)  *  COS( GAMMA  /  ANG)  -  WEIGHT  *  (GA  /  GC)  *  COS( GAMMA  /  ANG)))  *  DT  /  1 

658  ' - Altitude  Increment - 

659  15940  DELH(N)  *  (V0  *  SIN( GAMMA  /  ANG))  *  DT  /  1 

660  ' . . Range  increment . . -- 

661  DELR(N)  =  (V0  *  COS( GAMMA  7  ANG))  *  DT  /  1 

662 

663  ‘ - Go  back  for  2nd  pass  thru  equations - 

664  TIMEX  -  TIMEX  +  DT  7  2 

665  WEIGHT  =  WEIGHT  +  DELW(N) 

666  V0  =  V0  +  OELV(N) 

667  GAMMA  *  GAMMA  +  OELGAM(N) 

668  H0  -  H0  +  DELH(N) 

669  RANGE  =  RANGE  +  OELR(N) 

670  NEXT  N 


671  '  - After  two  passes  : 

672  ’ - Add  Increments  as  below  of  two - 

673  ’  Integration  steps  to  "old"  values  to 

674  1  obtain  “new"  values  of  Integration 

675  '  variables  at  “new"  time 

676  TIMEX  *  TIMEX0  ♦  DT 

677  WEIGHT  =  WEIGHT0  +  (0ELW(1)  ♦  DELW(2))  /  2 

678  V0  =■  V00  +  (DELV(l)  ♦  DELV(2))  /  2 

679  M0  »  V0  /  C0 

680  GAMMA  =  GAMMA0  +  (DELGAM(l)  +  QELGAM(2) )  7  2 

681  H0  *  H00  ♦  (DELH(l)  +  DELH( 2) )  /  2 

682  RANGE  *  RANGE0  +  (DELR(l)  +  DELR(2) )  7  2 

683  Q0OLD  *  Q0 :  H0OLD  =  H0 :  GAMMAOLD  -  GAMMA:  RH0OLD  *  RH0 

684 

685  ' - End  of  integration  routine-- . 

686  RETURN 

687 

688  'SPAGE 
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689  16000  ********************************************************************** 

690  '  TABLE  LOAOER  SUBROUTINE 

691  1  THIS  SECTION  PROVIDED  BY  KEN  WATSON 

692  I******************************************************************** 

693 

694  IL  «  IL  +  1:  IF  IL  *  1  THEN  ILTI  «  0:  'Initialize  table  location  Index 

695  1 1 ST ( IL)  »  2:  JIST(IL)  *  2 

696  1 0 ( IL)  «  ILTI  +  1:  ITO  =  IO(IL):  'Set  table  origin 

697  INPUT  #1,  T$( IL) ,  AT(ITO):  ‘Read  table  title  &  type 

698  IF  ATI ITO)  *  1  THEN  16090  ELSE  16150 

699  16090  '  1-0  TABLE 

700  IX  *  ITO  +  1:  INPUT  #1,  ATI IX) :  'Input  number  of  points 

701  IY  «  IX  +  AT ( I X ) 

702  FOR  I  -  1  TO  ATI  IX):  INPUT  #1.  AT(I  ♦  IX) :  NEXT  I:  'Input  x 

703  FOR  I  =  1  TO  ATI  IX):  INPUT  #1,  ATI  I  ♦  IY):  NEXT  I:  'Input  y 

704  ILTI  -  ATI  IX)  ♦  IY:  RETURN:  'Set  Index  of  last  point 

705  16150  '  2-0  TABLE 

706  IX  »  ITO  *  2:  INPUT  #1,  ATIITO  *  1),  AT(IX):  'Input  H  and  M 

707  IY  -  IX  +  ATIITO  +1):  IZ  *  IX  ♦  ATI  IX) 

708  FOR  I  *  1  TO  ATIITO  +  1):  INPUT  #1,  ATI  I  ♦  IX):  NEXT  I:  'Input  x 

709  FOR  1  =  1  TO  ATI  IX):  INPUT  #1,  AT(I  +  IY):  NEXT  I:  'Input  y 

710  FOR  J  «  1  TO  ATI  IX):  FOR  I  =  1  TO  ATIITO  ♦  1):  'Input  z 

711  K  -  I  ♦  ATIITO  ♦  1)  *  J  +  IZ:  INPUT  #1,  AT ( K ) :  NEXT  I:  NEXT  J 

712  ILTI  -  ATIITO  ♦  1)  *  (1  +  ATI  IX))  ♦  IZ:  RETURN 

713 

714 

715 

716  ' SPAGE 
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717  16230  ' ************************************************************** 

718  1  TABLE  PRINTER  SUBROUTINE 

719  '  PROVIOEO  BY  KEN  WATSON 

728  1  ************************************************************** 

721  LPRINT  :  LPR1NT  :  LPRINT  "  TJ(ITAB) :  ITO  «  lO(ITAB) 

722  LPRINT  "  *;  AT(ITO);  "-D  TABLE" ,  "STARTS  AT  *;  ITO:  LPRINT 

723  IF  AT ( ITO)  -  1  THEN  1629B  ELSE  1634* 

724  16298  IX  ■  ITO  +1:  IY  -  IX  ♦  AT(IX) 

725  LPRINT  TAB( 18);  "X";  TABI25);  *Y* 

726  FOR  I  *  1  TO  AT(IX) 

727  LPRINT  USING  “#######.#####";  AT ( I  ♦  IX),  AT(I  ♦  IY):  NEXT  I 

728  RETURN 

729  16348  IX  -  ITO  +2:  IY  -  IX  ♦  AT ( ITO  +1):  IZ  -  IX  ♦  AT( ITO  ♦  2) 

738  LPRINT  TA8(  18);  "X";  TABI25);  "Y";  TABI35);  "Z* 

731  FOR  J  -  1  TO  ATI ITO  *  2):  FOR  I  *  l  TO  ATI ITO  ♦  1) 

732  K  -  I  ♦  J  *  ATI ITO  ♦  I)  ♦  IZ 

733  LPRINT  USING  "#######.#####";  ATI  I  ♦  IX),  AT ( J  +  IY),  AT(K) 

734  NEXT  I:  NEXT  J:  RETURN 

735 

736 

737 

738  1 JPAGE 
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739  16466  1  ********************** . ******************************************** 

746  •  TABLE  INTERPOLATION  SUBROUTINE 

741  ’  PROVIDED  BY  KEN  WATSON 

742  •**•**••****••••***•****•••*••****•***•••••**•*•••*•*••***•*•*••„•***** 

743  ITO  ■  IO(ITAB):  IF  AT(ITO)  -  1  THEN  16446  ELSE  16596 

744  16446  *  1-D  table 

745  IX  -  ITO  ♦  1 :  N  =  AT(IX):  IY  *  IX  +  N 

746  IF  XARG  <*  AT( IX  +  1)  THEN  16476  ELSE  16496 

747  16476  I  »  2 

748  16486  12  *  I  ♦  IX:  II  *  12  -  1:  J2  «  I  ♦  IY:  J1  *  J2  -  1:  60T0  16576 

749  16496  IF  XARG  >*  AT{ IX  +  N)  THEN  16566  ELSE  16516 
756  16566  I  •  N:  GOTO  16486 

751  16516  I  «  IIST(ITAB) 

752  16526  12  *  I  ♦  IX:  II  *  12  -  1 :  J2  *  I  ♦  IY:  J1  *  J2  -  1 

753  IF  XARG  <*  AT(  12)  THEN  16546  ELSE  16566 

754  16546  IF  XARG  >  AT(I1)  THEN  16576  ELSE  16558 

755  16556  1*1-1:  GOTO  16528 

756  16566  1*1+1:  GOTO  16528 

757  16576  VALI  =  AT(J1)  +  ( AT { J2 )  -  AT(J1))  *  ((XARG  -  AT (II))  /  (AT( 12)  -  AT(I1))) 

758  IIST(ITAB)  =  I:  RETURN 

759 
766 

761 

762  16596  1  2-0  table 

763  IX  -  ITO  +  2:  N  *  AT( ITO  +  1) :  M  =  AT( ITO  +2):  IY  *  IX  ♦  N:  IZ  -  IX  +  M 

764  IF  XARG  <=  AT ( 1  ♦  IX)  THEN  16636  ELSE  16626:  ’Locate  II  &  12 

765  16626  IF  XARG  >*  AT(N  +  IX)  THEN  16646  ELSE  16656 

766  16636  1*2:  GOTO  16716 

767  16646  I  *  N:  GOTO  16716 

768  16656  I  *  IIST(ITAB) 

769  16666  12  *  I  ♦  IX :  II  *  12  -  1 

776  IF  XARG  <=  AT { 12 )  THEN  16686  ELSE  16766 

771  16686  IF  XARG  >  AT ( I 1 }  THEN  16716  ELSE  16696 

772  16696  1*1-1:  GOTO  16666 

773  16766  I  *  I  +  1:  GOTO  16666 

774  16716  IIST(ITAB)  *  I 

775  12  *  I  ♦  IX:  II  *  12  -  1 

776  IF  YARG  <*  AT ( 1  +  IY)  THEN  16756  ELSE  16746:  ’Locate  J1  &  J2 

777  16740  IF  YARG  >*  AT(M  +  IY)  THEN  16766  ELSE  16776 

778  16756  J  *  2:  GOTO  16836 

779  16766  J  *  M:  GOTO  16838 

786  16776  J  -  JIST(ITAB) 

781  16786  J2  *  J  +  IY:  J1  *  J2  -  1 

782  IF  YARG  <=  AT(J2)  THEN  16866  ELSE  16826 

783  16866  IF  YARG  >  AT(J1)  THEN  16836  ELSE  16816 

784  16816  J  *  J  -  1:  GOTO  16786 

785  16826  J  *  J  +  1:  GOTO  16786 

786  16836  JIST(ITAB)  =  J 

787  J2  *  J  ♦  IY:  J1  *  J2  -  1 

788  Kll  -  I  -  1  ♦  N  *  (J  -  1)  ♦  IZ:  K21  *  I  +  N  *  (J  -  1)  ♦  IZ 

789  K12  *  I  -  1  +  N  *  J  ♦  IZ:  K22  *  I  ♦  N  *  J  ♦  IZ 

796  DXI  *  (XARG  -  AT ( 1 1 ) )  /  (AT(  12)  -  AT (II)) :  ’Have  II  ,12 ,J1  ,J2 ,K11  ,K21  ,K12 ,K22 

791  OYJ  *  (YARG  -  AT ( J 1 ) )  /  ( AT ( J2 )  -  AT(J1)) 

792  VALI  -  AT (Kll)  +  (AT(K21)  -  AT(Kll))  *  DXI 

793  VAL2  *  AT(K12)  +  (AT(K22)  -  AT(K12))  *  DXI 

794  VALI  -  VALI  ♦  ( VAL2  -  VALI)  *  OYJ 

795  RETURN 

796 

797 

798  ' IPAGE 
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799  18001)  ********************************************************************* 

800  '  SUBROUTINE  TO  LOAD  SECOND  STAGE  DATA  FILE 

30 1  • ft******************************************************************* 

802  OPEN  "I",  #1.  IFN2$:  IL  »  0:  IFILE  *  1 

803  ' 

804  LINE  INPUT  #1.  TITLU 

805  LINE  INPUT  #1,  T I TL2 $ 

806 

807  GOSUB  16000:  'Input  airbreather  isp  vs.  Mach  and  PHI  as  table  fl 

808  GOSUB  16000:  'Input  airbreather  PHImax  vs.  Mach  as  table  #2 

809  GOSUB  16000:  'Input  Cdo  vs. Mach  as  table  13 

810  GOSUB  16000:  'Input  delta  Cd  vs.  Alt.  as  table  #4 

811  GOSUB  16000:  'Input  Clalpha  vs.  Mach  as  table  #5 

812  GOSUB  16000:  'INPUT  K  VS.  MACH  AS  TABLE  #6 

813  GOSUB  16000:  'Input  airbreather  Ao/Ac  vs.  Mach  and  alpha  as  table  #7 

814  GOSUB  16000:  'Input  rocket  Isp  vs.  Altitude  as  table  #8 

815  GOSUB  16000:  ‘Input  rocket  UOOTPMAX  vs. Mach  as  table  #9 

816 

817  INPUT  #1,  AC,  AREF 

818  INPUT  #1.  HEIGHT,  WFUEL.  UFINAL ,  VFINAL,  STAGE 

819  INPUT  #1,  DT,  OELPRINT 

820  INPUT  *1.  Q0COMD,  Q0MAX  ,  Q0MIN ,  ALPHAMAX  ,  ALPHAMIN 

821 

822  INPUT  #1,  LOADFAC,  ACCCOMD,  GAMMAMAX 

823  INPUT  #1.  FASTOIC 

824  INPUT  I I,  V0TEMP,  Q0FINAL 

825  INPUT  #1.  V0CRUISE,  SHITCH4 ,  GAMMA4MAX 

826  INPUT  #1,  FLAG5,  V05,  ALPHA5MAX ,  GAMMA5 ,  ACCC0MD5 

827  INPUT  #1,  Al,  A2 ,  B1 ,  Cl,  C2,  Dl,  DELQ,  V02,  H02 

828  ' 

829  PRINT 

830  PRINT  "INPUT  FILE  FOR  2ND  STAGE  HAS  BEEN  READ." 

831  PRINT  "LAST  ITEM  IN  2ND  STAGE  INPUT  FILE  IS  H02, VALUE  READ  FROM  FILE  IS";  H02 

832  PRINT  "IF  THIS  VALUE  IS  INCORRECT  THEN  INPUT  FILE  IS  IN  ERROR." 

833  PRINT 

834  INPUT  "Ready  to  run  2nd  stage: 

834  Do  you  want  a  listing  of  2nd  stage  tabular  Input  data(Y/N)";  HTINJ 

835  IF  HTIN$  *  "Y"  OR  WTIN$  =  "y"  THEN  18230  ELSE  18240 

836  18230  FOR  I TAB  *  1  TO  9:  GOSUB  16230:  NEXT  I  TAB 

837  18240  TIMEX2  *  TIMEX 

838  IF  EOF(l)  THEN  GOTO  18250 

839  18250  CLOSE  fl 

840  RETURN:  'Returns  to  main  routine 

841 

842 

843  ' SPAGE 
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844  21000  ' PRIN . BAS 

845 


POST-RUN  PRINTER  OUTPUT  SUBROUTINE 


'  8 MAY87/ 0030 /PROGRAM  "ET0TA8. BAS" .MERGED  FOR  USE  WITH  PROGRAM 
N  *  55:  'Number  of  lines  printed  per  page 

WIDTH  LPRMT  136 


LPRIWT  TAB(  40 )  : 
LPRIWT 

LPRIRT  T ABC  10) 
LPRIHT  T ABC  10) 
LPRIlfT 

LPRIWT  T ABC  10 ) 

LPRIWT 

LPRIWT 

LPRIWT  TABC10) 
LPRIWT  USIWG  " 
LPRIWT  USING  " 


DATE$;  T ABC  60);  TINES 

" IST/SINGLE  STAGE  INPUT  DATA  FILE: 
"2ND  STAGE  INPUT  DATA  FILE  : 

"PLOTTER/PRINTER  FILE  NAME  : 


LPRIWT  TABC10);  " . PROGRAM  ET02  RUN  SUMMARY . " 

LPRIWT  USING  "  LAUNCH  WEIGHT  =  tmttlt.tt  LBS.";  WLAUNCH 

LPRIWT  USING  ”  FINAL  WEIGHT  =  ittlttti.lt  LBS.";  WEIGHT 

LPRIWT  USING  "FUEL/PROPELLANT  REMAINING  =#f#0ffff.#f  LBS.";  WEIGHT  -  WFINAL 

LPRIWT 

IF  V0  >=  VORBIT  THEN  LPRIWT  TABC10) ;  "VEHICLE  REACHED  ORBITAL  VELOCITY 


IF  WEIGHT  <=  WFINAL  THEN  LPRIWT  T ABC  10) ; 
LPRIWT 

LPRIWT  USING  "FLIGHT  TIME  =  tltltlt.l 
LPRIWT  USING  "FINAL  VELOCITY  =  llllttl  .1 
LPRIWT  USING  "FINAL  ALTITUDE  =  # 111111.1 
LPRIWT  USING  "FINAL  Q0  =  »###### . # 

LPRIWT 

LPRIWT  USING 
LPRIWT  USING 


’RUN  TERMINATED  BY  FUEL  EXHAUSTION* 


1111111.11  SECS.";  TIMEX 

#*#«#*#.  It  FT/SEC";  V0 
Itltltt.lt  FT"  ;  H0 
####### . It  PSF"  ;  00 


"AVERAGE  ISP 
"MEAN  ISP 


nnnt.tt  sec 
nnitt.tt  sec 


ISPAVG 

ISPMEAN 


875  LPRIWT 

876  LPRIWT 

877  LPRIWT 

878  LPRIWT  "FINAL  CONDITIONS  ARE  SHOWN  BELOW:" 

879  LPRIWT 

880  WIDTH  LPRIWT  136 

881  LPRINT  CHRS  C  27 )  ;  "B";  CHRS  (  3 ) : '  Se  1  ec  ts  condensed  printCSRIS  printer) 

882  LPRIWT  CHRSC27) ;  CHRS (33) ;  CHRSC4) ;  :  'Selects  condensed  printCEpson  printer) 

883  LPRIWT  TABC5) ;  "TIME";  TABC15);  "V0";  TABC25) ;  "H0‘;  TABC35) ;  "R";  TABC45) ;  "ISP"; 
883  TAB! 55) ;  "ISPEFF";  TAB! 65) ;  "ISPMEAN";  TA8C75);  "WDOTF"  ;  TABC85) ;  "WDOTA" ;  TABC97) ;  "PHI" 


TABC85) 

"T/D" 


*•  TMBlKTMA  *  " 

"GAMMA";  TABC45);  "WEIGHT";  TABC55 
" RES1 " ;  TABC94);  "ISPAVG";  T ABC  102 


883  ;  TABC105);  "THRUST";  TABCl 12);  "THRUSTCOM";  TABC122) ;  "THRUSTMAX" 

884  LPRIWT  TABC15) ;  "Q0"  ;  TABC25) ;  "ALPHA";  TABC35) ;  "GAMMA";  TABC45) ;  "WEIGHT";  T 
884  );  "LIFT";  TABC65);  “ORAG"  ;  TABC75);  "WDOTP"  ;  TABC85) ;  "RES1";  TABC94);  "ISPAVG";  TA 

884  );  "ACCEL-G'S";  TABC115);  "ACCCOM"  ;  T ABC  128) ;  "T/D" 

885  LPRIWT 

886  LPRIWT  USING  “####### .  ##"  ;  TIMEX;  V0 ;  H0 ;  RANGE  /  6076;  ISP;  ISPEFF;  ISPMEAN; 

886  WD1TF;  WDOTA;  PHI;  THRUST;  THRUSTCOM;  THRUSTMAX 

887  LPRIWT  USING  "#######.##";  TIMEX;  Q0 ;  ALPHA;  GAMMA;  WEIGHT;  LIFT;  DRAG;  WDOTP; 

887  REf 1 ;  ISPAVG;  ACC;  ACCCOM;  THRUST  /  DRAG 

888  WIDTH  LPRINT  80 

889  LPRIWT  CHRSt 12) :  'Formfeed 

890  '  LPRINT  CHRS (27)  "B"  C HRS (  1)  : 'Resets  to  standard  print(SR15  printer) 

891  LPRIWT  CHRS (27) ;  CHRS (33);  CHRSC0);  :  'Resets  to  standard  printCEpson  printer) 

892 

893  CHOIS  =  "COPY"  ♦  CHRS! 32)  ♦  I FNS  ♦  CHRS 132)  +  " PRN" 

894  PRINT  CHOIS 

895  SHELL  CM01S 

896  LPRIWT  CHRSC 12) :  'Form  ‘eed 

897  IF  NSTAGES  *  1  THEN  GOTO  20820 

898  CMD2S  *  "COPY"  ♦  CHRSC 32 )  ♦  IFN2S  ♦  CHRSC 32)  ♦  "PRN" 

899  PRINT  CMD2S 

900  SHELL  CMD2S 

901  LPRIWT  CHRSC  12) 

902  20820  WIDTH  LPRIWT  136 

903  '"LPRINT  CHRS  ( 27 )  "B"  CHRS  ( 3 )  :  'Selects  condensed  print(SR15  printer) 

904  LPRINT  CHRS (27);  CHRS (33);  CHRS ( 4 ) ;  :  'Selects  condensed  printCEpson  printer) 


DIM  V( 25) 
COUNT  *  0 
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919 

918 

911 

912 

913 

914 

915 

916 

917 

918 

918 

919 

920 

921 

922 

923 

924 

925 

926 

927 

928 

929 

930 

931 

932 

933 

934 

935 

935 

936 

937 

938 

939 
948 

941 

942 

943 

944 

945 

946 

947 

948 

949 

950 

951 

951 

952 

953 

954 

955 

956 

957 

958 

959 
968 

961 

962 

963 

964 

965 

966 

967 

967 

968 

969 
978 

971 

972 

973 

974 

975 

976 


LPRINT  TAX  18)  ;  PFN$;  TAB(40);  DATE*;  TAB(60);  TIME! 

GOSUB  22460 

I 

OPEN  "I",  2,  PFN$ :  INPUT  #2,  VARNUM 
21140  FOR  I  *  1  TO  25 

INPUT  #2,  V(I) 

IF  EOF C  2 )  THEN  GOTO  21340 

NEXT  I 

COUNT  *  COUNT  +  1 

LPRINT  USING  "##«###*#*  .##" ;  V(l);  V { 3 ) ;  V ( 2 ) ;  V ( 1 4 ) ;  V(4);  V( 17 )  ;  V(8) ;  V(9) ;  V(18 
);  V  ( 20) ;  V  ( 2 1 ) 

IF  COUNT  >=  N  THEN  LPRINT  CHR$(12) 

IF  COUNT  >=  N  THEN  COUNT  *  0 
IF  COUNT  =  0  THEN  GOSUB  22460 
GOTO  21140 
21340  CLOSE  2 

t  I  •  I  I  t  «  I  ' 

GOTO  21700:  ’SKIPS  OVER  2ND  FIELD  FOR  PROGRAM  “ETO" 

LPRINT  CHR$( 12) :  ’FORM  FEEO 
COUNT  =  0 

GOSUB  22500 

OPEN  "I",  2,  PFN$:  INPUT  #2.  VARNUM 
21480  FOR  I  *  1  TO  25 

INPUT  #2,  V ( I ) 

IF  EOF ( 2 )  THEN  GOTO  21680 

NEXT  I 

COUNT  *  COUNT  +  1 

LPRINT  USING  ’## ffff* ##.##“  ;  V(l);  V(3);  V(2);  V(8) ;  V(9);  V( 10) ;  V(29);  V(28);  V(2 

0) 

IF  COUNT  >  =  N  THEN  LPRINT  CHR$ (12) 

IF  COUNT  >=  N  THEN  COUNT  =  0 
IF  COUNT  *  0  THEN  GOSUB  22500 
GOTO  21480 
21680  CLOSE  2 
21700  . 

LPRINT  CHR$(12):  ’ FORM  FEEO 
COUNT  *  0 

GOSUB  22540 

OPEN  "I*.  2,  PFN$ :  INPUT  #2.  VARNUM 
21800  FOR  I  =  1  TO  25 

INPUT  #2,  V ( I ) 

IF  EOF { 2 )  THEN  GOTO  22000 

NEXT  I 

COUNT  =  COUNT  +  1 

LPRINT  USING  “##ff f Iff f  .##" ;  V ( 1 ) ;  V ( 3 ) ;  V(2);  V(5);  V(22);  V(6);  V(7);  V(ll);  VC  12 
);  V( 13) 

IF  COUNT  >=  N  THEN  LPRINT  CHR$(12) 

IF  COUNT  >=  N  THEN  COUNT  =  0 
IF  COUNT  *  0  THEN  GOSUB  22540 
GOTO  21800 

22000  CLOSE  2 

1*1(11 

LPRINT  CHRX12)  :  'FORM  FEEO 
COUNT  »  0 
GOSUB  22580 

OPEN  "I".  2,  PFN$ :  INPUT  12,  VARNUM 
22120  FOR  I  -  1  TO  25 

INPUT  #2,  V(I) 

IF  EOF ( 2 )  THEN  GOTO  22320 

NEXT  I 

COUNT  -  COUNT  ♦  1 

LPRINT  USING  "ffMHIH  .#»“  ;  V(l);  V  ( 3 )  ;  V(2);  V(  14)  ;  V  ( 1 5 )  ;  V(I6);  V(  18)  ;  V(19); 

V ( 23 ) ;  V(24)  ;  V(25) 

IF  COUNT  >=  N  THEN  LPRINT  CHR»(12) 

IF  COUNT  >-  N  THEN  COUNT  *  0 
IF  COUNT  *  0  THEN  GOSUB  22580 
GOTO  22120 
22320  CLOSE  2 

«  (  #  I  I  I 
I  I  I  l  I  I  I  I 

LPRINT  CHR$(27) ;  "B";  CHRt(l):  'Resets  to  standard  print 

LPRINT  CHRX27);  CHR$(33);  CHR$(0);  :  '  Resets  to  standard  print(Epson  printer) 
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977  WIDTH  LPRIHT  80 

978  LPRIHT  CHR|( 12) :  'Formfeed 

979  60T0  22999 

980  22460  LPRIHT  TAB[6) ;  "TIME";  TAB(18) ;  “H0";  TAB(30);  “V0";  TAB(42);  "Q0";  TAB(54) ;  "R"; 

980  TAB(66);  "WEIGHT";  TA8(78);  "WDOTF  ";  TAB(90);  "  WDOTA  ";  TAB(102);  "PHI  ";  TABC114) 

980  ;  "  WDOTP  ";  TAB( 125) ;  "  RES1  " 

981  RETURH 

982  22500  LPRIHT  TAB(6) ;  "TIME";  TAB(18);  "H0";  TAB(30);  "V0";  TAB(42);  "WOOTF";  TAB(54);  "WD 

982  OTA";  TAB(66);  "PHI";  TAB(78);  "WDOTFROC"  ;  TABI90);  "WDOTLOX";  TABI102);  "WDOTP* 

983  RETURH 

984  22540  LPRIHT  TAB(6) ;  "TIME”;  TABI18);  "H0“;  TAB(30);  "V0";  TAB(42);  "ISP";  TABC54) ;  "ISPA 
984  VG";  TAB(66) ;  "ISPEFF";  TAB( 78) ;  "ISPMEAN";  TAB(90) ;  "THRUST";  TABI102);  "THRUSTCOM"; 

984  TABI114) ;  "THRUSTMAX" 

985  RETURH 

986  22580  LPRIHT  TABI6) ;  “TIME";  TAB(18);  “H0";  TA8( 30 ) ;  "V0“;  TAB(42);  "Q0";  TAB(54);  "ALPHA 
986  ";  TABI66);  "GAMMA";  TAB(78);  "LIFT";  TA8(90) ;  "DRAG";  TAB(102);  "ACCEL-G'S";  TAB(114);  " 

986  ACCCOM"  ;  TABU29)  ;  "T/D" 

987  RETURH 

988  22999  RETURH 

989 

990 
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5.0  CROSS  PEFEREtlCEL  LISTIMO  OF  PROGRAM  FTP  VARIABLES 


Each  variable  used  in  the  program  is  listed  alphabetically  along  with  the  sequential 
program  line  numbers  fcr  each  occurrence.  The  line  numbers  arc  consistent  with  4>e 
program  listing  in  Section  4.0. 
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A* 

192 

192 

226 

391 

AD 

558 

559 

560 

ABAC 

180 

553 

558 

598 

AI 

91 

268 

827 

A2 

91 

268 

827 

AC 

81 

558 

817 

ACC 

378 

399 

422 

633 

887 

ACCCOM 

221 

240 

241 

294 

303 

887 

ACCCOMO 

86 

221 

240 

294 

302 

ACCC0M05 

90 

357 

826 

ACCPERQ 

302 

303 

ALPHA 

83 

180 

184 

220 

223 

273 

278 

278 

283 

283 

336 

339 

342 

343 

344 

550 

553 

633 

654 

657 

ALPHA2MAX 

85 

239 

239 

ALPHA5MAX 

90 

358 

358 

826 

ALPHACRUISE 

339 

342 

ALPHAMAX 

84 

285 

285 

326 

345 

ALPHAMIN 

84 

286 

286 

397 

820 

ANG 

143 

550 

550 

550 

654 

AREF 

81 

185 

187 

238 

543 

AT 

41 

697 

698 

700 

701 

707 

707 

708 

708 

709 

722 

723 

724 

726 

727 

733 

733 

743 

745 

746 

757 

757 

763 

763 

764 

790 

790 

790 

791 

791 

B$ 

244 

391 

81 

91 

273 

827 

BIGG 

141 

203 

C$ 

315 

391 

C0 

525 

526 

679 

Cl 

91 

278 

827 

C2 

91 

278 

827 

CO 

186 

187 

332 

333 

544 

CDB 

176 

186 

327 

332 

538 

CL 

184 

185 

186 

331 

332 

CLALPHA 

178 

184 

238 

329 

331 

CMOl  $ 

893 

894 

895 

CM02* 

898 

899 

900 

COUNT 

149 

388 

388 

389 

404 

927 

934 

934 

936 

937 

953 

954 

959 

966 

966 

CPCV 

504 

OJ 

348 

349 

391 

01 

91 

283 

827 

OELCO 

177 

186 

328 

332 

539 

OELGAM 

43 

657 

667 

680 

680 

DELH 

43 

659 

668 

681 

681 

OELPRINT 

83 

389 

819 

OELQ 

91 

271 

272 

281 

282 

OELR 

43 

661 

669 

682 

682 

OELTHRUST 

586 

587 

OELV 

43 

654 

666 

678 

678 

DELW 

43 

651 

665 

677 

677 

DRAG 

187 

398 

399 

400 

422 

887 

OT 

83 

205 

208 

209 

209 

280 

281 

378 

651 

653 

819 

DV 

207 

208 

OX  I 

790 

792 

793 

OYJ 

791 

794 

E 

139 

140 

141 

503 

E$ 

363 

391 

FAMAX 

561 

563 

FASTOIC 

87 

561 

607 

613 

823 

FC 

142 

397 

528 

657 

FINISH 

458 

459 

FLAGS 

90 

261 

314 

826 

305 

303 

322 

822 

357 

398 

422 

550 

633 

238 

239 

239 

242 

268 

268 

273 

285 

285 

286 

286 

326 

331 

333 

345 

887 

358 

359 

393 

422 

542 

550 

358 

396 

820 

654 

657 

657 

657 

659 

661 

545 

817 

702 

702 

703 

703 

704 

706 

706 

709 

710 

710 

711 

711 

712 

712 

727 

729 

729 

731 

731 

732 

733 

749 

753 

754 

757 

757 

757 

757 

765 

770 

771 

776 

777 

782 

783 

791 

792 

792 

792 

793 

793 

793 

545 

544 

333 

542 

543 

544 

540 

542 

908 

917 

917 

919 

920 

920 

921 

937 

938 

943 

950 

950 

952 

953 

968 

969 

969 

970 

544 

827 

422 

545 

550 

633 

633 

654 

887 

226 

235 

243 

267 

270 

271 

277 

654 

656 

657 

659 

661 

664 

676 
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FLAGCRUISE 

324 

341 

FREQ 

119 

631 

FREQC 

631 

632 

G0 

529 

566 

GA 

42 

263 

GAMMA 

83 

220 

657 

657 

GAMMA0 

642 

686 

GAMMA4MAX 

89 

343 

GAMMA 5 

96 

361 

GAMMAMAX 

86 

822 

GAMMAOLO 

683 

GC 

138 

142 

654 

654 

GLMB 

44 

155 

H 

566 

509 

H0 

83 

142 

394 

421 

683 

869 

H00 

643 

681 

H02 

91 

95 

H0OLO 

312 

683 

HB 

44 

155 

HOEL 

509 

516 

I 

151 

196 

709 

709 

733 

734 

756 

758 

788 

788 

949 

962 

11 

568 

569 

526 

520 

1 1ST 

41 

695 

12 

748 

748 

IFILE 

52 

802 

IFNI 

47 

49 

IFN2$ 

50 

862 

IL 

52 

694 

I L  T I 

694 

696 

10 

41 

696 

ISP 

376 

396 

ISPA 

466 

556 

ISPAVG 

377 

399 

ISPEFF 

372 

373 

ISPMEAN 

374 

398 

ISPR 

461 

567 

ITAB 

76 

70 

539 

548 

751 

758 

ITCUM 

373 

373 

I  TO 

696 

697 

722 

722 

745 

763 

ITSUM 

376 

376 

IX 

766 

766 

707 

708 

729 

733 

764 

765 

IY 

701 

703 

763 

776 

!Z 

707 

711 

J 

710 

711 

784 

785 

J1 

153 

155 

787 

791 

J1ST 

41 

695 

J2 

748 

748 

K 

179 

186 

Kll 

788 

792 

K12 

789 

793 

K21 

788 

792 

K22 

789 

793 

634 

634 

378 

556 

550 

654 

657 

223 

343 

345 

361 

394 

659 

825 

826 

661 

667 

667 

680 

263 

37  2 

382 

502 

516 

657 

657 

657 

515 

519 

526 

516 

519 

177 

203 

220 

223 

237 

506 

566 

528 

539 

567 

886 

237 

241 

242 

827 

831 

509 

516 

519 

427 

762 

702 

762 

763 

709 

710 

711 

711 

726 

747 

748 

748 

750 

751 

766 

767 

768 

769 

772 

789 

789 

913 

914 

916 

963 

965 

511 

515 

516 

516 

516 

748 

752 

754 

757 

757 

751 

758 

768 

774 

752 

752 

753 

757 

769 

52 

853 

893 

854 

898 

694 

694 

695 

695 

696 

704 

712 

696 

721 

743 

421 

576 

595 

616 

633 

564 

598 

602 

604 

614 

422 

633 

872 

887 

397 

421 

633 

886 

421 

633 

873 

886 

576 

591 

619 

176 

177 

178 

179 

186 

541 

553 

554 

556 

567 

768 

774 

780 

786 

836 

374 

698 

766 

766 

766 

767 

723 

724 

729 

729 

729 

763 

763 

377 

761 

701 

762 

762 

763 

709 

710 

712 

724 

724 

745 

745 

745 

746 

7  48 

769 

775 

704 

707 

709 

724 

727 

777 

781 

787 

712 

729 

732 

763 

788 

711 

731 

732 

733 

734 

785 

786 

787 

788 

788 

155 

155 

155 

156 

748 

791 

786 

786 

752 

752 

757 

781 

781 

330 

332 

541 

544 

711 

792 

793 


657 

422  556  633  642  654 

683  887 


526  522  527  528  556 


241  242  312  328  382 
633  643  668  668  681 


703 

703 

708 

768 

768 

727 

727 

727 

731 

732 

752 

752 

755 

755 

756 

772 

773 

773 

774 

775 

936 

931 

933 

946 

947 

517 

519 

519 

519 

526 

769 

771 

775 

790 

796 

769 

776 

775 

775 

796 

696 

697  802 

886 

616 

618 

327 

328 

329 

336 

538 

568 

836 

662 

721 

721 

743 

768 

716 

711 

712 

721 

731 

731 

732 

743 

743 

764 

766 

766 

767 

707 

724 

726 

727 

729 

729 

749 

752 

763 

763 

763 

729 

733 

745 

748 

752 

788 

789 

789 

778 

779 

780 

781 

784 

789 

789 

752 

757 

757 

781 

783 

782 

787 

787  791 

711 

732 

733 

70 
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LIFT 

185 

220 

223 

344 

396 

LOADFAC 

86 

238 

822 

LOVRD 

333 

333 

334 

335 

LOVROOLD 

334 

335 

M 

763 

763 

777 

779 

H6 

83 

176 

178 

179 

180 

541 

553 

554 

556 

568 

HEARTH 

139 

203 

N 

647 

651 

654 

657 

659 

745 

745 

749 

750 

763 

848 

919 

920 

936 

937 

N1 

501 

508 

NSTAGES 

45 

46 

443 

444 

897 

P8 

523 

529 

PFN* 

100 

101 

856 

909 

912 

PHASE 

216 

230 

252 

321 

354 

PHI 

395 

421 

555 

556 

559 

PHIMAX 

554 

555 

598 

599 

PLOTS 

98 

99 

99 

446 

446 

629 

629 

PRCOM 

73 

73 

74 

407 

434 

PRCOMI 

72 

73 

73 

PRCOUNT 

150 

409 

410 

412 

418 

Q* 

194 

194 

197 

198 

Q8 

185 

187 

237 

238 

264 

277 

280 

281 

281 

282 

527 

543 

545 

633 

683 

Q0COM 

259 

262 

264 

270 

271 

oecoMD 

84 

237 

259 

262 

262 

Q8FINAL 

88 

262 

824 

Q0MAX 

84 

302 

305 

820 

Q8MIN 

84 

820 

QflOLD 

310 

683 

QOLD 

265 

267 

271 

275 

277 

R0 

516 

520 

522 

RANGE 

395 

421 

633 

644 

669 

RANGE0 

644 

682 

RATEDELQ 

270 

273 

280 

283 

RATEOQ 

267 

268 

271 

277 

278 

RE 

503 

506 

506 

REARTH 

140 

142 

203 

382 

382 

RES1 

393 

422 

633 

887 

RH0 

311 

522 

523 

527 

683 

RH0B 

44 

155 

516 

520 

RH0OLD 

311 

683 

STAGE 

82 

385 

444 

818 

START 

120 

459 

STARTTIME* 

120 

458 

SWITCH4 

89 

323 

349 

825 

T* 

41 

697 

721 

T0 

517 

519 

520 

523 

525 

TG 

42 

THRUST 

400 

400 

402 

421 

422 

604 

605 

609 

614 

633 

THRUSTAMAX 

564 

571 

590 

594 

THRUSTCOM 

401 

421 

550 

578 

586 

THRUSTMAX 

402 

421 

571 

572 

586 

THRUSTR 

590 

591 

594 

THRUSTRMAX 

570 

571 

587 

TIMEX 

83 

148 

218 

235 

392 

440 

445 

633 

639 

653 

TIMEX0 

639 

676 

TIMEX2 

445 

837 

TITL1* 

55 

412 

804 

TITL2* 

56 

412 

805 

TMB 

44 

155 

516 

517 

519 

TQ0 

235 

243 

243 

TSET5 

261 

359 

360 

361 

361 

TTAKEOFF 

218 

226 

226 

235 

V 

907 

914 

918 

918 

918 

918 

931 

935 

935 

935 

951 

951 

951 

951 

951 

398 

398 

422 

543 

633 

657 

887 

327 

329 

330 

393 

526 

538 

540 

602 

679 

661 

665 

666 

667 

668 

669 

670 

763 

765 

767 

788 

788 

789 

789 

952 

953 

968 

969 

929 

945 

961 

391 

391 

391 

391 

391 

561 

599 

602 

607 

613 

633 

886 

451 

451 

452 

453 

453 

455 

455 

437 

440 

446 

418 

425 

425 

265 

267 

270 

271 

271 

272 

275 

284 

291 

303 

305 

310 

395 

422 

870 

887 

272 

280 

281 

282 

291 

302 

303 

820 


281  284 

669  682  886 

281 

528 


529 


572 

576 

578 

594 

595 

597 

597 

633 

654 

657 

886 

887 

590 

605 

609 

633 

886 

633 

886 

421 

422 

433 

434 

436 

437 

439 

664 

664 

676 

837 

867 

886 

887 

528 


918 

918 

918 

918 

918 

918 

918 

935 

935 

935 

935 

935 

935 

947 

951 

951 

951 

951 

951 

963 

967 
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967 

967 

967 

967 

967 

V0 

142 

206 

207 

218 

237 

383 

385 

393 

421 

433 

528 

529 

633 

641 

656 

864 

868 

886 

VBB 

372 

378 

641 

678 

VB2 

91 

237 

241 

827 

V05 

90 

314 

826 

VBCRUISE 

89 

313 

313 

825 

VBTEMP 

88 

260 

262 

262 

824 

VAL1 

792 

794 

794 

VAL2 

793 

794 

VAL I 

176 

177 

178 

179 

180 

541 

553 

554 

556 

567 

VARNUM 

102 

103 

912 

929 

945 

VFINAL 

82 

385 

818 

VORBIT 

382 

383 

864 

VTAKEOFF 

85 

218 

MOOT 

575 

576 

593 

595 

615 

WDOTA 

394 

421 

560 

563 

607 

WDOTF 

393 

421 

573 

592 

593 

618 

633 

886 

WDOTFMAX 

563 

564 

573 

57  5 

592 

WDOTP 

394 

422 

574 

587 

588 

WOOTPMAX 

568 

570 

574 

575 

584 

WEIGHT 

142 

148 

205 

220 

223 

377 

384 

396 

396 

402 

528 

550 

633 

640 

653 

818 

861 

862 

865 

887 

WEIGHTS 

372 

373 

376 

640 

677 

WFINAL 

82 

205 

205 

384 

818 

WFUEL 

82 

818 

WLAUNCH 

82 

148 

374 

377 

653 

WTIN$ 

68 

69 

69 

834 

835 

X* 

391 

391 

391 

391 

391 

XARG 

176 

177 

178 

179 

180 

541 

553 

554 

556 

567 

764 

765 

770 

771 

790 

YARG 

180 

553 

556 

602 

776 

967 

967 

967 

967 

967 

241 

260 

262 

313 

314 

372 

378 

434 

436 

437 

439 

440 

526 

527 

657 

659 

661 

666 

666 

678 

679 

327 

328 

329 

330  538  539  540 

568 

961 

602 

757 

794 

651 

613 

633 

886 

604 

605 

605  607  613  614  615 

591 

593 

619 

633 

887 

238 

344 

372 

372 

373 

374 

376 

422 

433 

434 

436 

437 

439 

440 

654 

654 

657 

657 

665 

665 

677 

862 

865 

860 

835 

392 

327 

328 

329 

330 

538 

539 

540 

568 

602 

746 

749 

753 

754 

757 

777 

782 

783 

791 

fi.o  DEFINITION  nF  PROGRAM  LTO  VARIABL  FS 

The  program  variables  mi  t^is  section  ar<-  1  is*H  al  phabotical  1  y  hv  the 
■following  groupings: 


Simple  Peal  Variables 
Simple  String  Variable? 
Array  Real  Variables 
Array  String  Variable? 
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VARIABLE  LISTING  FOR  PROGRAM  “ETO" 


Simple  Real  Variables 

AO  Air  capture  streamtube.sq. ft , 

ABAC  Air  capture/cowl  area  ratio. Values  are 

input  as  table  6  of  input  data  file. 

A1  Adjustment  “knob"  for  controlling  Phase  3 

constant  q  path . Presently  set  to  8.1 

A2  Adjustment  “knob"  for  controlling  Phase  3 

constant  q  path. Presently  set  to  0.1 

AC  Air  inlet  cowl  area, sq.ft. 

ACC  Axial  acceleration, g's. Output  is  labeled 

ACCEL-G'S. 

ACCCOM  Commanded  axial  acceleration, g's 

ACCCOMD  Maximum  axial  acceleration^'  s 

ACCC0MD5  Commanded  acceleration  in  Phase  5 , g ' s 

ACCPERQ  Variable  used  to  adjust  acceleration  in 
Phase  3, constant  q  path. 

ALPHA  Angle  of  attack  .cegrees 

ALPHA2MAX  Maximum  angle  of  attack  in  Phase  2, degrees 

ALPHA5MAX  Maximum  angle  of  attack  in  Phase  5, degree 

ALPHACRUISE  Angle  of  attack  for  maximum  L/D  cruise 
in  Phase  4 

ALPHAMAX  Maximum  angle  of  attack  in  Phase  I  thru 
4  .degrees. 

ALPHAMIN  Minimum  angle  of  attack  .degrees 

ANG  57.3  degrees/radian 

AREF  Vehicle  lift  and  drag  reference  area, sq.ft. 

B1  Adjustment  "knob"  for  controlling  Phase  3 

constant  q  path .Presently  set  to  0.009 

BIGG  Universal  gravitational  constant 

C0  Velocity  of  sound, ft. /sec. 

Cl  Adjustment  "knob"  for  controlling  Phase  3 

constant  q  path.  Presently  set  to  0.1 

C2  Adjustment  "knob"  for  controlling  Phase  3 

constant  q  path .Presently  set  to  2.0 

CD  Vehicle  drag  coefficient 

CDB  Zero  lift  drag  coef f ic ient . Val ues  are 

input  as  table  3  of  input  data  file. 

CL  Lift  coefficient 

CLALPHA  Lift  coef f i c ient /degree . Val ues  are 

input  as  table  5  of  input  data  file. 

COUNT  Loop  counter  for  controlling  screen  and 

printer  interval. 

CPCV  Air  specific  heat  ratlo.Not  used 

01  Adjustment  “knob"  for  controlling  Phase  3 


Internal 

Input 

Input 

Input 

Input 

Output 

Output 

Input 

Input 

Internal 

Output 

Input 

Input 

Internal 

Input 

Input 

Internal 

Input 

Input 

Internal 

Internal 

Input 

Input 

Internal 

Input 

Internal 

Input 

Internal 

Internal 

Input 
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constant  q  path. Presently  set  to  0.1 


DELCO  Friction  increment  of  drag. Values  are  Input 

input  as  table  4  of  input  data  file. 

OELPRINT  Print  and  screen  output  interval  Input 

OELQ  Adjustment  “knob”  for  controlling  Phase  3  Input 

constant  q  path. A  value  equal  to  5%  of 
Q0C0MD  is  suggested. 

OELTHRUST  Iteration  tolerance  on  commanded  thrust  Internal 

DRAG  Vehicle  drag.lbf.  Output 

0T  Trajectory  integration  interval , sec.  Input 

DV  Reduced  integration  interval  used  near  Internal 

end  of  flight, secs. 

DXI  Variable  in  table  interpolation  routine.  Internal 

DYJ  Variable  in  table  interpolation  routine.  Internal 

FAMAX  Maximum  fuel/air  ratio  Internal 

FASTOIC  Stoichiometric  fuel/air  ratio  Input 

FC  Vehicle  centrifugal  force, lbf.  Internal 

FINISH  Clock  time  at  end  of  problem.hr/min/sec  Output 

FLAG5  Selector  for  Phase  S.ballistic  ascent.  Input 

Enter  t  if  not  used.l  If  ballistic  phase 
is  used. 

FLAGCRUISE  If  vehicle  is  to  cruise  at  maximum  L/D,  Internal 
this  variable  is  reset  from  zero  to  a 
value  of  1  to  bypass  optimum  alpha  search 
once  it  is  made. 

FREQ  Frequency  of  output  to  plotter  file  Input 

at  prompt 

FREQC  Plotter  file  counter  variable  Internal 

G!  Freestream  specific  airf low,lbm/(sec-sq  in)  Internal 

GA  Gravitational  acceleration, ft/sec-sec  Internal 

GAMMA  Vehicle  flight  path  angle.degrees  Internal 

GAMMA!  "Old”  flight  path  angle  at  start  of  Internal 

integration  interval .degrees 

GAMMA4MAX  Maximum  flight  path  angle  permitted  Input 

during  Phase  4  maximum  L/D  cruise  .degrees . 

A  small  value, less  than  !.5  degrees  is 
suggested . 

GAMMA5  Flight  path  angle  to  be  achieved  before  Input 

ballistic  path  begins, degrees 

GAMMAMAX  Maximum  allowed  flight  path  angle  for  all  Input 

phases  except  4  and  5, degrees. 

GAMMAOLO  Same  as  GAMMA! .degrees  Internal 

GC  32.174  lbm  ft  per  sec  sq/lbf  Internal 

H  Geopotential  altitude, ft.  Internal 

HB  Vehicle  flight  a  1 1 i tude  ,f t . An  initial  Input  and 

value  is  input. Values  during  run  are  Output 

output . 

HI!  "Old”  value  of  flight  altitude  at  start  Internal 

of  integration  interval, ft. 


75 


Phase  2  attitude  to  be  achieved  before  VB2  Input 
can  be  exceeded, ft. 


HI  2 


HIOLO 

Same  as  H0B 

Internal 

HOEl 

Altitude  variable  used  in  atmosphere 
routine , ft . 

Internal 

I 

Integration  loop  counter  and  index  variable 

Internal 

11 

Counter  and  index  variable 

Internal 

12 

Counter  and  index  variable 

Internal 

I  FILE 

Input  file  designator 

Internal 

IL 

Input  file  designator 

Internal 

ILTI 

Index  variable  in  table  loader 

Internal 

ISP 

Specific  impulse, secs. Sum  of  airbreathing 
and  rocket  thrust  divided  by  sum  of  air- 
breathing  fuel  and  rocket  propellant  flows. 

Output 

ISPA 

Specific  impulse  of  airbreather.sec. 

Values  are  input  as  table  1  of  input 
data  file. 

Input  and 
Output 

ISPAVG 

Running  average  specific  impulse  since 
start  of  flight, secs. 

Output 

ISPEFF 

Instantaneous  effective  specific  impulse 
defineu  by  (weight/gc)X(delta  vel oc ity/delta 
weight. 

Output 

ISPMEAN 

Running  mean  value  of  ISPEFF  since  start  of 
f  1  ight . 

Output 

I SPR 

Rocket  specific  impulse, secs. Values  are 
input  as  table  7  of  input  data  file. 

Input  and 
Output 

ITAB 

Input  table  number. Eight  tables  are  current¬ 
ly  used  in  input  data  file. 

Input 

ITCUM 

Cumulative  total  effective  impulse  since 
start  of  fl  ight , lb . secs . 

Internal 

I  TO 

Input  data  table  variable. 

Internal 

ITSUM 

Cumulative  total  propulsive  Impulse  since 
start  of  f  I  ight ,  lb . secs . 

Internal 

IX 

Variable  used  in  table  loader  and  lookup. 

Internal 

IY 

Variable  used  in  table  loader  and  lookup. 

Interna  1 

IZ 

Variable  used  in  table  loader  and  lookup. 

Internal 

J 

Variable  used  in  table  loader  and  lookup. 

Internal 

J1 

Variable  used  In  table  loader  and  lookup. 

Internal 

J2 

Variable  used  in  table  loader  and  lookup. 

Internal 

K 

Ratio  of  specific  heats. Also  a  variable 
used  in  table  loader  and  lookup. 

Internal 

Kll 

Variable  used  in  table  interpolation. 

Internal 

K12 

Variable  used  in  table  Interpolation. 

Internal 

K21 

Variable  used  in  table  interpolation. 

Internal 

*22 

Variable  used  in  table  Interpolation 

Internal 

LIFT 

Vehicle  lift.lbf. 
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Output 

LOADFAC 

Maximum  vertical  acceleration  allowed 
during  Phase  2  climb.g's. 

Input 

LOVRO 

L/D  value  during  Phase  4  cruise  alpha 
search 

Internal 

LOVROOtD 

"Old1'  L/D  value  used  in  Phase  4  cruise 
alpha  search 

Internal 

M 

Variable  used  in  table  interpolation 

Internal 

MB 

Flight  Mach  number . Input  a  small  starting 
value, say  .001. 

Input  and 
Output 

MB  6 

“Old"  value  of  Mach  number 

Internal 

MEARTH 

Mass  of  the  earth 

Internal 

N 

Integration  counter  for  equations  of  motion 

Internal 

N1 

Variable  in  atmosphere  routine 

Internal 

NSTAGES 

Number  of  vehicle  stages.l  or  2. 

Input  at 
prompt 

PB 

Atmospheric  static  pressure, psia. 

Internal 

PHASE 

Flight  phase  designator.Values  are  1  thru  5. 

Internal 

PHI 

Airbreathing  fuel  equivalence  ratio. 

Output 

PHIMAX 

Maximum  airbreathing  fuel  equivalence 
ratio. Values  are  input  as  table  2  of 
input  data  file. 

Input 

PRCOM 

Command  for  printer  output  during  run. 

It  is  set  to  1  by  responding  "V  at  prompt 
for  printer  output  decision. 

Interna  1 

PRCOUNT 

Counter  for  desnignating  number  of  output 
points  to  be  printed  per  page. 

Internal 

Qfl 

Flight  dynamic  pressure  ,  lbf/ sq . ft . 

Output 

00  COM 

Calculated  dynamic  pressure  command 

Internal 

QBCOMD 

Commanded  dynamic  pressure  at  start  of 

Phase  3,!bf/sq.ft. 

Input 

Q0 FINAL 

Desired  fimal  dynamic  pressure  at  orbital 
or  final  cond  it  ion  ,  lbf /sq . ft . 

Input 

Q0MAX 

Maximum  dynamic  pressure  allowed  during 
Phase  3,  lbf/sq. ft. 

Input 

Q0MIN 

Minimum  dynamic  pressure  allowed  during 
Phase  3, lbf/sq. ft. 

Input 

Q0OLO 

"Old"  value  of  dynamic  pressure, Ipf/sq, ft. 

Internal 

Q0L01 

Same  as  Q80LD. 

Internal 

RB 

Density  variable  in  atmosphere  routine. 

Internal 

RANGE 

Flight  range, nautical  miles. 

Output 

RANGES 

"Old"  value  of  range, nm. 

Internal 

RATEOELO 

Commanded  rate  of  change  of  dynamic 
pressure  in  Phase  3,psf/sec. 

Internal 

RATEOQ 

Current  rate  of  change  of  dynamic 
pressure  ,psf /sec . 

Internal 

RE 

Earth  radius 

Internal 

REARTH 


Same  as  RE 


Internal 


RE  S 1 


Not  used. 


RH8  Atmospheric  density , Ibm./cu . ft .  Internal 

RH8B  density  variable  used  in  atmosphere  routine  Internal 

RH801D  "Old"  value  of  density  ,lbm./cu. ft.  Internal 

STAGE  Identifier  in  data  table  for  first  stage  Input 

or  second  stage. Enter  1  or  2  in  input  data 
file. 

START  Clock  time  at  start  of  problem, hr/min/sec  Output 

SWITCH4  Used  for  Phase  4  cruise. Zero  selects  con-  Input 

slant  altitude  cruise.l  selects  constant  ' 

L/D  cruise. 

T0  Static  temperature  of  atmosphere  ,degrees-R.  Internal  * 


THRUST  Vehicle  thrust.lbf.  Output 

THRUSTAMAX  Maximum  available  airbreathing  thrust.lbf.  Internal 

THRUSTCOM  Commanded  thrust.lbf.  Output 

THRUSTMAX  Maximum  total  thrust  available, Ibf.  Output 

THRUSTR  Rocket  thrust.lbf.  Internal 

THRUSTRMAX  Maximum  rocket  thrust  ava i lable , Ibf .  Internal 

TIMEX  Might  time  since  takeoff  .sees .  Output  and 

Input  small  initial  value, say .881 .  Input 

TIMEX0  “Old"  value  of  time, secs.  Internal 

TIME X2  Flight  time  at  start  of  2nd  stage, secs.  Internal 

TQ0  Time  variable  for  testing  if  Phase  2  is  Internal 

completed. 

TSET5  Time  variable  in  Phase  5  to  indicate  that  Internal 


commanded  flight  p»fh  angle  has  been 
reached. When  value  becomes  1, ballistic 
ascent  can  begin. 

TTAKEOFF  Time  variable  used  to  test  if  takeoff  is  Internal 
completed. 


V  Identifier  for  printer  output  variables  Internal 

V6  Vehicle  flight  ve loc ity  ,  f t . /sec .  Output 

V08  "Old"  value  of  flight  velocity  , ft. /sec.  Internal 

VB2  Maximum  velocity  permitted  in  Phase  2  Input 

until  H02  altitude  is  reached ,ft . /sec . 

V85  Velocity  at  which  Phase  5  pullup  to  Input 

ballistic  ascent  is  to  begin, ft. /sec. 

V8CRUISE  Commanded  Phase  4  cruise  veloc i ty  ,ft . /sec .  Input 


V6TEMP  Velocity  at  which  commanded  dynamic  pressure  Input 
is  made  a  function  of  veloc  ity, ft. /sec . 

This  variable  is  used  with  Q8FINAL  to 
approximate  a  temperature  limited 
ascent  path. 

VAL1  Variable  in  table  i n t erpol at i on .  Internal 

VAL2  Variable  in  table  interpolat  ion.  Internal 

VAL I  variable  returned  from  table  interpolation  Internal 

VARNUM  Number  of  output  variables  written  to  Internal 

plotter  file. 
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VFINAL  Commanded  final  flight  velocity, ft. /sec.  Input 
For  a  single  stage  going  to  orbit, enter 
a  value  beyond  orbital  velocity, say  27060. 

For  a  single  stage  suborbital  vehicle, 
enter  the  desired  final  velocity. 

For  a  2-stage  veh ic le .enter  the  desired 
staging  velocity  in  the  1st  stage  data 
table  and  enter  27000  or  a  desired  sub- 
orbital  velocity  in  the  2nd  stage  data  table. 


VORBIT  Orbital  velocity , ft. /sec.  Output 

VTAKEOFF  Commanded  takeoff  velocity, ft. /sec.  Input 

WDOT  Total  airbreathing  fuel  and  rocket  prop-  Internal 

ellant  flow  rate, lbm. /sec. 

WOOTA  Propulsive  air  flow, lbm. /sec.  Output 

WOOTF  Airbreather  fuel  flow, lbm. /sec.  Output 

WOOTFMAX  Maximum  airbreather  fuel  flow, lbm. /sec .  Internal 

WOOTP  Rocket  propellant  flow  rate, lbm. /sec.  Output 


WDOTPMAX  Maximum  rocket  propellant  flow  rate, lbm. /sec. Input 


Values  are  input  as  table  8  of  input  data 
file. 

WEIGHT  Vehicle  weight, lbm.  Output 

WEIGHT0  "Old"  value  of  vehicle  weight, lbm.  Internal 

WFINAL  Vehicle  final  weight, lbm.  Input 

WFUEL  Vehicle  fuel+propel lant  weight, lbm.  Input 

Value  not  requ i red .enter  zero  if  desired. 

WLAUNCH  Starting  weight  of  vehicle, lbm.  Input 

XARG  Table  lookup  variable  Internal 

YARG  Table  lookup  variable  Internal 


Simple 

String  Variables 

At 

Phasel  .takeoff  rol 1 

Output 

B$ 

Phase  2, climb  at  constant  L/W 

Output 

Ct 

Phase  3, climb  at  commanded  Q 

Output 

CMDlt 

Command  string  to  print  1st  stage  input 
file  to  printer. 

Internal 

CM02t 

Command  string  to  print  2nd  stage  input 
file  to  printer. 

Internal 

Ot 

Phase  4, cruise  at  V0cruise 

Output 

Et 

Phase  5,pullup  and  ballistic  ascent 

Output 

I  FN$ 

Input  data  file  name  (single  stage  or  two 
stage)  read  from  screen  at  prompt. 

Input  at 
prompt 

I  FN2$ 

Input  data  file  name  for  2nd  stage. Name 
is  read  from  screen  at  prompt. 

Input  at 
prompt 

PFNt 

Plotter  file  name  entered  at  prompt 

Input  at 
prompt 

PLOTt 

Set  to  “ Y“  or  "N"  in  response  to  prompt 
for  decision  to  create  plotter  file. 

Input  at 
prompt 

PRCOMt 

Set  to  “Y”  or  "N"  in  response  to  prompt 
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Input  at 

for  decision  for  printer  output  during  run.  prompt 


0$ 

Used  to  halt/resume  program 

Internal 

STARTTIMES  Clock  time  at  start  of  program.hr/min/sec 

Output 

TITL1$ 

Text  of  first  line  of  input  data  file. 

Input 

TITL2$ 

Text  of  second  line  of  input  data  file. 

Input 

WTIN$ 

Set  to  "Y"  or  "N"  in  response  to  prompt 
for  listing  of  tabular  data  in  input  file. 

Input  at 
prompt 

X* 

Variable  which  is  set  to  A$  thru  E$  as 
each  flight  Phased  thru  5)  is  reached. 

Internal 

Array 

Real  Variables 

AT 

Sets  number  of  points  in  tables. 

Internal 

DELGAM 

Change  In  flight  path  in  one  integration 
pass,  degrees. 

Internal 

DELH 

Change  in  altitude  in  one  integration 
pass, feet. 

Internal 

OELR 

Change  in  range  in  one  integration  pass, 
feet. 

Internal 

DELV 

Change  in  velocity  in  one  integration 
pass, ft. /sec. 

Internal 

OELU 

Change  in  weight  In  one  Integration  pass, 
lbm. 

Internal 

GA 

Not  used  as  array  variable. 

Internal 

GLMB 

Variable  used  In  atmosphere  routine. 

Internal 

HB 

Variable  used  in  atmosphere  routine. 

Internal 

1 1ST 

Variable  used  In  table  loader. 

Internal 

10 

Variable  used  in  table  loader. 

Internal 

J1ST 

Variable  used  in  table  loader. 

Internal 

RHBB 

Air  density  variable  used  in  atmosphere 
routine. 

Internal 

TG 

Not  used. 

Internal 

TMB 

Temperature  variable  used  in  atmosphere 
routine. 

Internal 

Array 

String  Variables 

T$ 

Table  title  used  in  table  loader. 

Internal 

U. 5. Government  Printing  Office  19S9  -  '48  05b  24  1 48 


80 


